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HIGHLIGHT  SUMMARY 

Compaction,  the  mechanical  densif ication  of  soil,  is  an 
effective  means  of  improving  the  behavior  properties  of  soil  in  the 
field.   While  it  improves  the  overall  integrity  of  the  soil  mass,  com- 
paction can  also  induce  the  tendency  to  swell  in  many  fine  grained 
soils.   This  swelling  can  result  in  a  lowered  soil  strength  or  damage 
to  structures  placed  on  or  in  the  fill. 

It  is  the  desire  of  the  geotechnical  engineer  to  quantitatively 
predict  and  control  the  overall  performance  of  field  compacted  soil. 
The  state-of-the-art  suggests  that  the  behavior  of  laboratory  compacted 
soil  can  be  extrapolated  to  field  soil  compacted  at  a  similar  water 
content  (relative  to  optimum  water  content)  and  to  a  similar  dry  densi- 
ty.  Data  indicate  that  such  extrapolations  are  often  poor.   As  an 
alternative  to  this  method,  test  embankments  can  be  constructed,  but 
these  are  usually  too  costly  for  small  scale  projects. 

This  study  investigates  the  swell  pressure  response  of  a  medium 
plastic  clay  (AASHTO  A-6)  compacted  both  in  the  field  and  in  the  labor- 
atory.  Field  compaction  was  achieved  using  a  Caterpillar  Model  825 
compactor  and  a  RayGo  Rascal  Model  420C  vibratory  roller;  Proctor  im- 
pact and  Hveem-type  kneading  procedures  were  employed  in  the  laboratory. 
Three  energy  levels  and  various  molding  water  contents  were  used  for 
each  compactor  type  to  study  their  effects  on  swell  pressure  behavior. 


Constant  volume  swell  tests  were  performed  in  two  benchtop  consolido- 
meters. 

Statistical  regression  models  were  found  for  dry  density  and 
swell  pressure  response  for  each  compactor  type.   In  general,  these 
relationships  were  influenced  by  water  content  deviation  from  optimum 
moisture  content,  compactive  effort,  dry  density  (swell  pressure 
models  only),  and  various  interactions  of  these  terms.   The  field  com- 
pactors produced  dissimilar  swell  pressure  models;  the  laboratory 
methods  produced  very  similar  models,  but  these  were  quite  different 
from  the  field  models. 

Variability  of  the  dry  density  and  swell  pressure  magnitudes 
was  found  to  be  sizeable  but  predictable.   Field  samples  exhibited 
larger  variabilities  than  laboratory  compacted  samples.   Larger  un- 
certainty in  the  independent  regression  variables  caused  an  increase 
in  the  dry  density  and  swell  pressure  variabilities. 

Prediction  models  were  designed  to  provide  expected  average 
and  maximum  values  for  field  swell  pressure  response.   The  expected 
maximum  swell  pressure  is  the  sum  of  the  expected  mean  and  the  expected 
variability  of  swell  pressure.   The  prediction  models  are  presented  in 
graphical  form  for  use  by  the  design  engineer.   These  charts  enable  the 
engineer  to  select  compaction  variables  that  will  ultimately  control 
the  field  swell  pressure  to  a  predetermined  maximum  value. 

The  prediction  models  are  also  presented  in  computer  tabulation 
format.   The  tables  contain  values  of  the  compaction  variables  and 
their  uncertainty,  the  expected  mean  swell  pressure,  and  the  variabil- 
ity associated  with  it.   This  format  is  more  suitable  for  quality 
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control  work  and  evaluation  of  swell  pressure  behavior  for  compaction 
projects  under  construction. 

A  user's  guide  was  developed  and  is  entitled  "Application  of 
Results."   By  use  of  an  example  and  commentary,  it  indicates  how  the 
results  of  this  study  can  be  used  in  both  engineering  design  and 
quality  control  work.   Also,  a  method  is  proposed  for  applying  the  re- 
sults of  this  study  to  A-6  soils  that  may  have  somewhat  different  com- 
paction characteristics. 


INTRODUCTION 

The  design  of  many  engineering  structures  incorporates  the  use 
of  compacted  soil.   Some  projects  utilize  compacted  fills  to  improve 
foundation  conditions;  others,  such  as  dams  and  highways,  rely  on 
earth  embankments  to  achieve  their  intended  purpose.   Proper  use  of 
compacted  soil  ensures  both  safe  and  economical  design  of  long-life 
structures. 

The  compaction  process  itself  mechanically  densities  soil 
by  expelling  air  from  the  voids  at  an  essentially  constant  water  con- 
tent.  The  most  readily  apparent  product  of  this  process  is  an  in- 
crease in  soil  density.   However,  this  increased  density  is  not  the 
sole  or  ultimate  goal  of  the  engineer.   Rather,  he  intends  to  improve 
the  in-service  behavior  of  the  soil.   This  improvement  may  require 
an  increase  in  the  soil  strength,  a  decrease  of  its  compressibility 
or,  possibly,  a  reduction  in  the  permeability  of  the  soil  mass. 

In  current  practice  the  compaction  specification  is  written 
in  terms  of  the  end  result,  i.e.,  the  required,  as-compacted  soil 
density  achieved  within  a  limited  range  of  water  content.   Characteris- 
tics of  the  field  compaction  equipment  are  also  sometimes  specified. 
In  most  cases,  these  specified  requirements  are  based  on  results  of 
standard  laboratory  compaction  tests;  these  tests  provide  the  basis 
for  quality  control.  Unfortunately,  the  state-of-the-art  practice 


does  not  specify  the  actual  behavior  properties  to  be  achieved.   In- 
stead the  engineer  must  infer  these  properties  through  proper  selection 
of  as-compacted  density,  water  content  range,  and  field  compactor 
type. 

The  foregoing  specification  controls  only  the  as-compacted 
nature  of  the  soil.   It  is  only  by  additional  inference  that  the  engi- 
neer includes  the  environmental  changes  which  occur  over  the  life 
of  the  structure.   Such  service-life  changes  include  the  consolidation 
of  the  soil  mass  under  the  weight  of  successive  lifts  as  well  as  addi- 
tional volume  changes  induced  by  an  increased  saturation  of  the  soil 
(i.e.,  swell  or  collapse). 

An  effort  is  being  made  to  better  predict  the  behavior  of 
the  field  compacted  soil.   This  endeavor  includes  the  testing  of  the 
same  soil  compacted  both  in  the  laboratory  and  in  the  field.   The 
samples  are  tested  for  strength,  compressibility,  and  volume  change 
behavior,  and  a  correlation  between  the  results  obtained  from  lab 
and  field  compacted  samples  is  then  produced  for  each  type  of  behavior. 
It  is  intended  to  quantitatively  predict  the  field  behavior,  with 
data  acquired  from  laboratory  compacted  samples,  by  using  these  corre- 
lations.  In  addition,  an  effort  is  made  to  enable  the  selection  of 
the  correct  state-of-the-art  compaction  specification  which  produces 
a  field  mass  exhibiting  the  desired  behavior  parameters. 

Several  studies  have  already  been  undertaken  on  a  silty  clay 
and  a  medium  plastic  clay.   Essigmann  (1976) ,  Scott  (1977) ,  and  Price 
(1978)  have  conducted  studies  of  soil  strength  for  the  silty  clay.  For 
the  more  plastic  clay  (St.  Croix  clay) ,  DiBernardo  (1979)  has  studied 


the  compressibility  in  both  the  as-compacted  state  and  after  satura- 
tion.  As-compacted  shear  behavior  has  been  investigated  by  Weitzel 
(1979) ,  while  Johnson  (1979)  has  examined  the  shear  behavior  after 
saturation  and  consolidation;  both  of  these  studies  were  conducted 
at  various  confining  pressures.   The  compressibility  and  as-compacted 
shear  behavior  for  field  compacted  samples  are  currently  under  investi- 
gation.  White  (1980)  has  recently  completed  a  study  of  the  effects 
of  field  and  laboratory  compaction  on  the  differences  produced  in 
the  soil  fabric,  as  suggested  by  pore-size  distribution  parameters. 

The  study  reported  herein  concerns  itself  with  the  soil's 
tendency  to  swell,  a  behavior  which  the  compaction  process  induces. 
This  phenomenon  has  been  described  by  Ladd  (1959)  and  it  definitely 
affects  the  long-term  behavior  of  the  compacted  mass.   The  tendency 
to  swell  is  measured  in  terms  of  the  swelling  pressure  that  develops 
while  the  soil  remains  at  a  constant  total  volume  in  the  presence 
of   free  water.   This  investigation  examines  both  laboratory  and  field 
compacted  material.   Statistical  correlation  techniques  relate  the 
tendency  to  swell  with  the  respective  compaction  procedures  and  vari- 
ables. 

The  ultimate  goal  of  this  investigation  is  to  predict  the 
field  behavior  of  the  compacted  mass,  i.e.,  its  tendency  to  swell. 
It  also  investigates  the  selection  of  compaction  specifications  which 
explicitly  control  the  swelling  phenomenon  to  tolerable  levels.  It 
is  hoped  that  the  engineer  can  use  these  results  to  quantitatively 
specify  a  compacted  soil  mass  which  exhibits  the  desired  behavior. 
For  specifications  which  have  not  been  prepared  in  this  manner,  results 


are  presented  that  can  be  used  for  quality  control  testing  and  for 
prediction  of  swell  pressures  induced  by  the  compaction  process. 


1  -  LITERATURE  REVIEW 

Mechanisms  of  Compaction  and  As-Compacted 
Soil  Structure 

Proctor  (1933)  gave  one  of  the  first  descriptions  of  the  com- 
paction process,  the  variables  which  affect  this  process,  and  the 
changes  in  soil  structure  that  result  from  it.   He  postulated  that 
densification  occurred  when  the  compacting  force  overcame  the  fric- 
tional  resistance  between  soil  particles.   His  concept  of  structure 
dealt  with  particle  spacing,  rather  than  degree  of  particle  orienta- 
tion.  In  his  theory  water  acted  as  a  lubricant  between  particles; 
a  change  in  water  content  resulted  in  a  different  particle  spacing 
for  a  given  method  of  compaction. 

Proctor  also  noted  the  influence  of  compactive  effort  on  parti- 
cle spacing.   He  stated  the  need  to  select  the  proper  type  of  construc- 
tion equipment,  as  well  as  an  appropriate  moisture  content,  to  ensure 
proper  performance  of  a  compacted  soil  mass.   Proctor  developed  pro- 
cedures for  standard  compaction  tests  and  the  penetration  resistance 
needle  that  are  still  used  in  engineering  practice  today. 

Physico-chemical  concepts  were  used  by  Lambe  (1960)  to  explain 
the  structure  and  engineering  behavior  of  compacted  clays.   Under 
this  theory,  individual  clay  particles  govern  the  behavior  of  the 
compacted  soil.   Because  of  their  minute  size,  clay  particles  are 


dominated  by  electrical  forces  rather  than  mass  forces.   The  negatively 
charged  clay  particles  attract  cations  and  polar  water  molecules  to 

them.   The  swarm  of  counterions  and  water  is  called  the  diffuse  double 

* 
layer,   and  the  clay  particle  plus  the  double  layer  is  termed  the 

micelle. 

In  his  theory  of  compaction,  Lambe  (1960)  introduced  the  con- 
cept of  water  deficiency.   This  idea  stated  that,  "Any  given  soil 
particle  under  any  given  state  of  stress  requires  a  certain  amount 
of  water  to  develop  fully  its  double  layer."   This  deficiency  gives 
rise  to  negative  pore  water  pressures  in  compacted  soil,  and  the  soil 
will  try  to  imbibe  enough  water  to  satisfy  its  deficiency. 

For  compaction  at  low  moisture  contents  the  double  layers 
are  depressed.   This  reduces  interparticle  repulsions  and  flocculation 
occurs.   Therefore,  a  flocculated  structure  with  a  low  degree  of  parti- 
cle orientation  characterizes  clay  soils  compacted  at  low  moisture 
contents.   As  the  molding  water  content  is  increased,  Lambe  postulates 
that  double  layers  expand  and  interparticle  repulsions  increase.   This 
causes  a  more  orderly  (parallel)  arrangement  of  particles  and  therefore 
a  higher  density.   Above  OMC  the  particles  achieve  an  even  more  orderly 
arrangement,  but  the  double  layers  are  large  enough  to  reduce  the 
concentration  of  soil  particles  per  unit  of  volume.    This  causes 
a  decrease  in  dry  density. 

An  increase  in  compactive  effort  results  in  higher  densities 
by  causing  a  more  parallel  arrangement  and  closer  spacing  of  the  clay 
particles.   At  high  molding  water  contents,  however,  the  effect  may 


Diffuse  double  layer,  diffuse  layer,  and  double  layer  are  used 
interchangeably  in  the  literature  and  also  in  this  report. 


be  to  merely  align  the  particles  without  significantly  changing  the 
particle  spacing. 

The  descriptions  and  hypotheses  given  above  are  based  on  the 
Gouy-Chapman  theory  for  electrical  potential  in  the  double  layer. 
Lambe  (1960)  noted  that  real  soils  do  not  conform  to  assumptions  made 
in  the  Gouy-Chapman  theory,  and  this  limits  its  application.   He  also 
pointed  out  that  the  variation  of  structure  with  molding  water  content 
depends  on  the  particular  soil.   Some  soils  may  exhibit  only  a  slight 
improvement  of  orientation  with  increasing  water  content,  while  others 
may  become  fully  dispersed  at  low  moisture.   Soils  with  intermediate 
responses  also  exist. 

Seed  and  Chan  (1961)  agreed  in  principle  to  Lambe' s  (1960) 
hypotheses  about  the  effects  of  a  soil's  structure  on  its  engineering 
behavior  and  to  his  theory  of  compaction.   After  conducting  their 
own  testing  programs,  however,  Seed  and  Chan  offered  some  additions 
and  modifications  to  Lambe 's  proposals.   For  compaction  at  low  moisture 
contents,  they  agreed  with  Lambe 's  hypothesis.   At  higher  water  con- 
tents, however,  Seed  and  Chan  stated  that,  "In  many  cases  it  is  the 
progressive  increase  in  shear  deformation  (for  a  constant  compaction 
effort)  as  the  water  content  is  increased  which  is  largely  responsible 
for  the  progressive  increase  in  the  degree  of  orientation  of  the  clay 

particles. " 

The  results  of  their  testing  program  led  to  the  following 

conclusions : 

1)   For  samples  prepared  dry  of  OMC,  all  methods  of  compaction 
produced  no  appreciable  shear  deformations,  and  essentially 
flocculated  structures  resulted  in  all  cases. 


2)  For  samples  of  the  same  composition  prepared  wet  of  OMC,  the 
method  of  compaction  had  a  pronounced  influence  on  the  result- 
ing structure. 

3)  Kneading  compaction  wet  of  OMC  produced  the  highest  shear 
strain  during  compaction  and,  therefore,  the  most  dispersed 
structure.   An  increase  in  the  foot  penetrations  resulted 
in  an  even  higher  degree  of  orientation. 

4)  Impact  compaction  wet  of  OMC  caused  slightly  less  shear  strain 
and  a  lower  degree  of  dispersion  than  kneading  compaction.   How- 
ever, these  differences  were  small. 

5)  Static  compaction  produced  little  shear  deformation  during  com- 
paction both  wet  and  dry  of  optimum.   A  relatively  flocculated 
structure  resulted  in  all  cases. 

6)  Vibratory  compaction  also  produced  small  shear  strains  during 
compaction.   The  structure  was  similar  to  that  for  static  com- 
paction, but  slightly  more  dispersed.   Apparently  the  vibrations 
enabled  the  particle  to  achieve  a  higher  degree  of  orientation. 

7)  The  magnitude  of  structural  changes  resulting  from  an  increase 
in  molding  water  content  was  dependent  on  the  soil  type. 
Lambe  (1960)  also  noted  this. 

An  effective  stress  theory  of  compaction  has  been  proposed  by 
Olson  (1963)  for  kneading  compaction.   This  theory  successfully  ex- 
plained the  densification  process  in  terms  of  changes  in  effective 
stresses,  pore  water  pressures,  and  shearing  strength  of  the  soil. 
Olson  believed  that  physico-chemical  variables  controlled  the  shape 
of  the  compaction  curve,  and  therefore  the  resulting  as-compacted 
soil  structure.   He  also  noted  reasons  why  impact  compaction  can  be 
used  well  above  OMC  but  kneading  methods  cannot;  this  latter  effect 
was  observed  in  this  study. 

Nalezny  and  Li  (1967),  Barden  and  Sides  (1970),  and  Hodek 
(1972)  offer  a  deformable  aggregate  model  to  explain  the  compaction 
process.   Due  to  electric  forces  between  individual  clay  particles, 


* 
the  particles  combine  to  form  aggregates.    These  authors  state  that 

the  behavior  of  these  aggregates  governs  the  soil's  response  to  the 

compaction  process.   The  particles  within  each  aggregate  may  assume 

a  flocculated  or  dispersed  structure,  however,  depending  on  water 

content,  type  of  soil,  type  of  compaction,  and  other  factors  previously 

mentioned. 

At  low  moisture  contents,  the  individual  aggregates  are  quite 
strong  and  do  not  readily  deform  into  the  available  pore  space  when 
the  compactive  effort  is  applied.   This  results  in  many  large  pores 
between  the  aggregates  (raacropores) .   As  the  water  content  is  increased, 
the  aggregates  become  softer  and  can  deform  more  easily  into  the  spaces 
between  them.   This  increases  the  density  and  the  macropores  become 
smaller.   Near  optimum  moisture  content  these  factors  produce  their 
maximum  effect;  maximum  dry  density  is  achieved,  the  macropores  are 
nearly  filled  with  water,  and  the  air  voids  become  discontinuous. 
Above  OMC  an  increase  in  water  content  enlarges  the  pore  volume  and 
dry  density  decreases. 

The  authors  postulate  that  the  structure  within  aggregates 
changes  from  a  flocculated  to  a  more  dispersed  structure  as  water 
content  is  increased.   Barden  and  Sides  (1970)  used  the  scanning  elec- 
tron microscope  to  study  this  concept  for  kaolin  and  a  naturally  occur- 
ring clay.   For  both  clays  at  low  magnifications  (110X) ,  individual 
peds  were  visible  at  moisture  contents  below  optimum  but  were  not 
observable  at  water  contents  above  OMC.   At  high  magnifications  (9500X) 


The  literature  also  refers  to  these  groups  of  particles  as 


clusters,  peds,  macropeds,  and  soil  crumbs. 
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both  clays  exhibited  an  essentially  dispersed  structure  both  wet  and 
dry  of  OMC;  the  natural  clay  appeared  to  be  somewhat  more  flocculated 
dry  of  OMC,  however. 

Young  and  Sheeran  (1973)  stated  that  "Examination  of  electron 
micrographs  and  soil  performance  shows  that  individual  particles  rarely 
act  as  single  particle  units.  .  .  .    The  different  sizes  and  arrange- 
ments of  particle  groups  observed  in  fabric  viewing  suggest  that  re- 
sponse behavior  might  be  controlled  by  the  kinds  and  arrangements  of 
these  particle  groups."   These  observations  lend  support  to  the  deform- 
able  aggregate  model  of  soil. 

Pore  size  distribution  measurements  of  laboratory  compacted 
clays  were  conducted  by  Diamond  (1970  and  1971),  Sridharan  et  al. 
(1971),  Ahmed  et  al.  (1974),  Bhasin  (1975),  Reed  et  al.  (1979),  and 
Garcia-Bengochea  et  al.  (1979)  .   These  tests  were  performed  on  kaolin, 
illite,  Boston  blue  clay,  and  other  types  of  clay.   Results  from  these 
studies  were  in  general  agreement  and  found  that  dry  of  OMC,  more 
pore  space  is  found  in  the  relatively  large  pores.   Bhasin  (1975) 
found  that  an  increase  in  compactive  effort  dry  of  OMC  decreased  the 
porosity  and  quantity  of  large  pores,  while  wet  of  OMC  it  had  little 
effect  on  the  porosity  or  pore  size  distribution.   These  results  are 
in  agreement  with  the  deformable  aggregate  soil  model. 

White  (1980)  conducted  pore  size  distribution  measurements 
on  the  St.  Croix  clay  used  in  this  study.   Laboratory  samples  were 
prepared  by  impact  and  kneading  compaction,  and  field  samples  were 
compacted  using  a  Caterpillar  Model  825  roller  and  a  RayGo  Rascal 
Model  420C  roller.   The  following  conclusions  were  made  from  his  study: 

1)   More  pore  space  is  found  in  the  relative  large  pores  dry  of 
OMC  for  all  compactors. 
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2)  Laboratory  fabric  and  field  fabric  are  different. 

3)  Impact  and  kneading  compaction  produce  the  same  fabric.   Ahmed 
et  al.  (1974)  also  found  this  to  be  true  for  Grundite  samples. 

4)  Differences  in  fabric  exist  on  the  dry  side  only  between  the 
Rascal  and  Caterpillar  field  compacted  samples. 

These  conclusions  are  in  general  agreement  with  the  deformable 
aggregate  model.   It  should  also  be  pointed  out  that  White  (1980)  con- 
ducted the  first  known  study  of  pore  size  distributions  in  field  com- 
pacted soil. 

From  the  literature  it  can  be  seen  that  every  compaction  vari- 
able can  have  an  effect  on  the  resulting  soil  structure.   The  micro- 
and  macrostructure  of  compacted  clay  appear  to  be  intimately  related 
to  each  other,  as  well  as  to  the  subsequent  behavior  of  the  compacted 
soil.   Factors  causing  changes  in  structure  should  ultimately  be  ob- 
served to  cause  changes  in  the  measured  behavior  response. 

The  molding  water  content  is  very  important  in  controlling 
the  as-compacted  structure.   Differences  in  structure  have  been  shown 
to  be  most  pronounced  when  dry-side  structures  are  compared  with  those 
produced  by  wet-side  compaction.   Therefore,  it  is  reasonable  to  expect 
different  swell  pressure  relationships  for  samples  compacted  dry  of 
OMC  when  compared  to  those  for  wet-side  samples. 

Seed  and  Chan  (1961)  found  that  different  methods  of  laboratory 
compaction  can  result  in  different  structures  and  behavior.   Based 
on  this,  one  can  expect  field  and  laboratory  compaction  to  likewise 
produce  different  structures  and  behavior.   White  (1980)  found  this 
to  be  true  for  soil  fabric. 
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Based  on  the  information  presented  in  the  literature,  a  com- 
prehensive study  of  field  swell  pressure  behavior  should  contain 
several   features.    The  range  in  water  content  should  produce  suffici- 
ent data  for  wet-side  and  dry-side  relationships  to  be  investigated. 
As  noted  by  Price  (1978)  variations  in  compaction  variables  and  behavior 
response  should  be  measured  and  accounted  for  in  the  analysis.   The 
effects  of  different  types  of  field  compaction  equipment  should  be 
investigated.    Field  relationships  should  be  compared  with  those  of 
laboratory   samples  to  determine  if  a  large  difference  in  their  as- 
compacted  structures  exists . 

Mechanisms  of  Swelling 

Theories  of  swelling  usually  account  for  this  phenomenon  at 
the  particle  level.   For  a  saturated  clay,  Ladd  (1959)  proposed  that 
osmotic  pressures  between  clay  particles  are  largely  responsible  for 
clay  swelling.   Osmotic  pressures  exist  when  an  aqueous  solution  is 
separated  from  pure  water  (or  a  solution  of  different  concentration) 
by  a  semi-permeable  membrane.   The  water  tends  to  pass  through  the 
membrane  to  dilute  the  solution.   The  pressure  required  to  prevent 
this  flow  in  termed  the  "osmotic  pressure"  of  the  solution. 

For  clay  particles,  the  concentration  of  cations  in  the  double 
layer  is  higher  than  in  the  surrounding  free  water.   These  cations 
are  held  in  the  double  layer  by  the  negative  charge  on  the  clay  parti- 
cles.  This  negative  electric  field  acts  as  a  semi-permeable  membrane 
by  allowing  water  to  enter  the  double  layer  and  by  preventing  the 
ions  from  leaving  it. 
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Free  water,  if  present,  attempts  to  enter  the  double  layer 
in  order  to  neutralize  the  difference  in  ion  concentration.   If  water 
enters,  the  double  layer  expands  and  the  interparticle  spacing  in- 
creases; i.e.,  swelling  occurs.   An  effective  stress  equal  to  the 
osmotic  pressure  is  required  to  prevent  this  swelling,  and  is  termed 
the  swelling  pressure.   The  osmotic  pressure  can  be  calculated  from 
the  van't  Hoff  equation 

P   =  RT(C  -  C  )2 
o         CO 

where  P  =  osmotic  pressure,  R  =  gas  constant,  T  =  absolute  tempera- 
ture, C  =  ion  concentration  in  the  free  pore  water,  and  C  =  ion 
o  c 

concentration  in  the  central  plane  between  particles.   Bolt  (1956) 
used  this  equation  and  the  Gouy-Chapman  theory  of  the  electric  double 
layer  to  calculate  interparticle  spacings  for  illite  and  montmorillonite 
suspensions.   Ladd  (1959)  cautioned  that  real  soils  are  far  from  the 
ideal  system  used  in  these  theories,  and  the  theories  should  be  used 
only  in  a  qualitative  manner  for  real  soils.   He  noted  that  the  influ- 
ence of  the  negative  electric  field  and  van  der  Waals  forces  on  the 
double  layer  water  could  also  be  partially  responsible  for  swelling 
in  clay-water  systems. 

Seed,  Mitchell,  and  Chan  (1961)  observed  that  mechanical  com- 
ponents of  swelling  can  also  exist.   While  the  elastic  compression 
of  solid  clay  particles  is  probably  insignificant  within  the  practical 
range  of  loadings,  volume  changes  resulting  from  the  bending  of  platy 
clay  particles  may  be  appreciable.   They  stated  that,  "By  making  rea- 
sonable assumptions  concerning  moduli  of  elasticity,  it  may  be  shown 
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that  a  typical  clay  plate,  acting  as  a  simple  beam  between  other  par- 
ticles, may  be  held  by  menisci  in  a  deformed  state  where  the  deflection 
is  about  10  percent  of  the  span." 

To  prove  the  existence  of  these  effects,  they  compacted  several 
identical  samples  at  a  high  water  content  and  allowed  them  to  swell 
in  calcium  acetate  solutions  of  various  concentrations.   At  high 
electrolyte  concentrations  a  constant  value  of  about  1.1%  swell  occurred. 
This  swelling  was  attributed  to  mechanical  factors  since  osmotic  pres- 
sures and  water  adsorptive  forces  were  negligible  for  these  test  condi- 
tions. 

Ladd  (1959)  and  Seed  et  al.  (1961)  suggested  that  air  can 
also  play  a  role  in  the  swelling  behavior  of  compacted  clays.   If 
the  air  is  initially  at  atmospheric  pressure  and  there  are  many  inter- 
connected tubular  air  voids,  water  entering  the  soil  from  all  directions 
will  cause  the  air  to  compress  and  exert  a  pressure  along  the  walls 
of  the  voids.   If  the  soil  structure  cannot  withstand  this  pressure 
an  expansion  may  occur.   These  conditions  can  exist  for  clays  compacted 
dry  of  OMC.   Terzaghi  and  Peck  (1967)  explained  that  this  increase 
in  air  pressure  is  at  least  partially  responsible  for  slaking  of  dry 
soils  immersed  in  water. 

From  this  discussion  it  appears  that  the  physico-chemical 
model  provides  a  better  explanation  of  mechanisms  of  swelling  than 
does  the  deformable  aggregate  model.   However,  the  interaction  between 
micro-  and  macrostructure  can  become  important  at  low  molding  water 
contents. 
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Compositional  Factors  Affecting  the  Swelling 
Behavior  of  Compacted  Clays 

Because  swelling  is  governed  primarily  by  physico-chemical 
forces,  many  factors  indeed  can  influence  the  magnitude  of  swell  or 
swell  pressure  induced  in  the  soil.   The  composition  of  the  clay  and 
its  adsorbed  cations  is  of  fundamental  importance.   Mitchell  (1976) 
noted  that  water  can  be  adsorbed  between  structural  sheets  in 
montmorillonite.   This  results  in  larger  amounts  of  swell  than  for 
minerals  in  which  this  does  not  occur  (e.g.,  illite).   He  also  noted 
that,  since  clay  minerals  govern  swelling  behavior,  an  increase  in 
clay  content  will  logically  increase  the  amount  of  swell. 

Lambe  (1960)  stated  that,  ".  .  .a  decrease  in  the  double 
layer  thickness  reduces  the  electrical  repulsion,  which,  in  turn, 
causes  a  tendency  toward  f locculation, "  and,  "...  that  a  tendency 
toward  flocculation  is  usually  caused  by  increasing  electrolyte  concen- 
tration,  ion  valence,  and  temperature;  and  decreasing  dielectric 
constant,  size  of  hydrated  ion,  pH,  and  anion  adsorption."   If  any 
of  these  variables  in  the  soil-water  system  changed,  as  would  occur 
in  the  presence  of  free  water,  Mitchell  (1976)  stated  that,  "In  gen- 
eral, the  thicker  the  double  layer  the  less  the  tendency  for  particles 
in  suspension  to  flocculate  and  the  higher  the  swelling  pressure  in 
cohesive  soils."   Since  a  flocculated  structure  results  in  more  parti- 
cle-to-particle contacts  and  smaller,  more  depressed  double  layers, 
it  is  logical  to  expect  a  higher  degree  of  swelling  from  this  type 
of  particle  arrangement.   That  is,  a  highly  depressed  double  layer 
can  undergo  a  larger  expansion  to  an  equilibrium  condition  in  the 
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presence  of  free  water  than  a  double  layer  that  was  initially  less 
depressed. 

Seed  et  al.  (1961)  and  Ladd  (1959)  demonstrated  that  an  increase 
in  the  pore  fluid  electrolyte  concentration  caused  a  decrease  in  swell- 
ing.  This  is  consistent  with  the  osmotic  pressure  and  double  layer 
theories  previously  discussed. 

In  addition  to  these  purely  physico-chemical  factors,  the 
compaction  variables  also  play  an  important  role  in  the  swelling  be- 
havior of  compacted  clays.   These  variables  include  the  molding  water 
content,  type  and  amount  of  compaction,  dry  density,  degree  of  satura- 
tion, and  the  resulting  as-compacted  soil  structure.   It  must  be  remem- 
bered, however,  that  these  variables  are  not  completely  separable 
from  physico-chemical  influences. 

Several  investigators  have  shown  that  molding  water  content 
has  a  dominant  influence  on  the  swelling  behavior  of  compacted  clays. 
These  investigators  include  Holtz  and  Gibbs  (1956),  Ladd  (1959),  Seed 
and  Chan  (1961),  Seed  et  al.  (1961),  Parcher  and  Liu  (1965),  and  others. 
This  should  be  expected  since  it  has  previously  been  shown  that  molding 
water  content  controls  the  double  layer  size,  as-compacted  soil  struc- 
ture, and  dry  density  for  a  given  method  of  compaction.   These  investi- 
gatiors  have  shown  that,  in  general,  a  decrease  in  molding  water  con- 
tent causes  an  increase  in  the  magnitude  of  swelling  or  swell  pressure 
induced  by  a  given  method  of  comapction. 

This  occurrence  can  be  attributed  to  the  following  factors. 
At  low  moisture  contents  the  double  layer  is  quite  depressed  and  a 
large  water  deficiency  exists  according  to  Lambe  (1960) .   In  the 
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presence  of  free  water  this  results  in  a  larger  intake  of  water,  larger 
expansion  of  the  double  layer,  and  larger  magnitude  of  swelling  or 
swell  pressure  than  for  samples  compacted  at  a  higher  water  content. 
Ladd  (1959)  noted  that  pore  water  tensions  exist  in  compacted  clay, 
particularly  if  compacted  dry  of  OMC.   As  water  enters  the  compacted 
soil  it  destroys  the  capillary  menisci,  thus  decreasing  the  effective 
stresses.   This  permits  the  double  layers  to  expand  until  a  new  equilib- 
rium condition  is  reached.   Over-all  swelling  of  the  sample  can  be 
prevented  by  applying  an  effective  stress  equal  to  the  swell  pressue. 

Molding  water  content  also  has  an  influence  on  the  as- 
compacted  soil  structure  which,  in  turn,  has  an  effect  on  the  swelling 
behavior  of  compacted  clays.   Seed  and  Chan  (1961)  concluded  that 
"Samples  compacted  dry  of  optimum  (tending  to  have  more  flocculated 
structures)  exhibited  greater  swell  pressures  than  samples  of  the  same 
final  composition  compacted  wet  of  optimum  (tending  to  have  more 
dispersed  structures)." 

Seed  and  Chan  (1961)  also  studied  the  effects  of  different 
methods  of  compaction  on  swelling  behavior.   For  samples  compacted 
to  the  same  composition  dry  of  OMC,  kneading  and  static  compaction 
produced  essentially  flocculated  structures  and  exhibited  similar  swell- 
ing characteristics.   Wet  of  OMC,  however,  statically  compacted  samples 
retained  an  essentially  flocculated  structure  and  produced  considerably 
higher  swell  pressures  than  did  the  relatively  dispersed  samples  pre- 
pared by  kneading  compaction. 

Holtz  and  Gibbs  (1956),  Seed  et  al.  (1961),  and  others  have 
shown  that  an  increase  in  dry  density  for  a  constant  water  content 
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results  in  larger  magnitudes  of  swelling  and  swell  pressure.   This, 
in  effect,  states  that  an  increase  in  compactive  effort  at  a  constant 
water  content  results  in  higher  tendencies  to  swell.   Changes  in  parti- 
cle spacings  and  soil  structure  are  responsible  for  the  increased  swell- 
ing. 

In  this  investigation,  compaction  variables  were  the  only 
compositional  factors  that  were  allowed  to  vary.   The  effects  of  purely 
physico-chemical  factors  was  reduced  substantially  by  maintaining  a 
nearly  constant  temperature  during  testing  and  by  using  only  demineral- 
ized  water.   In  this  way,  the  effects  of  the  compaction  variables  on 
swell  pressure  were  not  clouded  by  unrelated  factors. 

The  Influence  of  Test  Conditions 
on  Swelling  Behavior 

Several  methods  of  evaluating  the  swelling  characteristics 
of  expansive  soils  have  been  proposed  in  the  literature.   Holtz  and 
Gibbs  (1956)  showed  that  free-swell  tests,  colloid  content,  plasticity 
index,  and  shrinkage  limit  are  good  indicators  of  the  degree  of  expan- 
sion that  can  be  anticipated  from  different  soils.   Seed,  Woodward, 
and  Lundgren  (1962)  used  clay  content,  plasticity  index,  shrinkage 
limit,  and  activity  to  predict  the  swelling  potential  of  compacted 
clays. 

These  tests  work  well  for  their  intended  purpose;   that  is, 
to  identify  potentially  expansive  soils.   However,  these  methods  can- 
not measure  the  effects  of  different  compaction  conditions  on  the  in- 
duced tendency  to  swell.   Because  of  this  they  are  not  applicable  to 
the  intended  purposes  of  this  research. 
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In  recent  years,  several  investigators  have  studied  soil  suction 
as  a  means  of  predicting  the  swelling  behavior  of  clays.   Measured 
values  of  suction  must  be  correlated  with  meaningful  engineering  be- 
havior, such  as  swelling  or  swell  pressure,  before  they  can  be  used. 
These  correlations  can  be  quite  inaccurate  at  low  suctions  according 
to  Johnson  (1977).   For  these  reasons,  measurement  of  soil  suction 
was  not  employed  in  this  study. 

Tests  that  measure  swell  pressure  or  the  magnitude  of  swelling 
have  been  used  by  many  investigators.   Meaningful  information  about 
the  swelling  characteristics  of  a  soil  can  be  obtained  directly  from 
the  data  produced  by  these  tests.   Seed  et  al.  (1961),  Fredlund  (1969), 
and  others  have  shown  that  test  conditions  can  have  a  large  influence 
on  the  interpretation  of  data  obtained  from  these  tests.   The  following 
is  a  discussion  of  procedural  factors  that  affect  swell  pressure  and 
swelling  tests. 

Lambe  (1960)  demonstrated  the  effects  of  temperature  changes 
on  a  sample  under  constant  load.   An  increase  in  temperature  depressed 
the  double  layer  and  the  sample  compressed.   The  effects  were  reversed 
for  a  temperature  decrease.   Seed  et  al.  (1961)  showed  that  temperature 
changes  can  influence  the  components  of  the  measuring  system.   Tempera- 
ture variations  of  10°  to  15°  F  can  change  the  observed  swelling  by 
0.5  percent.   Therefore,  temperature  conditions  must  be  held  constant. 

Time  can  have  significant  effects  on  the  measured  value  of 
swell  response.   Seed  et  al.  (1961),  Fredlund  (1969),  and  many  others 
have  shown  that  several  days  may  be  required  for  the  full  swelling 
response  to  develop.   Sample  size  and  method  of  introducing  free  water 
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have  an  influence  on  the  time  required  for  testing.   By  using  a  rela- 
tively small  sample  size  and  allowing  water  to  enter  both  ends,  testing 
time  was  dramatically  reduced  for  this  study. 

Barber  (1956)  observed  that  as  the  time  between  compaction 
and  testing  increased,  the  measured  swell  under  loading  decreased. 
Nalezny  and  Li  (1967)  observed  that  samples  stored  for  three  days  ex- 
hibited less  swelling  and  lower  swell  pressures  than  samples  tested 
immediately  after  compaction.   They  attributed  this  to  thixotropic 
hardening  that  occurred  during  storage. 

Kassif  and  Baker  (1971)  studied  the  aging  effects  of  compacted 
clay  samples  stored  for  as  long  as  90  days.   They  noted  that  thixotropy 
in  a  compacted  clay  is  associated  with  an  increase  of  the  negative 
pore  water  pressures  (i.e.,  they  become  more  negative),  and  the  forma- 
tion with  time  of  electrochemical  bonds  between  clay  particles.   Their 
test  results  showed  that  swell  pressures  increased  for  short  storage 
durations  (less  than  about  10  days)  then  decreased  with  storage  time. 
The  effects  of  increased  negative  pore  water  pressures  dominated  the 
swelling  behavior  for  short  periods  of  storage.   As  storage  time  in- 
creased, more  bonds  formed  between  particles  and  swell  pressures  were 
lower.   To  minimize  these  time-dependent  effects,  storage  time  should 
be  the  same  for  all  samples  tested.   This  time  should  be  sufficient 
to  permit  elastic  rebound  and  water  equilibration  of  the  compacted 
sample. 

The  amount  of  expansion  measured  in  a  swelling  test  is  dependent 
on  the  surcharge  load  and  stress  history  during  the  test.   Gibbs  and 
Holtz  (1956)  showed  a  systematic  variation  between  swelling  and 
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surcharge  load.   As  would  be  expected,  lighter  surcharge  loads  permitted 
more  swelling,  but  the  relationship  between  swelling  and  surcharge 
was  not  linear.   Gibbs  and  Holtz  used  these  relationships  to  estimate 
actual  movements  that  might  occur  under  various  load  conditions  in 
the  field. 

Seed  et  al.  (1961)  showed  the  effects  of  varying  the  stress 
history  during  a  swelling  test.   Samples  placed  under  a  constant  1  psi 
surcharge  load  swelled  more  than  those  subjected  to  an  initially  higher 
load  that  was  subsequently  decreased  to  1  psi.   This  effect  can  be 
explained  in  terms  of  particle  reorientations  that  occur  in  the  sample. 
Under  low  surcharge  loads,  sample  expansion  can  take  place  with  ease. 
Higher  surcharge  pressures  cause  more  internal  particle  reorientations 
to  occur.   Lambe  (1960)  suggested  that  the  particles  tend  to  align 
with  their  flat  surfaces  perpendicular  to  the  direction  of  applied 
stress.   This  rearrangement  results  in  an  essentially  more  dispersed 
structure;  Seed  and  Chan  (1961)  have  shown  this  type  of  structure  to 
exhibit  lower  swelling  characteristics  than  the  initially  more  floccu- 
lated one.   Because  of  these  particle  reorientations,  samples  used 
in  constant  volume  swell  tests  were  not  subsequently  used  to  measure 
the  magnitude  of  swelling. 

Seed  er  al.  (1961)  have  demonstrated  that  a  volume  expansion 
of  as  little  as  0.1  percent  can  cause  a  large  error  in  the  magnitude 
of  the  observed  swell  pressure.   A  larger  amount  of  expansion  will 
obviously  result  in  a  lower  observed  swell  pressure.   This  phenomenon 
is  conversely  related  to  the  amount  of  swelling  that  will  occur  under 
different  static  surcharge  loads. 
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In  order  to  obtain  more  consistent  measurements  of  swell  pres- 
sures, constant  volume  swell  (CVS)  tests  should  be  conducted.   Several 
factors  must  be  considered  in  order  to  minimize  volume  changes.   Seed 
et  al.  (1961)  observed  that  the  faces  of  compacted  samples  developed 
a  curvature  after  being  trimmed.   They  attributed  this  to  elastic  re- 
bound.  When  these  samples  were  placed  in  the  testing  apparatus,  only 
a  small  portion  of  each  face  was  in  contact  with  the  loading  platens. 
This  permitted  a  volume  expansion  to  occur  before  an  entire  face  was 
in  contact  with  the  end  plates.   As  a  result  of  this  expansion,  measured 
swell  pressures  were  considerably  lower  for  these  samples  than  for  sam- 
ples trimmed  flush  immediately  before  testing. 

Most  CVS  tests  are  conducted  in  ordinary  consolidometers. 
Fredlund  (1969)  outlined  the  many  procedural  factors  that  can  affect 
the  results  of  such  tests.   Of  these  factors,  compressibility  of  the 
apparatus  can  result  in  the  largest  measurement  errors.   Fredlund  noted 
that  the  loading  mechanism,  porous  stones,  and  filter  paper  all  con- 
tribute to  the  apparatus  compressibility.   Filter  paper  can  contribute 
a  significant  proportion  of  this  compressibility;  it  not  only  has  an 
instantaneous  compression,  but  also  a  time-dependent  component.   He 
suggested  that  calibration  tests  be  conducted  by  replacing  the  soil 
specimen  with  a  steel  block.   Results  of  these  tests  can  then  be  used 
to  ensure  that  a  constant  volume  is  indeed  maintained  during  CVS  tests 
on  soil.   He  also  cautioned  that  machine  parts  (porous  stones,  loading 
platens,  etc.)  should  never  be  interchanged.   If  this  occurs,  calibra- 
tion tests  should  be  conducted  using  the  new  parts. 
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Seating  of  the  porous  stones  and  the  soil  sample  must  also 
be  accounted  for.   Fredlund  noted  that  most  seating  occurs  under  low 
pressures.   Therefore,  a  small  initial  seating  load  should  be  used 
in  order  to  get  a  representative  initial  deflection  reading.   He  also 
observed  that  friction  in  the  mechanical  components  of  the  apparatus 
is  negligible  for  most  conso Udometers. 

For  the  intended  purposes  of  his  study,  it  is  believed  that 
CVS  tests  provide  the  most  reliable  measurement  of  swelling  behavior 
and  also  the  most  meaningful  data.   The  procedures  set  forth  by 
Fredlund  (1969)  were  followed  in  this  investigation,  and  filter  paper 
was  replaced  with  thin,  polycarbonate  film. 

Escario  (1969) ,  GrOmko  (1969) ,  and  others  have  attempted  to 
measure  the  in  situ  swelling  behavior  of  expansive  clays.  However, 
the  soil  samples  compacted  in  the  field  and  collected  for  the  swell 
pressure  measurements  presented  in  this  report  are,  to  the  author's 
knowledge,  the  first  data  of  this  type  reported. 

As  can  be  seen  from  this  literature  review,  the  mechanisms 
controlling  the  compaction,  structure,  and  swelling  behavior  of  fine- 
grained soils  are  very  complex  and  not  completely  understood.   It  is 
hoped  that  future  research  will  result  in  clearer,  more  comprehensive, 
and  better  substantiated  theories  regarding  these  mechanisms. 
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2  -  APPARATUS  AND  EXPERIMENTAL  PROCEDURE 

Soil  Studied 
The  soil  used  in  both  the  field  and  laboratory  phases  of  this 
investigation  was  a  medium  plastic  clay  called  St.  Croix  clay.   It 
was  obtained  from  a  cut  and  fill  area  made  as  part  of  a  realignment 
project  of  State  Road  37,  approximately  four  miles  south  of  St.  Croix, 
Indiana.   It  was  a  shale  and  sandstone  residual  of  tan  color  with 
gray  and  red  mottling.   Other  pertinent  data  are  given  in  Table  2.1 
and  the  grain  size  distribution  curve  is  shown  in  Figure  2.1 

Test  Pad  Construction 

In  June  of  1978,  test  pads  were  constructed  in  the  fill  area 
of  the  previously  mentioned  SR  37  realignment  project.   The  purpose 
of  the  test  pad  was  to  create  field-compacted  soil  for  subsequent 
investigations  of  mechanical  properties  and  fabric  descriptors. 

Ten  test  pads  were  constructed  14  ft  (4.3  m)  wide  and  116  ft 
(35.4  m)  long.   Their  layout  is  shown  in  Figure  2.2.   A  subbase  of 
varying  thickness  (4  to  5  ft  or  more)  was  constructed  to  provide  a 
level  base  for  the  test  pads.   The  pads  were  placed  as  8  in  (20  cm) 
loose  lifts  and  a  tractor-drawn  disk  was  used  to  break  up  large  clumps 
of  soil. 

Five  test  pads  were  compacted  by  a  Caterpillar  Model  825  tamp- 
ing foot  roller  and  five  were  rolled  with  a  RayGo  Rascal  Model  420C 
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Table  2.1   Index  Properties  and  Classification  of  St.  Croix  Clay 


Atterberg  Limits  (%) 
Mean  (low,  high) 


Shrinkage  Limit,  w   (%) 

s 


Specific  Gravity,  G 

s 

Percent  Finer  Than  0.002  mm 
(clay  content) 


Field 

Laboratory 

WL 

40.1    (30.0,    53.2) 

38.3    (37.0, 

39.2) 

WP 

18.4    (16.7,    21.3) 

16.3    (15.6, 

16.7) 

h 

22.2    (16.7,    29.0) 

22.0    (20.9, 

22.7) 

11.8 
2.78 

34 


Percent  Finer  Than  0.001  mm 
(colloid  content) 
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Skempton's  Activity,  A 


0.65 


Unified  Soil  Classification 


CL 


AASHTO  Classification 


A-6 


Swell  Potential  High 

(Based  on  Gibbs  and  Holtz,  1956) 
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Figure  2.2   Test  Pad  Layout  Showing  Pad  Number,  Compactor  Type,  and 
Water  Content  Level 


Table  2.2   Field  Compaction  Equipment  Specification 


Caterpillar  Model  825 

Length  23  ft  -  4  in 

Width  13  ft  -  8  in 

Wheelbase  11  ft  -  8  in 

Weight  63,000  lb 


No.  of  Drums 

4 

No.  of  Pads/Drum 

65 

Drum  Length 

44.5  in 

Top  Speed  in  1st 

Gear 

3  mph 

RayGo  Rascal  Model  420C 

Length  18  ft  -  9  in 

Width  9  ft  -  0  in 

Wheelbase  9  ft  -  0  in 

Weight  25,160  lb 


No.  of  Drums  1 

No.  of  Pads /Drums  140 

Drum  Length  84  in 

Top  Speed  in  1st 

Gear  4  mph 


Vibration 
Frequency 
Dynamic  Force 


hydraulic,  direct  drive 
1100  -  1500  vpm 
32,000  lb 
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padded  vibratory  drum  compactor.    These  compactor  types  were  desig- 
nated "C"  and  "R"  respectively  for  sample  identification  purposes; 
their  specifications  are  given  in  Table  2.2.   These  rollers  were  selec- 
ted by  the  contractor  because  of  their  capability  to  perform  satis- 
factorily with  the  intended  soil. 

For  each  roller  type,  the  five  test  pads  were  brought  to  dif- 
ferent average  water  contents  prior  to  compaction.   The  range  in  water 
contents  from  pad  to  pad  was  selected  to  include  points  both  wet  and 
dry  of  what  was  believed  to  be  the  optimum  moisture  content  for  the 
respective  rollers.   For  identification  these  were  designated  "1" 
to   "5"  from  lowest  to  highest  moisture  level.   Within  each  pad  a 
uniform  moisture  distribution  was  desired.   Although  special  measures 
were  taken  in  the  field,  it  can  be  seen  from  Table  2.3  that  some  varia- 
tion did  result.   A  sizeable  variation  appears  to  be  inherent  in  compac- 
tion projects,  and  control  of  the  water  content  variation  has  been 
associated  with  proper  performance  of  the  embankment  (Price,  1978). 

The  desired  water  content  for  each  test  pad  was  achieved  in 
the  field  in  different  ways.  First,  moisture  contents  were  taken  at 
four  different  locations  in  the  loose  pad  using  the  Speedy  Moisture 
Tester  apparatus,  and  the  average  was  tabulated.   If  the  pad  was  in- 
itially too  wet,  the  tractor-drawn  disk  was  employed  to  facilitate 
drying  and  to  help  produce  moisture  uniformity.   Additional  water 
contents  were  taken  periodically  during  disking  until  the  desired 
magnitude  was  achieved.   If  the  pads  were  initially  too  dry  then  water 


* 
Unless  otherwise  noted,  future  references  to  "tamping  foot", 
"padded  drum",  and  "segmented  pad"  will  be  used  interchangeably  and 
will  apply  to  both  compactor  types. 
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Table  2.3  Density  and  Moisture  Content  Values  Obtained  from  Field 
Nuclear  Gage  Measurements 


No.  of  Passes 

4 

8 

16 

Pad  Compactor 
No.    Type 

Moisture 
Level 

Yd 

w 

Yd 

w 

Yd 

w 

108.2    22.5 
95.0    26.0 
101.7    23.1 
R5A 


106.9   23.6 

106.3  19.4 

103.4  19.2 
R5B 


111.9   17.0 

102.6  28.1 

105.7  21.8 
R5C 


97.2    26.1 
106.9   23.8 
100.0    20.0 
C5A 


101.2    25.1 

97.7    27.4 

103.8   21.1 

C5B 


106.0    22.6 

106.5    21.8 

106.3    20.7 

C5C 


103.1    21.2 
108.0    19.2 
100.0    22.1 
C4A 


106.6    20.5 

100.0    22.7 

105.5    19.4 

C4B 


111.6    17.0 
104.0    21.2 
99.4   21.4 
C4C 


105.0    10.3 
104.9    16.3 
107.4    12,7 
C3A 


109.1    14.7 
110.3    12.4 
103.8    41.0 
C3  B 


107.4  15.4 
105.9    15.6 

107.5  14.0 
C3C 


99.5  20.6 
99.4  21.1 
96.9   19.2 

R4A 


101.2    17.9 

102.5    25.7 

101.2    17.9 

R4B 


103.3    19.0 

95.6    22.8 

105.9    17.0 

R4C 


103.2  17.1 

105.3  15.7 
105.8    13.9 

C2A 


108.4  13.7 

111.5  16.1 
114.2    14.7 

C2B 


102.1    15.2 

108.4    17.2 

103.9   15.5 

C2C 


98.0    14.8 

104.0    14.9 

96.0    16.9 

CIA 


113.2  15.4 

109.3  14.5 
111.1    16.5 

C1B 


110.2    15.6 

112.5    16.0 

112.8    16.8 

C1C 


104.3    17.7 
98.0    16.6 
101.8    14.2 
R2A 


111.9   13.8 

106.8    12.5 

102.3    14.9 

R2B 


105.9   15.1 

110.7    14.2 

108.7    15.6 

R2C 
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Table  2.3   (continued) 


No.  of  Passes  4  8  16 

Pad  Compactor   Moisture 
No.    Type        Level 


10 


Yd 

w 

Yd 

w 

Yd 

w 

98.9 

15.7 

103.2 

16.8 

114.5 

13.6 

101.7 

16.4 

107.6 

12.9 

112.5 

14.7 

104.2 

14.6 

107.4 

14.4 

106.3 

14.5 

R1A 

RIB 

R1C 

103.5 

14.0 

107.6 

17.0 

111.9 

16.5 

98.3 

13.0 

104.8 

15.7 

105.8 

15.0 

95.9 

14.6 

106.4 

12.9 

108.6 

15.9 

R3A 

R3B 

R3C 

31 


was  added  by  passage  of  a  calibrated  water  truck;  the  soil  was  then 
mixed  again  by  disking.  Those  pads  which  were  at  the  desired  water 
content  were  briefly  disked  for  mixing  purposes  and  then  compacted. 

Once  the  desired  average  water  content  was  achieved  in  a  pad, 
soil  compaction  began  immediately.   Each  pad  was  sampled  after  4, 
8,  and  16  passes  of  the  compaction  eguipment  (see  section  titled  "Field 
Sampling"  for  discussion  of  sampling  techniques) .   These  samples  were 
labeled  "A",  "B",  and  "C",  respectively,  for  identification  of  energy 
levels.   Several  steps  were  taken  during  the  compaction  process  to 
reduce  variability  in  applied  energy.   Equipment  was  operated  at  top 
speed  in  first  gear  (about  3  mph)  for  all  passes.   Ample  room  was 
allowed  for  the  roller  to  reach  this  top  speed  before  entering  the 
actual  test  pad  area,  and  equipment  was  carefully  guided  to  prevent 
overlap  at  the  interface  of  adjacent  coverages.   The  same  equipment 
operator  was  used  for  all  compaction  work. 

Despite  these  precautions,  it  appeared  that  the  segmented 
pad  roller  itself  provided  a  major  source  of  input  energy  variability. 
During  a  single  pass  of  the  roller  drum  only  the  soil  under  the  raised 
pads  received  the  maximum  applied  energy;  areas  surrounding  the  pads 
might  have  had  no  compaction  at  all.   On  subsequent  passes  the  compac- 
tor pads  could  have  run  over  either  previously  compacted  soil  or  rela- 
tively uncompacted  areas.   As  the  number  of  passes  increased,  so  did 
the  overall  uniformity  for  the  entire  test  pad.   The  variation  in 
applied  energy  was  believed  to  be  one  cause  of  the  variability  in 
field  dry  density  as  shown  in  Table  2.3  and  the  resulting  compaction 
curves  presented  in  Figures  2.3  and  2.4. 
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The  compactors  themselves  could  have  provided  a  source  of 
variation;  their  configurations  differed,  as  did  the  ways  in  which 
they  applied  energy  to  the  soil.   In  a  single  pass,  the  Rascal  equip- 
ment applied  compactive  effort  through  one  wide  drum  followed  by  two 
rear  rubber  tires  spaced  3  ft  (0.9  m)  apart.   Higher  confinement, 
and  thus  higher  compaction,  might  have  been  achieved  near  the  center 
of  this  drum,  with  more  "squeezing-out"  having  occurred  near  its  ends. 
Applied  energy  was  divided  almost  equally  between  static  weight  and 
dynamic  vibratory  force;  the  effect  of  the  rubber  tires  was  unknown. 

The  Caterpillar  equipment  used  four  tamping  drums  to  apply 
compactive  effort.   For  each  pass  of  this  equipment,  two  drums  traveled 
over  a  given  area  of  soil  which  increased  the  probability  of  pad- 
to-soil  contact  for  a  given  soil  location;  the  Caterpillar  equipment 
also  had  approximately  twice  as  many  pads  as  the  Rascal  compactor. 
The  effect  of  confinement  was  probably  not  as  large  for  the  narrower 
Caterpillar  drums,  and  neglecting  travel  speed,  all  compaction  energy 
was  supplied  by  static  weight.   Actual  weight  distributions  for  the 
two  compactors  were  not  known.   In  view  of  the  above  discussion,  it 
would  be  suggested  that  the  Caterpillar  equipment  would  produce  more 
uniform  compaction,  but  this  suggestion  is  not  evident  in  the  data 
of  Table  2.3  and  Figures  2.3  and  2.4.   Price  (1978)  also  provided 
a  discussion  of  field  compaction  variability. 

Field  Sampling 
Sampling  of  the  test  pad  began  immediately  after  the  required 
number  of  passes  had  been  completed.   The  test  pad  was  laid  out  in 
a  2  ft  by  2  ft  (0.6  m  by  0.6  m)  grid  pattern  and  marked  around  the 
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perimeter  with  consecutively  numbered  stakes.   The  required  number 
of  sampling  tubes  were  lightly  lubricated  with  silicone  oil  and  were 
distributed  to  their  proper  locations  along  with  labels  and  plastic 
bags.   These  locations  had  been  determined  using  a  random  number 
process  in  an  earlier  phase  of  this  project.   A  typical  sample  label, 
and  its  explanation,  is  shown  in  Figure  2.5.   For  this  study,  six 
samples  were  taken  for  each  combination  of  compactor  type,  water  con- 
tent level,  and  energy  level.   A  total  of  840  samples  were  collected 
for  the  entire  project,  of  which  180  were  taken  for  this  work  on  swell- 
ing tendency. 

Sampling  tubes  were  driven  into  the  ground  with  a  drop  hammer, 
and  all  tubes  were  driven  at  the  bottoms  of  impressions  left  by  the 
segmented  pads  of  the  compactors;  this  was  done  in  an  attempt  to  re- 
duce the  variability  in  applied  energy  as  well  as  to  increase  the 
useable  sample  height.   After  being  driven,  the  tubes  were  dug  out 
of  the  ground,  wrapped  in  the  plastic  bags,  labeled,  and  carefully 
transported  to  the  extruding  area.   Extrusion  and  other  field  opera- 
tions are  to  be  explained  later  in  this  section. 

Sampling  tubes  and  drop  hammers  were  manufactured  in  the 
Central  Machine  Shop  at  Purdue  University.   The  swell-sample  tubes 
were  made  from  steel  tubing  with  an  internal  diameter  of  2.51  in 
(6.38  cm),  an  external  diameter  of  2.75  in  (6.99  cm),  and  were  5.0 
in  (12.7  cm)  long.   One  end  was  machined  to  form  a  cutting  edge.   The 
drop  hammer  consisted  of:   (1)  a  driving  head  with  a  recess  and  four 
set  screws  used  to  steady  the  sampling  tube;  (2)  a  pipe  threaded  into 
the  driving  head  to  guide  the  falling  weight;  (3)  a  disk  on  the  end 
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Figure  2.5   Typical  Field  Sample  Label  With  Explanation 
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of  the  pipe  to  use  as  a  handle  and  to  control  the  height  of  drop; 
(4)  the  falling  weight  itself;  and  (5)  two  wood  cushions,  one  between 
the  sampling  tube  and  driving  head,  and  one  between  the  driving  head 
and  falling  weight.   The  driving  assembly,  with  tube,  is  shown  in 
Figure  2.6.   The  16.7  lb  (7.6  kg)  weight  fell  28  in  (71  cm)  producing 
about  40  ft «lb  (54  N«m)  of  energy  per  blow,  with  several  blows  (usually 
10  to  25)  being  required  to  drive  each  tube.   The  complete  unit  weighed 
less  than  30  lb  (13.6  kg). 

In  place  density  and  water  content  measurements  were  made 
during  the  sampling  phase  by  ISHC  personnel.   The  Troxler  Model  3401 
nuclear  gage  was  used  at  various  locations  in  each  test  pad  for  these 
measurements.   Values  obtained  from  the  nuclear  gage  were  given  in 
Table  2.3  and  the  corresponding  field  compaction  curves  were  plotted 
in  Figures  2.3  and  2.4.   These  plots  indeed  showed  the  large  variabil- 
ity within  the  soil  mass.   Some  of  the  factors  which  contributed  to 
this  variability  have  previously  been  discussed.   Those  which  occurred 
during  the  sampling  phase  were:   (1)  the  weather  during  sampling  was 
warm  (90°F  +  daily  temperatures);  and  (2)  four  to  six  hours  elapsed 
from  start  to  finish  of  sampling  for  a  single  test  pad.   Although 
efforts  were  made  throughout  the  sampling  program  to  mitigate  moisture 
losses,  some  losses  probably  did  occur. 

Shortly  after  being  removed  from  the  ground,  samples  were 
extruded  from  the  tubes  using  either  an  hydraulic  jack  (see  Figure  2.7) 
or  an  electrically  driven  loading  press.   Samples  were  then  wrapped 
in  plastic,  covered  with  cheese-cloth,  and  heavily  coated  with  paraf- 
fin.  Next,  they  were  placed  into  plastic  bags  along  with  their 
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Figure  2.6   Driving  Assembly  with  Sampling  Tube 
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Figure  2.7   Hydraulic  Jack  Used  for  Extrusion  of  Tube  Samples 
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labels,  packed  into  5  gallon  metal  cans,  and  padded  with  shredded 
foam.   All  samples  of  similar  origin  were  packed  together,  and  the 
cans  were  then  transported  to  and  stored  in  the  Purdue  Geotechnical 
Laboratories. 

At  the  time  of  extrusion  several  observations  were  made  about 
the  samples.   Rock  fragments  appeared  in  all  samples  and  some  samples 
were  destroyed  because  large  rocks  or  gravel  pockets  were  present; 
other  samples  were  scored  by  rocks  being  dragged  along  their  periphery 
by  the  sampling  tubes.   Many  dry-side  samples  were  quite  brittle, 
broke  easily,  and  their  usable  height  was  generally  shorter  than  the 
average.   Some  samples  showed  horizontal  cracks  while  others  showed 
distinct   interfaces  between  individual  lumps  of  soil  which  had  not 
been  intermixed  by  the  compaction  process.   Samples  in  the  low  and 
intermediate  water  content  ranges  sometimes  showed  distinct  zones 
of  different  water  content.   This  was  probably  due  to  the  incomplete 
mixing  of  soil  lumps  which  were  at  different  initial  water  contents. 
High  water  contents  produced  the  longest,  most  uniform  samples. 

Preparation  of  Soil 

Soil  used  for  the  laboratory  compaction  phase  of  this  study 
had  been  obtained  directly  from  test  pads  six  and  eight.   Bag  samples 
had  been  transported  along  with  the  field  samples  and  stored  in  large 
metal  cans  at  the  Purdue  Geotechnical  Engineering  Laboratories.  While 
in  storage,  these  bag  samples  came  to  an  equilibrium  air-dry  water 
content  of  about  2.5%  (±  0.3%). 

In  this  air-dry  condition  the  soil  formed  very  hard  lumps, 
with  rock  fragments  of  varying  sizes  also  being  observed.   A  small 
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jaw  crusher  was  used  to  break  the  lumps  which  could  not  easily  be 
broken  by  hand.   In  the  preparation  process  the  bag  samples  were  first 
passed  through  the  No.  4  sieve,  and  rock  fragments  were  removed  from 
both  the  retained  material  as  well  as  from  that  which  passed  through 
the  sieve.   Soil  retained  on  the  sieve  was  then  transferred  to  the 
jaw  crusher.   After  crushing,  the  material  was  sieved  and  rock  pieces 
again  removed.   This  crushing  and  sieving  operation  was  repeated  for 
each  batch  of  soil  until  the  breakdown  was  accomplished.   After  siev- 
ing, the  soil  was  stored  in  two  large  metal  cans  which  were  lined 
with  polyethylene  bags. 

The  object  of  the  crushing  operation  was  to  produce  a  soil 
mixture  which  could  be  used  with  ASTM  Standard  Test  Designations  D  698 
(Standard  Proctor)  and  D  1557  (Modified  Proctor)  for  compaction  of 
soil  passing  the  No.  4  sieve.   To  avoid  crushing  rock  fragments  em- 
bedded within  soil  lumps,  the  jaws  on  the  crusher  had  been  set  wider 
than  the  No.  4  sieve  opening.   It  is  felt  that  the  crushing  process 
did  not  significantly  change  the  grain  size  distribution;  the  process 
prevented  crushing  of  most  larger  rocks  while  it  removed  only  a  small 
portion  of  those  fragments  passing  the  No.  4  sieve.   This  method  of 
crushing  proved  helpful  since  370  lb  (168  kg)  of  soil  was  prepared. 

Since  the  bulk  soil  remained  near  2.5%  water  content,  mixing 
proportions  were  easily  and  accurately  determined.   Soil  required 
for  compaction  (about  5.3  lb  or  2.4  kg)  was  placed  in  a  mixing  pan, 
and  demineralized  water  was  added  with  a  hand-operated  atomizer  until 
the  desired  water  content  was  reached.   Frequent  mixing  prevented 
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the  formation  of  large  aggregations  and  helped  to  produce  a  relatively 
uniform  moisture  distribution. 

After  mixing,  the  soil  was  put  into  two  layers  of  plastic 
bags,  closed  securely  with  twist  ties,  and  labeled.   The  bags  of  soil 
were  then  placed  in  a  humid  plastic  barrel  and  allowed  to  cure  a  min- 
imum of  3  days  prior  to  compaction.   Humidity  in  the  barrel  was  main- 
tained by  using  a  bottom  layer  of  water  and  gravel  and  a  tightly  fitted 
cover.   Previous  investigators  (see  Weitzel,  1979)  found  this  3-day 
period  to  be  appropriate  for  ensuring  moisture  homogeneity. 

Laboratory  Compaction 

Laboratory  compacted  samples  used  in  this  study  were  prepared 

by  both  kneading  and  impact  methods.   Kneading  compaction  was  used 

for  two  reasons:   (1)  companion  studies  in  this  project  used  the 

method;  and  (2)  it  was  believed  that  this  method  produced  shearing 

strains  and  loading  patterns  more  similar  to  those  produced  by  field 

compaction  (Weitzel,  1979).   Impact  compaction  was  used  to  provide   a 

standard  basis  for  comparison.    Laboratory  samples  were  labelled 
in  a  manner  similar  to  that  used  for  field  samples;  "K"  and  "I"  were 

used  to  designate  kneading  and  impact  samples,  respectively. 

Kneading  Compaction 
The  kneading  compactor  was  manufactured  by  the  August  Manu- 
facturing Company  of  Oakland,  California  and  is  shown  in  Figure  2.8. 
It  is  an  automatic  apparatus  with  mold  rotation  and  foot  tamping  being 
driven  by  an  electric  motor;  a  pneumatic-hydraulic  system  supplies 
the  foot  pressure.   The  mold  rotates  60°  between  tamps  that  are  applied 
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Figure  2.8   Kneading  Compactor 
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at  a  rate  of  30  per  minute,   with  full  face  coverage  of  the  soil  layer 
being  achieved  by  six  tamps  of  the  compactor  foot.   A  more  complete 
description  of  this  compactor  is  given  by  Gaudette  (1960) . 

Samples  were  compacted  at  three  energy  levels  (foot  pressures) 
chosen  to  simulate  the  compaction  curves  produced  by  impact  compaction. 
Figure  2.9  shows  the  relationship  between  kneading  and  impact  compac- 
tion obtained  by  DiBernardo  (1979).   This  figure,  along  with  prelimin- 
ary compaction  tests,  was  used  to  determine  the  appropriate  foot  pres- 
sures shown  in  Table  2.4. 


Table  2.4  Energy  Levels  for  Kneading  Compaction 


Energy 
Level 


Foot  Pressure 
(psi) 


Corresponding 
Impact  Test 


4 
7 

26 


Low  Energy  Proctor 
Standard  Proctor 
Modified  Proctor 


Energy  levels  were  labelled  "A",  "B",  and  "C"  for  identification  pur- 
poses similar  to  the  method  used  for  field  samples. 

To  compact  a  sample,  a  Standard  Proctor  split  mold  of  1/30 

3        3 
ft   (944  cm  )  with  a  collar  (see  Figure  2.10)  was  positioned  and 

securely  bolted  to  the  rotating  table  of  the  compactor.   The  mold 

and  collar  had  been  lightly  lubricated  with  silicone  oil  to  facilitate 
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Figure   2.9     Comparison  of   Impact   and  Kneading  Compaction  Curves    (after 
Di   Bernardo,    1979) 
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Figure  2.10   Standard  Proctor  Split  Mold  with  Collar 
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future  trimming  and  sampling  operations.   A  batch  of  soil  which  had 
been  allowed  to  cure  was  sampled  for  water  content  and  spooned  into 
the  mold.   The  soil  was  compacted  in  five  layers  of  approximately 
equal  thickness  and  each  layer  received  30  tamps  from  the  foot.   The 
top  of  each  layer  was  scarified  before  adding  more  soil  to  ensure 
sample  homogeneity.   This  was  most  important  for  the  drier  samples. 
Upon  comple  tion  of  compaction  the  mold  was  unbolted  from  the  com- 
pactor table,  the  collar  was  removed,  excess  soil  was  screeded  from 
the  top  of  the  sample,  and  mold  and  soil  were  weighed  for  the  density 
determination.   Compaction  curves  resulting  from  these  measurements 
are  shown  in  Figure  2.11. 

Impact  Compaction 
Impact  compaction  was  also  performed  for  three  different  energy 
levels.   The  Low  Energy  Proctor  samples  were  prepared  as  specified 
in  the  U.S.  Army  Engineering  Manual  EM  1110-2-1906  (1970)  and 
labelled  "A"  for  identification  purposes.   The  Standard  Proctor  (B) 
and  Modified  Proctor  (C)  were  specified  in  ASTM  D  698  and  D  1557, 
respectively.   Information  regarding  these  procedures  is  given  in 
Table  2.5.   As  with  kneading  compaction,  a  Standard  Proctor  split 
mold  was  also  used.   Procedures  for  soil  handling,  sampling  for  water 
content,  scarifying,  trimming,  and  weighing  were  the  same  as  for  knead- 
ing samples.   The  Impact  Compaction  curves  are  shown  in  Figure 
2.12. 
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Figure   2.12      Impact  Compaction  Curves    for  Compaction  Mold  Values 


50 


Table  2.5  Energy  Levels  for  Impact  Compaction 


Energy   Hammer    Drop      Number      Blows/ 
Level    Weight   Height    of  Layers    Layer 


Compaction 
Method 


A  5.5  lb  12  in 
B  5.5  lb  12  in 
C      10   lb    18  in 


3 

15 

3 

25 

5 

25 

Low  Energy  Proctor 
Standard  Proctor 
Modified  Proctor 


Tube  Sampling 

After  compaction  it  was  desired  to  treat  laboratory  samples 
in  a  manner  similar  to  that  used  for  field  samples.   For  this  reason, 
tube  sampling  was  performed  on  the  compacted  soil  in  the  mold. 

With  the  halves  of  the  split  mold  securely  tied  together  with 
heavy  latex  surgical  tubing  (designed  for  use  as  a  tourniquet) ,  the 
four  screws  holding  the  halves  together  were  unscrewed,  and  the  mold 
was  removed  from  its  base.   The  mold  was  then  placed  on  the  ram  of 
an  hydraulic  jack.  A  sampling  tube  (the  kind  used  for  field  sampling) 
was  lightly  lubricated  with  silicone  oil  and  centered  on  top  of  the 
soil  in  the  mold.   The  tube  was  then  jacked  in  at  a  constant  rate 
through  the  full  height  of  soil  in  the  mold.   Figure  2.13  shows  these 
features  with  the  sampling  tube  partially  inserted  into  the  soil. 

Excess  soil  was  removed  from  around  the  sampling  tube  and 
some  was  used  for  water  content  determinations.   The  sample  was  then 
extruded  from  the  tube  with  the  hydraulic  jack.   This  same  operation 
was  shown  for  field  samples  in  Figure  2.7.   After  extrusion,  each 
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Figure  2.13   Tube  Sampling  of  Laboratory  Compacted  Samples 
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sample  was  wrapped  in  a  double  thickness  of  polyethylene  bags,  cush- 
ioned with  paper  towels,  and  stored  in  a  humid  container.   Storage 
time  was  set  at  three  to  four  days  to  allow  equilibration  of  moisture 
within  the  sample. 

The  tube  sampling  procedures  were  the  result  of  several  trials 
at  obtaining  quality  samples.   During  the  sampling  operation,  the 
latex  tubing  provided  enough  lateral  support  to  prevent  the  soil  in 
the  mold  from  cracking  as  the  tube  was  inserted,  but  allowed  enough 
expansion  to  accommodate  the  added  volume  of  the  sampling  tube.   Pre- 
vious investigators  on  this  project  (Di  Bernardo,  1979;  Johnson,  1979; 
Weitzel,  1979)  used  a  somewhat  different  method  of  tube  sampling. 
They  compacted  soil  into  a  Standard  Proctor  mold  and  extruded  it  into 
a  mold  which  was  split  on  one  side.   Tube  sampling  was  then  done  in 
this  semi-split  mold.  Preliminary  work  done  for  this  study  indicated 
that  the  semi-split  mold  could  not  accommodate  the  volume  of  the  sam- 
pling tube  without  causing  excess  soil  disturbance. 

Trimming  of  Soil  Samples 
Constant  volume  swell  (CVS)  tests  were  conducted  in  Geonor 

bench  top  consolidometers.   Test  specimens  for  this  apparatus  were 

3         3 
cylindrical  with  a  2.44  in   (40.0  cm  )  volume,  a  1.99  in  (5.05  cm) 

diameter,  and  a  0.181   in  (2.0  cm)  height.   This  required  the  use  of 
a  special  Geonor  simulated  piston  sampling  apparatus. 

The  Geonor  sampling  apparatus  was  designed  for  use  with  soft 
clays;  its  purpose  was  to  cut  a  soil  sample,  of  the  correct  diameter 
and  with  minimum  disturbance,  directly  into  the  stainless  steel  confin- 
ing ring  used  in  the  consolidometer.   This  could  not  be  achieved 
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satisfactorily  with  the  relatively  stiff  clay  used  in  this  study.   Pre- 
liminary trials  resulted  in  test  specimens  which  were  loose  in  the 
confining  ring.   This  condition  was  undesirable  for  constant  volume 
swell  tests;  other  investigators  (Dawson,  1956;  Gromko,  1974)  have 
shown  that  allowing  even  a  very  small  amount  of  expansion  can  greatly 
reduce  the  measured  swell  pressure.   To  avoid  this  situation,  stainless 
steel  cutting  rings  were  manufactured  with  an  inside  diameter  0.039 
in  (1.0  mm)  larger  than  that  of  the  confining  ring.   These  cutting 
rings  produced  specimens  which  could  be  trimmed  to  fit  tightly  in 
the  confining  ring. 

Sample  number,  field  location,  and  initial  sample  height  were 
recorded  for  each  sample.   Other  notes  were  made  throughout  the  trim- 
ming process  that  indicated  such  things  as:   color,  homogeneity,  con- 
sistency, water  content,  plasticity,  defects,  irregularities  of  shape, 
and  the  number,  size,  and  location  of  rocks  and  voids. 

The  first  step  in  the  trimming  process  for  field  samples  was 
to  remove  the  wax  and  cheese  cloth  covering.   This  was  accomplished 
using  a  sharp  utility  knife  to  cut  through  the  covering,  and  care 
was  taken  to  prevent  damage  of  the  soil  sample  inside.   Laboratory 
samples  were  simply  removed  from  their  polyethylene  bags. 

The  next  step  was  to  prepare  tube  samples  for  use  in  the  Geonor 
sampling  apparatus.   The  ends  of  the  sample  were  trimmed  perpendicular 
to  the  axis  of  the  sample  cylinder,  leaving  a  maximum  height  of  about 
2.5  in  (6.4  cm).   The  ends  were  trimmed  with  a  wire  saw  for  the  wetter 
samples  and  with  a  sharp  pen  knife  for  the  drier  ones.   These  end 
trimmings  were  saved  for  the  companion  study  of  pore  size  distribution. 
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Care  was  taken  to  trim  all  test  specimens  from  as  close  to  the  mid- 
height  of  tube  samples  as  possible.   This  minimized  the  effects  of 
systematic  density  variation  in  a  compacted  soil,  as  described  by 
Gau  and  Olson  (1971). 

After  this  initial  trimming,  the  soil  sample  was  centered 
on  the  platen  of  the  Geonor  sampling  apparatus,  and  the  top  piston 
was  lowered  into  position;  the  piston  prevented  movement  of  the  sample 
as  the  cutting  ring  was  being  advanced.   This  situation  is  illustrated 
in  Figure  2.14.   The  cutting  ring  was  advanced  by  rapidly  pushing 
the  top  handle  down,  and  Figure  2.14  also  shows  the  cutting  ring  fully 
advanced  into  the  soil  sample.   The  cutting  ring  was  held  in  the  appar- 
atus by  spring  clips  such  that  when  the  handle  was  subsequently  raised, 
the  cutting  ring  remained  in  place  on  the  newly  trimmed  soil  cylinder. 
Trimmings  produced  by  this  cutting  action  were  then  used  for  water 
content  determination.   Figure  2.15  shows  the  newly  formed  soil  cylin- 
der and  cutting  ring,  with  trimmings  removed  from  the  perimeter.   The 
cutting  ring  was  then  advanced  completely  through  the  sample  by  hand 
with  the  aid  of  a  stainless  steel  tube. 

This  trimming  procedure  was  slightly  modified  for  laboratory 
samples  compacted  relatively  dry  at  the  highest  (C)  energy  level. 
The  cutting  ring  was  started  into  the  sample  with  the  Geonor  device. 
When  further  advance  could  not  be  achieved  by  hand,  the  sample  was 
placed  in  an  hydraulic  jack.   The  cutting  ring  was  then  advanced  through 
the  sample  with  the  jack,  aided  by  the  previously  mentioned  stainless 
steel  tube. 
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Figure  2.15   Cutting  Ring  on  Sample  with  Trimmings  Removed 
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After  removing  the  cutting  ring,  any  voids  around  the  perimeter 
of  the  cylinder  were  filled  with  soil  shavings.   These  voids  were 
usually  caused  by  rock  fragments  being  dislodged  from  the  soil  by 
the  cutting  ring;  the  size  and  number  of  these  voids  were  recorded 
in  the  notes  for  each  test. 

Next,  the  confining  ring  was  installed  by  hand  near  the  sample 
mid-height.   This  ring  was  advanced  in  the  same  direction  as  the  cut- 
ting ring,  and  much  care  was  taken  to  advance  it  parallel  to  the  axis 
of  the  soil  cylinder.   As  the  confining  ring  was  advanced,  it  shaved 
off  a  thin  layer  of  soil  from  the  circumference  and  thereby  produced 
a  tight  (but  not  excessively  tight)  fit  in  the  ring. 

Excess  soil  was  then  removed  from  the  ends  of  the  test  speci- 
mens.  A  wire  saw  was  used  for  wetter  specimens  while  drier  ones  re- 
quired cutting  and/or  scarifying  and  scraping  with  a  sharp  pen  knife. 
These  operations  were  performed  rather  gently  to  minimize  sample  dis- 
turbance.  Rock  fragments  which  protruded  beyond  the  ends  of  the  speci- 
men were  removed;  all  surface  voids  were  filled,  and  their  sizes  and 
locations  were  recorded.   Faces  of  the  test  specimen  were  made  flush 
using  a  straightedge,  and  then  they  were  covered  with  glass  plates 
to  prevent  drying  and  damage.   This  assembly  (confining  ring,  soil 
specimen,  and  glass  plates)  was  weighed  for  subsequent  water  content 
and  density  determinations,  and  then  transported  to  the  testing  area. 

It  should  be  noted  that  several  field  samples  were  discarded 
during  the  trimming  phase.   Samples  which  were  drier  than  about  15% 
moisture  content  were  usually  quite  brittle.   Some  of  these  were  of 
insufficient  height  due  to  defects  or  poor  recovery  from  the  field, 
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while  others  fractured  during  the  trimming  process.   Some  were  discard- 
ed because  a  large  volume  of  rock  fragments  had  to  be  removed  from 
the  test  specimen,  which  probably  made  the  sample  unrepresentative. 
This  latter  occurrence  was  independent  of  moisture  content. 

Field  samples  were  also  observed  to  vary  in  plasticity  from 
sample  to  sample.   For  this  reason,  liquid  limit  tests  were  performed 
for  each  field  sample  tested,  and  results  are  presented  in  Appendix 
A.   Plastic  limit  tests  were  performed  for  one  sample  from  each  series 
(i.e.,  each  combination  of  compactor  type,  water  content  level,  and 
energy  level);  the  resulting  average  and  range  of  Atterberg  limits 
are  presented  in  Table  2.1.   Results  of  these  tests  showed  that  the 
soil  was  somewhat  variable  but  the  average  values  remained  within 
the  bounds  of  the  AASHTO  A-6  classification. 

Constant  Volume  Swell  Test 

Constant  volume  swell  tests  (CVS)  were  conducted  on  both  field 
and  laboratory  compacted  samples.   Field  samples  were  randomly  selected 
for  testing  with  the  aid  of  a  random  number  table;  this  helped  reduce 
the  amount  of  bias  in  the  testing  procedure.   Laboratory  compacted 
samples  were  similarly  chosen  in  a  random  fashion.   A  total  of  149 
samples  were  tested,  with  93  of  these  compacted  in  the  field  and  56 
compacted  by  laboratory  methods. 

Two  Geonor  bench  top  consolidometers  were  located  in  a  con- 
trolled temperature  room  and  used  for  the  CVS  tests.   These  devices 
applied  load  to  the  soil  specimen  through  a  system  of  levers  and 
weights,  which  was  slightly  modified  for  the  CVS  tests.   The  original 
weight  hanger  was  replaced  with  a  plastic  bucket,  and  small  glass 
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beads  were  used  to  apply  load  instead  of  the  conventional  metal  weights. 
The  glass  beads  were  chosen  for  the  following  reasons:   (1)  they  easily 
provided  the  range  of  weights  (very  small  to  very  large)  which  were 
required  in  this  application;  (2)  they  applied  load  in  a  more  continu- 
ous manner  than  the  metal  weights  as  swell  pressure  developed;  (3) 
they  permitted  more  precise  control  of  the  amount  of  weight  applied, 
and  thus  the  measured  swell  pressure. 

The  consolidometers  were  calibrated  before  any  CVS  tests  were 
conducted,  as  Fredlund  (1969)  concluded  that  the  compressibility  of 
the  consolidometer  can  have  a  significant  effect  on  the  interpretation 
of  swell  test  data.   A  stainless  steel  cylinder  was  used  to  replace 
the  soil  specimen  for  these  tests.   Loads  were  applied,  deflections 
were  measured,  and  calibration  curves  were  drawn.   Results  of  these 
tests  were  used  during  CVS  tests  to  ensure  that  the  soil  specimen 
remained  nearly  constant  in  volume.   As  a  check,  calibrations  were 
performed  periodically  during  the  testing  phase. 

Preliminary  calibrations  were  conducted  with  cellulose  fiber 
filter  paper  placed  between  the  porous  stones  and  steel  cylinder. 
This  paper  produced  large  def elections,  portions  of  which  were  time- 
dependent.   The  filter  paper  was  replaced  with  thin  polycarbonate 
filter  membranes  (manufactured  by  the  Nuclepore  Corp.,  California) 
and  the  undersirable  effects  were  eliminated. 

After  the  trimming  process  was  completed,  the  soil  specimen 
and  confining  ring  were  placed  in  the  consolidation  chamber.   Porous 
stones  in  a  saturated  surface  dry  condition  and  dry  filter  membranes 
were  placed  over  the  ends  of  the  specimen.   A  collar  with  bolts  was  in- 
stalled to  hold  the  confining  ring  in  position,  and  a  loading  platen 


60 


was  placed  over  the  top  stone.   Figure  2.16  shows  the  consolidation 

chamber  both  before  and  after  assembly. 

The  chamber  was  then  placed  on  the  loading  frame  of  the  consol- 

idometer.   In  this  condition,  the  apparatus  was  balanced  to  apply 
no  load  to  the  specimen.   A  seating  load  of  1.3  psi  (9.0  kPa)  was 
applied  in  order  to  get  a  reliable  initial  reading  of  vertical  deflec- 
tion.  While  this  load  was  sufficient  to  seat  the  soil,  porous  stones, 
and  other  parts  of  the  consolidometer,  its  magnitude  and  duration 
were  not  large  enough  to  cause  any  appreciable  sample  disturbance. 
It  was  essential  to  get  a  reliable  initial  reading  because  all  subse- 
quent measurements  of  volume  change  were  based  upon  this  value. 

Next,  380  ml  of  demineralized  water  was  added  to  the  chamber, 
which  reduced  the  seating  load  to  1.0  psi  (6.9  kPa) .   The  water  was 
permitted  to  reach  the  soil  through  both  porous  stones,  and  swelling 
commenced  immediately  upon  inundation.   Loads  were  applied  as  swell 
pressures  developed,  and  values  were  tabulated.   Deflections  were 
measured  after  each  load  application  with  a  dial  gage  that  read  to 
0.0001  in  (0.0025  mm).   These  measured  values  were  then  compared  with 
the  deflections  required  to  maintain  a  constant  volume.   Time  and 
temperature  were  also  recorded.   The  consolidometer  with  chamber 
in  place  is  shown  in  Figure  2.17. 

During  the  test,  loads  were  added  only  if  the  deflection  mea- 
surements indicated  that  the  specimen  was  still  expanding;  maximum 
swell  pressure  was  achieved  when  the  soil  specimen  produced  no  further 
expansion  under  the  sustained  load.   Laboratory  samples  compacted 
very  dry  of  OMC  exhibited  a  sudden  reduction  in  volume  after  reaching 
maximum  swell  pressure.   When  this  occurred,  weights  were  removed 


61 


I 

en 
< 

u 

01 
4J 

< 

c 
<ti 

0) 

M 

o 

aa 


t: 

u 

c 

o 

•H 

4J 

n 


o 
to 

c 
o 
u 


kO 


a 

•H 


62 


u 
u 

V 

B 

c 

T3 


0 
ui 
c 

0 

u 

ft 

o 

u 

c 
<u 

ca 

M 

O 

c 

0 
(D 

c 


CN 


■H 


63 


until  no  further  volume  change  took  plase ,  and  thus  an  equilibrium 
swell  pressure  was  achieved.   The  maximum  or  equilibrium  load  was 
then  maintained  for  a  minimum  of  18  hours. 

After  completion  of  the  CVS  test,  the  chamber  was  removed  from 
the  consolidometer  frame  and  the  water  was  drained.   The  soil  sample 
was  then  extruded  from  the  confining  ring,  weighed  for  water  content 
and  density  determinations,  oven-dried  for  a  minimum  of  24  hours,  and 
weighed  again.   Compaction  curves  resulting  from  these  oven-dry  test 
specimen  values  are  presented  in  Figures  2.18  through  2.21  inclusively. 

For  laboratory  compacted  samples,  these  curves  are  very  similar 
to  the  ones  produced  from  compaction  mold  values.   This  indicates  that 
the  trimming  and  testing  procedures  cause  little  change  in  these 
measured  quantities.   For  the  field  compacted  samples,  however,  large 
differences  exist  between  field  nuclear  gage  and  test  specimen  values. 
This  is  evident  by  comparing  Figures  2.3  and  2.4  with  Figures  2.18 
and  2.19,  respectively.   The  figures  show  that  higher  densities,  less 
data  scatter,  and  more  well-defined  curves  appear  for  the  specimen 
values.   These  discrepancies  can  be  attributed  to  several  factors. 
Some  drying  may  have  occurred  during  storage,  thus  causing  increases 
in  dry  densities.   As  with  the  laboratory  compacted  samples,  the  trim- 
ming and  testing  procedures  probably  had  some  effect  on  w  and  y , .   All 
field  samples  were  taken  from  the  bottoms  of  impressions  left  by  the 
compactor.   Nuclear  gage  measurements  were  made  at  random  points  in 
the  fill  and  represented  a  volume  somewhat  different  from  that  of  the 
field  test  specimens.   As  a  result  there  is  a  large  probability  that 
test  specimens  soil  received  a  larger  applied  energy  than  did  the 
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Figure  2.18  Caterpillar  Compaction  Curves  for  Test  Specimen  Values 
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volume  surveyed  by  the  nuclear  gage.   In  addition,  because  the  two  sets 
of  data  were  taken  by  different  groups  of  people,  there  is  a  strong 
possibility  that  the  nuclear  gage  data  reported  herein  have  not  had 
applied  to  them  the  usually  routine  corrections  based  upon  sand-cone 
density  data. 

Figure  2.22  shows  how  swell  pressure  increases  with  time  for 
typical  samples.   The  pressure  builds  rapidly  upon  inundation,  and 
most  samples  achieve  maximum  values  in  from  two  to  six  hours.   Other 
investigators  (Seed  et  al.,  1961;  Parcher  and  Liu,  1965)  have  found 
swell  pressure  as  well  as  expansion  to  continue  increasing  for  several 
days,  which  was  not  observed  in  this  study.   This  contradiction 
probably  occurs  for  the  following  reasons:   (1)  the  samples  in  this 
study  are  smaller  in  size  than  those  used  by  others;  (2)  this  study 
permits  water  to  enter  the  test  specimen  through  both  ends,  while 
others  may  not;  (3)  the  soils  studied  may  have  idfferent  tendencies 
to  swell. 
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3  -  ANALYSIS  OF  DATA 


Introduction 


The  purpose  of  this  study  was  to  find  the  best  prediction  model 
for  the  tendency  to  swell  induced  by  compaction,  i.e.,  the  swell  pres- 
sure.  As  the  data  were  gathered,  it  became  apparent  that  variation 
existed  in  all  the  measured  parameters.   This  precluded  the  use  of 
a  functional  model,  i.e.,  one  in  which  all  data  points  lie  on  the  locus 
of  the  function.   Rather,  a  statistical  relation  was  derived  by  linear 
regression  analysis.   As  Neter  and  Wasserman  (1974)  stated, 

A  regression  model  is  a  formal  means  of  expressing  the  two  essen- 
tial ingredients  of  a  statistical  relation: 

(1)  A  tendency  of  the  dependent  variable   Y   to  vary  with  the 
independent  variable  or  variables  in  a  systematic  fashion; 

(2)  A  scattering  of  observations  around  the  curve  of  statisti- 
cal relationship. 

The  regression  equation  produces  a  curve  or  surface  of  "best 
fit!"  through  the  data  points.   These  "fitted  values"  are  not  adequate 
for  predictive  purposes,  as  the  analysis  must  also  assess  the  vari- 
ability of  the  data  set.   In  some  cases  the  variability  can  be  quite 
high  with  respect  to  the  expected  value.   As  the  amount  of  data  scatter 
increases,  the  usefulness  of  a  fitted  value  model  decreases.   The  pre- 
dictive models  for  swell  pressure  are  based  upon  a  regression  equation 
to  which  is  added  the  appropriate  amount  of  variability  (the  way  in 
which  this  variability  was  assessed  will  be  explained  in  a  later 
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section) .   The  resulting  model  is  then  expressed  as  a  fitted  or  ex- 
pected value  plus  a  variability  term.   This  then  creates  a  "ceiling 
value"  for  prediction  purposes.   This  ceiling  value  represents  the 
predicted  maximum  swell  pressure  that  should  be  observed  in  the  embank- 
ment.  The  engineer  can  use  this  to  guide  his  design  and/or  compaction 
specification  with  assurance  that  he  has  control  over  the  tendency 
to  swell  which  is  induced  by  compaction. 

This  analysis  paralleled  that  done  by  Price  (1978)  but  con- 
tained the  following  differences. 

(1)  This  analysis  was  extended  to  include  the  prediction  of  field 
swell  behavior  from  laboratory  compaction  variables. 

(2)  Some  of  the  independent  variables  were  represented  in  a  dif- 
ferent manner  (see  "Selection  of  Variables"  section). 

(3)  Both  sides  of  optimum  moisture  content  were  achieved  in  the 
field  and  were  accounted  for  in  this  analysis. 

(4)  The  Ridge  Regression  computer  program  was  not  used;  instead, 
variables  were  selected  using  hypothesis  testing. 

Selection  of  Variables 
The  basic  independent  variables  chosen  for  the  analysis  were 
water  content  as  a  deviation  from  the  optimum  moisture  content  (Aw) 
and  energy  ratio  (E  )  based  on  the  lowest  gross  input  energy  for  each 
compactor  type;  dry  density  (y    )  was  also  used  for  swell  pressure  pre- 
diction.  Aw  was  chosen  because  some  investigators  say  it  is  the  con- 
trolling variable  for  behavior  (Price,  1978).   Weitzel  (1979)  showed 
that  it  was  difficult  at  best  to  determine  the  amount  of  input  energy 
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which  was  actually  used  to  compact  the  soil  in  a  laboratory  kneading 
compactor.   The  largest  problem  came  from  the  inability  to  measure  the 
soil  rebound  as  the  tamping  foot  was  raised.   Similar  problems  occur 
with  laboratory  impact  compaction.   In  addition,  the  applied  field  en- 
ergy also  escapes  exact  measurement.   Therefore,  energy  ratio  was 
chosen  because  it  could  easily  be  quantified  for  each  type  of  compactor 
and  used  for  comparisons  (Table  3.1).   It  must  be  remembered,  however, 
that  energy  ratios  were  not  exactly  comparable  as  each  set  (lab  and 
field)  was  based  on  a  different  datum. 

Table  3.1  Energy  Ratios 


Compactor 
Type 


Quantity 
Measured 


Energy 
Level 


Caterpillar 

No.  of  Passes 

s 

4 

1.0 

Rascal 

8 

2.0 

16 

4.0 

Impact 

Applied  Energy 
(ft  •  lbs) 

247.5 

1.00 

412.5 

1.67 

1875.0 

7.58 

Kneading 

Gage  Pressure 
(psi) 

4 

1.00 

7 

1.75 

26 

6.50 

A 
E 
C 


A 
B 
C 


A 
B 
C 
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Analytic  Procedures 

The  first  step  in  the  analysis  was  to  plot  dry  density  against 
absolute  water  content  for  each  energy  level  and  type  of  compactor  to 
obtain  the  corresponding  compaction  curves.   Maximum  dry  density  and 
optimum  moisture  contents  (OMC)  were  selected  for  each  curve.   Dry 
density  was  then  plotted  against  water  content  from  optimum  (Aw)  and 
energy  ratio  (E  ) .   Swell  pressure  (SP)  was  plotted  against  water  con- 
tent from  optimum  (Aw),  dry  density  (y,),    and  energy  ratio  (E  ). 

The  curvilinear  nature  of  these  scatterplots  was  clearly  evi- 
dent.  Therefore,  second  order  terms  of  the  basic  independent  variables 
were  also  included  in  the  regression  analysis.   In  order  to  account 
for  possible  interactions  between  independent  variables,  they  were  com- 
bined two  at  a  time  to  a  maximum  power  of  two  for  each  independent  vari- 

2     2 
able  (e.g.,  Aw  x  y        represents  the  highest  order  used). 

After  selecting  the  appropriate  set  of  independent  variables, 

the  Purdue  computer  program  DRRSQU  was  employed.   This  gave  values  of 

2 

the  coefficient  of  multiple  determination  (R  )  in  ranking  order  for 

2 
all  possible  combinations  of  independent  variables.   R*~  measured  the 

proportion  of  total  variation  that  was  "explained"  or  "accounted  for" 

2 
by  the  set  of  variables  used  in  that  regression  model.   While  R  was 

the  only  statistical  descriptor  used  in  this  phase  of  the  analysis, 

it  was  not  the  only  criterion  used  for  selecting  the  final  prediction 

model.   Other  selection  criteria  which  were  used  will  be  explained  later 

in  this  section. 

2 

Equations  which  showed  a  reasonably  high  R  value  were  sub- 
jected to  additional  analysis  using  the  REGRESSION  computer  program 
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(an  SPSS  program  developed  by  Nie,  et  al. ,  1975).   Only  those  equations 
with  three  or  fewer  independent  variables  were  considered.   Increasing 

the  number  of  independent  variables  beyond  three  did  not  appreciably 

2 
affect  the  R  values,  but  it  did  serve  to  complicate  the  model.   Plots 

of  the  residuals  were  made  and  studied,  and  no  trends  were  observed. 
This  indicated  that  the  residuals  behaved  in  a  random,  independent  man- 
ner and  that  the  error  variance  was  constant;  meeting  these  conditions 
were  basic  requirements  of  least  squares  regression  analysis.   For  each 
of  the  models  selected,  the  analysis  then  determined  if  the  entire  equa- 
tion, as  well  as  each  term  contained  therein,  significantly  contributed 
toward  the  description  of  the  true  relationship.   All  such  hypothesis 
testing  used  a  95%  confidence  coefficient,  which  implies  that  the  risk 
of  making  a  Type  I  statistical  error  was  controlled  at  5%.   A  Type  I 
error  in  this  context  would  be  to  claim  that  a  regression  equation, 
or  variable (s)  in  an  equation,  did  not  contribute  toward  the  true 

description  when  in  fact  it  did  contribute. 

2 
It  was  mentioned  earlier  that  a  high  value  of  R   alone  did 

not  provide  an  adequate  criterion  for  selecting  the  final  prediction 

model.   The  analysis  described  in  the  preceding  paragraph,  along  with 

the  criteria  in  Table  3.2  were  used  to  select  the  final  prediction 

model.   Each  category  of  data  required  somewhat  different  treatment. 

Therefore,  more  detailed  descriptions  follow. 

Dry  Density  and  Swell  Pressure  Magnitude 
Swell  pressure  magnitude  was  analyzed  first.   Data  were  first 
separated  by  compactor  type  and  whether  compaction  was  performed  in 
the  lab  or  in  the  field.   Water  content  and  dry  density  had  been 
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Table  3.2  Criteria  Used  for  Selecting  Prediction  Models 

1.  The  Model  must  be  relatively  simple  in  order  to  encourage  use 

(this  includes  using  the  minimum  number  of  independent  variables 
and  equations) . 

2 

2.  The  model  must  provide  a  good  statistical  fit  (this  refers  to  R  , 

residual  analysis,  and  hypothesis  testing). 

3.  The  shape  of  the  regression  surface  should  reflect  the  trends 
shown  by  the  data  points. 

4.  The  model  must  describe  the  variability  associated  with  the  data. 

5.  The  model  must  be  able  to  predict  field  response  from  laboratory 
data. 
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measured  in  two  ways  for  each  compactor  type,  but  only  one  set  of  mea- 
surements could  be  used  for  the  analysis.   The  set  of  measurements 
chosen  for  each  case  was  based  on  the  ability  to  describe  the  statis- 
tical relationship  as  well  as  the  ease  of  gathering  data  for  future 
application  of  this  research. 

For  laboratory  compaction,  water  content  and  dry  density  were 
measured  for  both  the  compaction  mold  and  the  test  specimen.   These 
compaction  curves  are  shown  in  Figures  2.11  and  2.12,  and  Figures  2.20 
and  2.21,  respectively.   From  these  curves  it  can  be  seen  that  the  place 
of  measurement  caused  little  difference  in  the  measured  values.   Since 
the  measurements  were  nearly  equal  compaction  mold  values  were  used 
for  the  analysis.   It  was  felt  that  little,  if  any,  accuracy  was  lost 
by  choosing  these  instead  of  test  specimen  values. 

Plots  of  swell  pressure  versus  water  content  deviation  from 
optimum  (Aw)  for  impact  and  kneading  compaction  are  given  in  Figures 
3.1  and  3.2  respectively.   The  data  suggest  that  swell  pressure  is  a 
continuous  function  of  Aw   for  each  energy  level.   As  a  result,  the 
data  sets  for  laboratory  compaction  did  not  require  further  breakdown 
into  subsets.   Resulting  regression  equations  are  presented  in  Table 
3.3.   It  is  interesting  to  note  that  the  same  variables  appear  in  the 
equations  for  both  types  of  laboratory  compaction,  but  with  different 
regression  coefficients. 

For  field  compaction,  water  content  and  dry  density  measure- 
ments were  taken  in  the  field  by  nuclear  devices  and  also  from  test 
specimens  in  the  laboratory.   Initial  analyses  were  performed  using 
the  field  measurements.   As  Figures  2.3,  2.4,  2.18  and  2.19  showed, 
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large  differences  existed  between  laboratory  and  field  measurements 
for  each  field  compactor  type.   Probable  causes  for  these  discrepancies 
were  previously  discussed.   Variations  in  water  content  and  dry  density 

existed  within  each  combination  of  water  content  level,  energy  level, 

* 
and  compactor  type.    This  was  evidenced  for  both  field  and  test  speci- 
men values  in  Table  2.3  and  Appendix  A,  respectively.   These  variations, 
in  turn,  produced  variability  in  swell  pressure  measurements  within 
a  series.   It  was   impossible  to  directly  relate  a  specific  field  mea- 
surement with  any  given  test  specimen,  and  therefore  average  values 
of  field  water  content  and  dry  density  had  to  be  used  for  each  sample 
series.   These  averages  limited  the  analysis  to  only  five  discreet  com- 
binations of  water  content  and  dry  density  per  energy  level  when,  in 
fact,  the  test  specimens  produced  a  minimum  of  14  such  combinations. 
At  each  average  field  water  content  the  swell  pressure  varied  greatly 

(see  Figures  3.3  and  3.4),  the  statistical  relationship  was  obscured 

2 
and  poor  R  values  resulted.   Because  of  this,  it  was  decided  to  investi- 
gate swell  pressure  as  a  function  of  test  specimen  water  content  and 
dry  density. 

Figure  3.5  shows  swell  pressure  plotted  against  test  specimen 
water  content  deviation  from  OMC  for  the  Caterpillar  compactor.   A  dis- 
tinct cusp  is  evident  near  the  OMC  with  different  trends  in  the  data 
points  on  each  side  of  it.   Polynomial  regression  equations  do  not  ade- 
quately describe  these  features  over  all  water  contents  ;   therefore  the 
data  are   broken  into  two  subsets.   All  data  points  dry  of  OMC  are 
grouped  together  for  analysis.   The  resulting  equation,  shown  in 


A  combination  of  water  content  level,  energy  level,  and  com- 
pactor type  will  also  be  referred  to  as  a  series  or  sample  series. 
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Figure  3.3   Caterpillar  Swell  Pressure  Data  for  Field  Nuclear 
Gage  Values 
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Table  3.3,  reflects  the  linear  trend  in  these  data.  For  points  wet  of 
OMC,  similar  analysis  results  in  a  quadratic  regression  equation.   Peak 
values  of  swell  pressure  occur  slightly  dry  of  OMC,  and  are  also  in- 
cluded in  the  analysis  of  wet  side  data.   They  serve  to  raise  this  re- 
gression line  such  that  it  intersects  the  dry  side  regression  line 
closer  to  OMC  and  at  a  higher  swell  pressure.   This  is  believed  to  be 
appropriate  and  desirable  for  the  given  data  set.   It  does  not  appreci- 
ably change  any  statistical  parameters.   It  should  be  noted  that  water 
content  from  optimum  (Aw)  is  the  only  significant  descriptor  of  swell 
pressure  induced  by  the  Caterpillar  compactor. 

Swell  pressure  data  for  the  Rascal  compactor  are  shown  in  Fig- 
ures 3.6  and  3.7.   More  data  scatter  is  apparent  here  than  for  the 

2 
Caterpillar  data  which  resulted  in  consistently  lower  R   values.   Sim- 
ilar trends  in  the  data  occur  for  the  A  and  B  energy  levels.   To 
adequately  describe  them,  they  had  to  be  separated  from  the  C  energy 
level.   The  data  from  the  C  energy  level  show  a  flat  response  with 
much  scatter  and,  as  a  result,  a  significant  regression  equation  could 
not  be  found  for  this  energy  level.   The  regression  equations  are  again 
given  in  Table  3.3. 

Dry  density  was  a  significant  variable  in  many  of  the  regression 
relations  for  swell  pressure.   Because  of  this  it  was  also  subjected 
to  regression  analysis.   Water  content  and  dry  density  measurements 
were  taken  from  the  compaction  mold  for  laboratory  compacted  samples 
and  from  test  specimens  for  field  compaction.   Data  from  each  compactor 
type  were  again  analyzed  separately. 
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After  compaction  curves  were  drawn  and  optimum  moisture  contents 
selected,  each  data  set  was  analyzed.   Initial  analyses  showed  that 

quadratic  models  poorly  described  the  data  over  all  water  contents, 

2 
but  still  produced  a  high  value  of  R  .   For  this  reason  the  data  sets 

were  subdivided  into  dry  of  OMC  and  wet  of  OMC  subsets.   The  results 

of  the  regression  analysis  are  shown  in  Table  3.3.   It  should  be  noted 

that  both  the  impact  and  kneading  compaction  data  sets  contained  one 

point  wet  of  OMC  which  was  not  tested  for  swell  pressure.   These  data 

points  were  needed  to  help  define  the  compaction  curves  but  not  the 

swell  pressure  curves. 

Variability  of  Dry  Density  and  Swell  Pressure 
The  statistical  analysis  for  variabilities  associated  with  dry 
density  and  swell  pressure  also  paralleled  that  of  Price  (1978) .   The 
variability  in  water  content  and  dry  density  for  all  of  the  compaction 
conditions  studied  influenced  other  soil  properties  such  as  strength 
and,  in  this  case,  swell  pressure.   Therefore,  swell  pressure  variabil- 
ity was  not  determined  for  specific  values  of  water  content  and  dry 
density,  but  rather  for  a  range  of  these  values  at  a  particular  com- 
pactive  effort. 

Two  ways  of  measuring  the  variability  of  a  data  set  are  of  in- 
terest in  this  report.   The  variance  of  n  observations  (a  data  set) 
essentially  measures  the  average  deviation  from  their  mean,  Y.   The 
standard  deviation  is  the  square  root  of  this  variance.   Another  mea- 
sure, the  error  mean  square  (MSE) ,  assesses  the  deviation  of  observed 
values  (Y.)  around  the  fitted  regression  line.   The  difference  between 
these  two  measures  is  illustrated  in  Figure  3.8.   The  first  approach 
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measures  the  scatter  around  a  single  value.   This  study,  however,  indi- 
cates that  dry  density  and  swell  pressure  vary  in  a  systematic  fashion 
with  other  variables.   This  makes  it  necessary  to  determine  the  variabil- 
ity at  several  response  levels.   Therefore,  the  error  mean  square  method 
is  appropriate  for  application  in  this  analysis. 

All  regression  models  estimate  the  true  relationship  among  the 
variables  in  a  data  population;  this  produces  an  inherent  uncertainty 
in  the  accuracy  of  the  model.   As  a  result,  these  estimates  are  most 
accurate  near  the  average  values  of  the  dependent  variable  and  all  inde- 
pendent variables.   Accuracy  of  estimation  decreases  with  distance  from 
these  mean  values.   Therefore,  for  a  given  degree  of  confidence,  the 
confidence  interval  must  increase  in  size  with  distance  from  these  means. 

For  a  general  linear  regression  model,  the  following  formula 
is  used  to  give  an  interval  estimation  (confidence  interval)  of  the 
mean  response  E(Y  )  at  one  point: 

Yu  -  t(l  -  a/2;  n  -  p)  •  S(Y,  )  <  E(Y,  )  <  Y,  +  t(l  -  a/2;  n  -  p)  •  S(Y,  ) 
h  h  —    n  —  n  n 


Where i 


E(Y.  )  =  the  expected  mean  response  for  the  given 
h 

levels  of  the  independent  variables  [XI 
Y  =  the  fitted  regression  value  for  the  given 
values  of  [)t  ] 
n  =  the  number  of  observations  in  the  data  set 
p  -  1  =  the  number  of  independent  variables  in  the 
regression  equation 
a  =  probability  of  committing  a  Type  I  error 
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t(l  -  a/2;  n  -  p)  =  appropriate  t-statistic 

S(Y  )  =  square  root  of  S  (Y  ) 

S2(Y  )  =  MSE[X"  (X'X)"^,  ]   =  estimated  variance  of  1 
h        h       h  h 

MSE  =  error  mean  square.   This  value  was  obtained 
from  the  REGRESSION  computer  program  output 
[X]  =  matrix  of  observed  values  of  the  indepen- 
dent variables  in  the  regression  model 
[X] '  =  [X]  transposed 

[X^]  =  column  vector  containing  values  of  the  in- 
dependent variables  for  the  response  level 
under  investigation 
[X^] '  =  [X  ]  transposed. 

This  formula  for  variability  is  usually  evaluated  for  E (Y  )  by  using  a 
fixed  set  of  independent  variables,  [X,] .   This  produces  a  point  on  the 
response  surface,  Y  ,  and  also  the  confidence  interval  at  that  re- 
sponse level.   Although  this  is  sufficient  for  most  applications  of 
this  formula,  it  does  not  suffice  for  studies  of  soil  compaction. 

The  state-of-the-art  of  compaction  cannot  guarantee  precise 
values  of  [)L]     (i.e.,  w,  y    ,   E  ).   Rather,  these  quantities  must  be 
viewed  as  expected  values,  at  each  response  level,  with  some  inherent 
variability.   These  additional  deviations  must  also  be  included  in  the 
total  variability  of  the  dependent  variable. 

In  this  study,  although  each  test  pad  was  constructed  to  have 
a  nearly  constant  water  content,  a  relatively  large  variation  did  occur. 
The  range  in  water  content  was  calculated  for  each  level  of  compactive 
effort  and  compactor  type.   The  average  range  for  both  compactors  was 
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approximately  2.8%  for  the  test  specimens  and  3.5%  for  field  nuclear 
gauge  measurements.   Since  these  values  will  not  apply  to  all  compaction 
projects,  dry  density  and  swell  pressure  variabilities  were  calculated 
for  several  different  water  content  ranges.   For  illustrative  purposes 
herein  an  expected  range  of  2.8%  was  used. 

The  amount  of  variability  associated  with  compactive  effort 
could  not  be  quantified  and,  therefore,  was  assumed  to  be  zero  for  this 
analysis.   The  regression  equations  presented  in  the  previous  section 
were  used  to  determine  the  expected  dry  density  for  each  level  of   Aw 
and   E  .   A  95%  confidence  criterion  was  arbitrarily  chosen  for  deter- 
mining  dry  density  variability.   For  each  value  of   Aw  chosen,  the  95% 
confidence  interval  of  dry  density  was  calculated  for  Aw  plus  and  Aw 
minus  the  half-range  variation  in  water  content  (equal  to  1.4%  here). 
The  larger  of  these  two  values  was  then  chosen  as  the  expected  dry 
density  variability,  V(v •  ) .   This  method  was  also  used  to  determine 
the  swell  pressure  variability  for  the  Caterpillar  compactor,  since   Aw 
was  the  only  significant  variable  in  the  regression  model. 

The  following  method  was  used  to  estimate  swell  pressure  vari- 
ability for  the  other  compactor  types.   Since   Aw   and  dry  density  ap- 
pear in  these  models,  the  variability  of  both  was  included  in  the 
analysis.   It  was  assumed  that   Aw   and  dry  density  varied  independently 
at  each  response  level.   Therefore,  swell  pressure  variability,  V(SP) , 
was  calculated  for  the  four  conditions  shown  below: 
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Aw                   y,               V(SP) 
—  a  


Aw  +  1.4           y,  +  V(Y,)  V(SP)n 

a      a  1 

Aw  +  1.4          y      -  V(y  J  V(SP)n 

a      a  2 

Aw  -  1.4          Y,  +  V(Y,)  V(SP), 

a      a  3 


Aw  -  1.4  y,  -  V(Y,)  V(SP)„ 

a      a  4 


The  largest  value  of  V(SP)  obtained  from  the  four  combinations  was 
chosen  as  the  expected  swell  pressure  variability  for  that  response 
level.   This  analysis  was  also  done  for  a  95%  confidence  criterion. 

A  two-sided  confidence  interval  was  used  for  calculation  of  dry 
density  variability.   The  appropriate  t-statistic  was  t(l  -  a/2;  n  -  p) 
or  t (0.975;  n  -  p) .   The  purpose  of  the  swell  pressure  variability  was 
to  create  a  "ceiling  value"  for  swell  pressure,  which  required  the  use 
of  a  one-sided  confidence  interval.   In  this  case  however,  four  values 
of  variability  were  calculated  at  each  response  level.   In  order  to 
maintain  an  overall  a-level  of  0.05,  each  variability  was  evaluated 
at  a  level  of  a/4.   The  appropriate  t-statistic  was  t(l  -  a/4;  n  -  p) 
or  t(0.9875;  n  -  p) .   This  was  analogous  to  using  a/2  in  a  two-sided 
test  at  a  single  point.   As  previously  mentioned,  swell  pressure 
variability  for  the  Caterpillar  compactor  needed  to  be  evaluated  at 
only  two  points.   The  appropriate  t-statistic  therefore  was  t (0.975; 
n  -  p)  . 

The  analysis  presented  here  was  evaluated  for  one  water  content 
range  (2.8%)  and  confidence  coefficient  (95%).   The  computer  programs 
presented  in  Appendix  B  can  be  adjusted  to  accomodate  the  expected 
conditions  of  future  projects. 
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4-DISCUSSION  OF  RESULTS 

Dry  Density  and  Swell  Pressure  Magnitude 
Dry  density  is  a  significant  prediction  variable  in  swell  pres- 
sure equations  for  three  of  the  four  compactor  types  used  in  this 
study;  thus,  prediction  models  for  dry  density  are  also  developed 
and  presented.   For  these  three  compactors,  models  for  dry  density 
are  needed  to  predict  both  the  magnitude  and  the  variability  of  swell 
pressure.   These  models  also  enable  the  prediction  of  field  dry 
density  for  inspection  testing,  as  outlined  in  the  "Application  of 
Results".   The  dry  density  models  can  be  directly  substituted  into 
those  for  swell  pressure,  but  this  does  not  simplify  the  calculation 
of  SP.   Such  a  substitution  requires  the  use  of  two  swell  pressure 
models  (one  for  dry  of  OMC  and  one  for  wet  of  OMC)  where  one  had  pre- 
viously sufficed. 

Caterpillar  Compactor 
Regression  results  show  that  dry  density  is  not  a  function 
of  energy  ratio  for  points  dry  of  OMC.   The  slope  of  this  regression 
line  (the  regression  coefficient  of   Aw)  is  shown  in  Table  3.3  and 
Figure  4.1,  and  is  approximately  three  times  steeper  than  those  for 
the  other  compaction  types.   These  features  could  be  a  result  of: 
(1)  a  relatively  small  number  of  data  points  dry  of  OMC  for  this 
compactor  (only  12  points);  (2)  the  narrow  range  of  water  contents 


94 


i"_ 


in 
o 

ii 

O 

a 


125 


120 


115 


no  - 


>" 

CO 

z 

UJ 
Q 

> 
cr 
Q 

Q 
UJ 

1- 

o 

UJ 

a. 
x 

UJ 


105- 


100 


95 


ENERGY 
LEVEL 

O  A 
A  B 
□    C 


t 


r  =  5.863Aw  +  118.3 
d 


— i— 
-6 


- 1.609  Aw  +  1.006  ER+  114.3 


T" 

-4 


-2 


T" 

4 


-T- 

6 


"T" 

8 


— i 
10 


WATER   CONTENT   DEVIATION     FROM    OMC,      Aw(%) 


Figure  4.1   Caterpillar  Dry  Density  Relationships  for  Test  Specimen 
Values 


95 

achieved  dry  of  OMC  in  the  field  (all  Aw  wetter  than  -2%);  (3)  an 
exceptionally  large  amount  of  variability  in  the  data;  (4)  the  inabil- 
ity of  this  compactor  to  achieve  consistent  densities  in  this  soil 
when  compacted  dry.   Wet  of  OMC,  dry  density  is  a  function  of  both 
Aw  and  ER,  and  regression  results  are  comparable  to  those  for  the 
other  compactor  types.   Regression  relations  for  the  Caterpillar  com- 
pactor dry  density  are  shown  in  Figure  4.1. 

Dry  of  OMC,  swell  pressure  is  a  function  of  Aw  only  and 
reaches  a  maximum  near  OMC.   This  trend  is  contrary  to  results  obtained 
for  the  other  compactor  types.   In  general,  swell  pressure  has  been 
observed  to  increase  to  a  plateau  as  water  content  decreases.   This 
discrepancy  could  be  a  result  of  the  same  four  factors  mentioned  in 
the  preceding  paragraph,  especially  item  (4) .    Nalezny  and  Li  (1967) 
postulated  that  if  dry  density  becomes  quite  low  at  low  molding  water 
contents,  many  large  pores  exist.   As  water  is  imbided,  the  aggregates 
expand  into  the  pore  space,  which  reduces  the  measured  swell  pressure. 

Wet  of  OMC,  swell  pressure  is  a  curvilinear  function  of  Aw. 
Again,  dry  density  and  energy  ratio  are  not  significant  variables. 
The  swell  pressure  decreases  as  water  content  increases,  which  is 
consistent  with  other  compactor  types  used  in  this  study,  as  well 
as  with  results  presented  by  other  investigators.   However,  the  re- 
gression relation  increases  in  magnitude  beyond  approximately  Aw  = 
7%.   This  is  believed  to  be  a  consequence  of  the  parabolic  equation 
rather  than  of  the  actual  data  points,  and  the  equation  should  not  be 
used  beyond  this  value  of  Aw.   The  horizontal  asymptote  shown  in  Figure 
3.5  is  believed  to  be  more  appropriate  in  this  range. 
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Rascal  Compactor 

Regression  analyses  for  dry  density  reveal  that  Aw  and  E 

R 

are  significant  prediction  variables  both  wet  and  dry  of  OMC.   Table 
3.3  and  Figure  4.2  present  these  results.   The  figure  reveals  that 
more  data  scatter  exists  for  points  dry  of  OMC  than  for  those  wet 
of  OMC. 

Swell  pressure  is  broken  into  two  data  sets  for  analysis; 
the  A   and   B   energy  levels  form  one  set  while  the   C   level  com- 
prises the  other.   Figure  3.6  shows  that  the  A  and   B   energy  levels 
exhibit  similar  trends,  with  generally  higher  swell  pressures  being 
displayed  for  the  B   level.   This  swell  pressure  regression  rela- 
tion decreases  as   Aw   increases,  reaching  a  minimum  value  slightly 
wet  of  OMC;  at  larger  values  of   Aw,  swell  pressure  again  increases. 
The  increase  in  swell  pressure  wet  of  OMC  is  contrary  to  results  of 
studies  on  lab  compacted  soils  found  in  the  literature  (Seed  and  Chan, 
1961;  Nalezny  and  Li,  1967).   This  phenomenon  could  be  the  result 
of:   (1)  the  large  amount  of  data  scatter  produced  by  this  compactor, 

as  discussed  in  the  "Analysis  of  Data"  chapter  and  evidenced  by  the 

2 
low   R   values;  (2)  changes  in  soil  fabric  that  may  occur  at  certain 

water  contents  for  this  vibratory  compactor.   It  should  be  noted  that 

White  (1980)  found  no  significant  differences  in  pore  size  distribution 

for  samples  compacted  wet  of  OMC  by  this  roller.   Table  3.3  indicates 

2 
that   Aw   and  E    are  not  significant  variables  in  the  equation 

presented,  but  these  variables  are  retained  in  the  model  for  the  fol- 

2 

lowing  reasons:   (1)   Aw  provides  the  curvature  of  the  model  and 

this  is  believed  to  be  appropriate;  (2)   E  produces  two  distinct 

R 

curves,  and  the  data  trends  suggest  this. 
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For  the  C  energy  level,  a  significant  regression  equation 
for  swell  pressure  could  not  be  found.   This  is  due  to  the  large  amount 
of  data  scatter  and  the  relatively  flat  (horizontal)  trend  shown  in 
Figure  3.7.   Possible  causes  for  the  data  scatter  are  discussed  in 
preceding  sections.   The  flat  trend  in  swell  pressure  is  an  interesting 
occurrence,  but  attempts  at  explanation  only  produced  speculations. 

It  should  be  reiterated  that  a  significant  regression  equation 
for  swell  pressure  could  not  be  found  for  the  Rascal   C   energy  level. 
Thus,  the  equation  presented  herein  should  be  used  with  caution.   It 
is  believed,  however,  that  this  equation,  when  used  with  its  associated 
variability,  provides  a  realistic  model  for  the  swell  pressure. 

Laboratory  Compactors 

Laboratory  impact  and  kneading  compaction  methods  produced 
similar  results.   The  comapction  curves,  shown  in  Figures  2.11  and 
2.12,  are  nearly  identical,  as  are  their  maximum  dry  densities  and 
optimum  moisture  contents  for  corresponding  energy  levels  (see  Table 
4.1).   In  general,  the  slopes  of  the  kneading  compaction  curves  are 
somewhat  steeper  than  those  of  the  impact  compaction  curves.   All  of 
these   features  also  occur  for  the  regression  results  presented  in 
Figures  4.3  and  4.4,  and  Table  3.3. 

Swell  pressure  responses  are  also  very  similar  for  the  two 
compactor  types  and  are  shown  in  Figures  3.1  and  3.2.   The  regression 
equations  presented  in  Table  3.3  include  the  same  independent  vari- 
ables, but  with  different  coefficients.   These  similarities  are  not 
surprising,  as  White  (1980)  found  no  significant  differences  in  pore 
size  distributions  of  the  fabric  produced  by  the  two  methods  of 
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Table  4.1   Summary  of  Optimum  Moisture  Contents  and  Maximum  Dry 
Densities  Obtained  from  Compaction  Curves 


Location 

Compactor 

of 

Energy 

OMC 

Maximum  Dry 

Type 

Measurement 

Level 

(%) 

Density  (pcf) 

Caterpillar 

Test 

A 

15.0 

114.5 

Specimen 

B 

13.6 

118.0 

C 

13.4 

118.5 

Field 

A 

15.7 

105.5 

Nuclear  Gage 

B 

15.6 

110.9 

C 

15.5 

111.5 

Rascal 

Test 

A 

15.3 

114.0 

Specimen 

B 

14.8 

115.0 

C 

14.0 

117.0 

Field 

A 

16.2 

103.7 

Nuclear  Gage 

E 

15.4 

107.2 

C 

15.0 

110.6 

Impact 

Compaction 

A 

17.8 

109.5 

Mold 

B 

16.6 

113.8 

C 

12.3 

125.0 

Test 

A 

16.6 

111.1 

Specimen 

B 

15.6 

114.2 

C 

12.5 

122.7 

Kneading 

Compaction 

A 

18.6 

110.7 

Mold 

E 

16.5 

114.9 

C 

12.6 

124.1 

Test 

A 

17  7 

109.8 

Specimen 

E 

16.2 

113.9 

C 

12.8 

122.9 
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compaction.   Seed  and  Chan  (1961)  also  found  that  differences  in  struc- 
ture, strength,  shrinkage,  and  swell  resulting  from  impact  and  knead- 
ing methods  are  small. 

The  data  trends  and  regression  curves  are  nearly  identical 
for  the  A  and  B  energy  levels  for  the  two  compactors.   These  re- 
lationships are  of  the  same  general  shape  as  those  found  in  the  liter- 
ature (Nalezny  and  Li,  1967).   For  the  C  energy  level,  however, 
some  discrepancies  appear.   The  data  points  suggest  a  sharp  peak  near 
a  water  content  2%  dry  of  OMC;  swell  pressure  drops  rapidly  with  dis- 
tance from  this  peak,  especially  on  the  wet  side.   These  peaks  resemble 
that  found  near  OMC  for  the  Caterpillar  compactor  data,  but  the  poly- 
nomial regression  equations  found  for  laboratory  compacted  samples 
could  not  simulate  these  sharp  peaks.   Another  discrepancy  arises 
from  the  swell  pressure  values  near  this  peak.   Dry  of  this  peak, 
impact  compaction  produces  swell  pressures  which  are  10  psi  to  35 
psi  (69  kPa  to  240  kPa)  higher  than  those  for  kneading  compaction. 

Samples  compacted  very  dry  of  OMC  exhibit  a  tendency  to  col- 
lapse in  the  CVS  test.   For  these  samples  the  swell  pressure  increases 
with  time  in  a  normal  fashion,  reaches  some  maximum  value,  then  begins 
to  rapidly  reduce  in  volume  under  this  sustained  load.   At  this  time, 
the  load  on  the  sample  is  decreased  until  no  further  value  change 
occurs,  and  thus  an  equilibrium  swell  pressure  is  achieved.   In  Fig- 
ures 3.1  and  3.2,  collapsed  samples  are  indicated  by  two  symbols  con- 
nected by  a  line.   The  open  symbol  shows  the  maximum  swell  pressure 
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prior  to  collapse,  while  the  darkened  symbol  represents  the  equilibrium 
value  after  collapse.   These  figures  also  show  that  the  onset  of  col- 
lapse occurs  slightly  wetter  for  kneading  samples  than  for  impact 
samples.   This  difference  is  also  probably  due  to  small  differences 
in  dry  density  and  soil  structure  produced  by  the  two  methods  of  com- 
paction.  An  explanation  of  this  collapse  phenomenon  is  given  below. 

Barden  and  Sides  (1969)  suggested  that  the  collapse  phenomenon 
occurs  primarily  on  a  macroscopic  level,  and  is  not  due  to  micro- 
scopic changes  from  a  floculated  to  a  dispersed  soil  structure.   Com- 
paction on  the  dry  side  produces  an  open  structure  for  the  macroscopic 
particles  (macropeds) .   These  macropeds  resist  distortion  during  com- 
paction because  of  their  high  shear  strength,  thus  resulting  in  large 
air-filled  interpedal  channels.   On  the  wet  side,  however,  the  macro- 
peds have  a  lower  shear  strength  and  easily  deform  into  a  more  homogen- 
eous mass  with  few  large  voids.   Observations  made  during  this  study 
also  revealed  the  presence  of  individually  visible  macropeds  for  dry- 
side  laboratory  compaction,  while  wet-side  samples  were  of  a  much 
more  homogeneous  nature.   White  (1980)  also  confirmed  the  higher  fre- 
quency of  large-size  pores  for  dry-side  compaction  of  this  soil. 

As  water  enters  the  relatively  dry  macropeds,  it  breaks  down 
the  surface  tension  menisci  which  results  in  a  higher  (less  negative) 
pore  water  pressure  and  a  lower  effective  stress.   As  the  swell  pres- 
sure develops,  the  higher  applied  load  eventually  causes  the  now  softer 
macropeds  to  displace  into  the  available  pore  space,  thus  causing 
collapse.   Mishu  (1963)  concludes  that  collapse  occurs  when  the  combin- 
ation of  internal  and  applied  stresses  are  such  that  slippage  between 
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particles  occurs.   As  the  effective  stress  decreases  upon  wetting 
a  reduced  resistance  to  particle  slippage  results.   When  the  applied 
load  is  large  enough,  slippage  occurs;  at  lower  applied  loads  no 
particle  slippage  takes  place  and  swell  can  result.   For  samples  com- 
pacted very  dry  of  OMC,  the  combination  of  internal  and  applied 
stresses  in  the  CVS  test  eventually  reaches  the  onset  of  particle 
slippage,  and  collapse  occurs. 

Variability  of  Dry  Density  and  Swell  Pressure 
Variations  in  the  results  of  compaction  are  expected  due  to 
inherent  variations  in  soil  charcteristics,  compaction  processes, 
and  testing  procedures.   Possible  causes  of  these  variations  have 
been  discussed  in  preceding  sections  of  this  work.   These  variations 
in  turn,  produce  variability  in  the  mechanical  properties  of  the  soil, 
such  as  swell  pressure  and  strength. 

Before  the  amount  of  variability  in  swell  pressure  was  deter- 
mined for  the  field  compacted  soil,  the  variations  in  the  compaction 
variables  were  assessed.   The  compaction  variables  used  were  water 
content  deviation  from  OMC,  input  energy  ratio,  and  dry  density  where 
applicable.   The  amount  of  swell  pressure  variability  caused  by  input 
energy  variations  could  not  be  quantified.   The  variations  in  Aw 
were  expressed  as  half-ranges,  with  average  values  found  for  field 
nuclear  gage  and  test  specimen  measurements.   These  half-ranges  were 
then  used  to  determine  the  amount  of  variability  associated  with  dry 
density.   Finally,  swell  pressure  variability  was  assessed  using  the 
computer  programs  given  in  Appendix  B.   A  more  thorough  discussion 
of  these  procedures  is  given  in  the  "Analysis  of  Data"  chapter. 
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Other  Associated  Relationships 
Figures  4.5  through  4.14  show  some  interesting  relationships 
derived  from  the  data.   Such  relationships  may  help  to  predict  compac- 
tion characteristics  and  swell  pressure  behavior  with  a  minimal  amount 
of  physical  testing.   Figure  4.5  shows  swell  pressure  plotted  against 
dry  density  for  impact  compaction.   It  can  be  seen  that  a  "loop"  func- 
tion appears;  dry  of  OMC  swell  pressure  increases  with  increasing 

density  to  a  point  near   Yj    >    and  then  decreases  as   v -.  also  de- 

amax  d 

creases  wet  of  OMC.   This  relationship  shows  the  interdependence  be- 
tween water  content  and  dry  density,  and  its  effect  on  swell  pressure. 
This  interdependence  gives  insight  as  to  why   Aw   and  y        appear 
in  an  interaction  term  in  several  of  the  regression  equations  for 
swell  pressure.   Similar  trends  in  data  are  found  for  all  compactor 
types  used  in  this  study,  but  they  are  most  clearly  illustrated  by 
the  impact  compaction  results. 

Figure  4.6  shows  that  maximum  swell  pressure  decreases  as 
OMC  increases.   For  a  given  soil,  OMC  can  increase  only  if  input  energy 
decreases,  which  also  results  in  lower  dry  densities.   This  combination 
of  lower  densities  and  higher  water  contents  can  only  reduce  the  tend- 
ency to  swell.   This  is  consistent  with  evidence  presented  in  the 
literature  (Seed  et  al.,  1961). 

Figure  4.7  presents  the  relationship  between  the  maximum  swell 
pressure  and  that  which  develops  at  OMC.   The  Caterpillar  compactor 
produces  maximum  swell  pressures  nearly  equal  to  those  at  OMC.   The 
other  three  compactor  types  develop  maximum  swell  pressures  which 
are  approximately  twice  the  magnitude  produced  at  OMC. 
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Figures  4.8  and  4.9  show  the  relationship  between  OMC  and 
the  water  content  at  which  the  maximum  swell  pressure  occurs.   Two 
distinct  relationships  can  be  seen  -  one  for  laboratory  samples  and 
one  for  field  samples.   Figure  4.9  indicates  that  maximum  swell  pres- 
sures develop  dry  of  OMC,  but  nearer  to  it  for  field  compacted  samples 
than  for  laboratory  compacted  samples.   In  combination  with  Figure 
4.10,  these  graphs  also  indicate  that  maximum  swell  pressure  develops 
nearer  to  OMC  as  input  energy  increases. 

Figure  4.11  reveals  four  distinct  relationships  for  dry  density 
as  functions  of  water  content  at  maximum  swell  pressure;  these  rela- 
tionships roughly  parallel  the  zero  air  voids  curve  (100%  S  curve). 
However,  when  dry  density  at  maximum  swell  pressure  is  replotted  versus 
its  respective  OMC,  the  points  form  a  single  relationship  which  also 
parallels  the  zero  air  voids  curve  (see  Figure  4.12). 

Figure  4.13  shows  the  relationship  between  swell  pressure 
at  OMC  and  energy  ratio,  while  Figure  4.14  relates  OMC  to  energy  ratio. 
It  is  interesting  to  note  that  these  relationships  are  nearly  linear 
for  each  compactor  type  when  plotted  on  semi- logarithmic  axes. 

It  would  be  of  interest  to  determine  which  of  these  ten  rela- 
tionships are  similar  (or  identical)  for  other  compaction  equipment 
and  soil  types.   If  similarities  exist,  they  could  provide  a  simple 
basis  for  forecasting  certain  compaction  characteristics  and  swell 
tendencies.   Future  investigations  might  also  disclose  the  mechanisms 
which  produce  the  relationships  shown  in  Figures  4.5  through  4.14. 
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Figure  4.6  Maximum  Swell  Pressure  Vs.  Optimum  Moisture  Content 
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Figure  4.7   Maximum  Swell  Pressure  Vs.  Swell  Pressure  at  OMC 
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Figure  4.10  Water  Content  Deviation  from  OMC  at  SP    Vs.  Enerqv  Ratio 
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Figure  4.11   Dry  Density  at  SP    Vs.  Water  Content  at  SP 
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Figure  4.12   Dry  Density  at  SP    Vs.  Opt_mum  Moisture  Content 
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Prediction  of  Field  Response  from 
Laboratory  Compacted  Samples 

One  of  the  major  goals  of  this  study  is  the  prediction  of 
field  swell  pressure  response  from  laboratory  compacted  samples.   Sev- 
eral methods  of  prediction  exist;  however,  not  all  appear  to  be  appli- 
cable to  this  research.   The  following  is  a  discussion  of  these  methods 
and  their  potential  for  use  in  this  work. 

An  ideal  situation  would  occur  if  laboratory  and  field  com- 
pacted samples  produced  identical  response  curves.   This  is  not  the 
case,  as  evidenced  by  Figures  4.15  and  4.16.   These  figures  represent 
the  laboratory  impact  response  curves  plotted  with  those  for  the 
Caterpillar  and  Rascal  compactors,  respectively.   Similar  results 
occur  when  the  kneading  curves  are  substituted  for  the  impact  curves. 
Price  (1978)  has  also  shown  that  lab  and  field  responses  are  not  iden- 
tical . 

Figures  4.15  and  4.16  also  show  that  direct  graphical  compari- 
sons cannot  be  done.   Instead,  one  family  of  curves  must  be  trans- 
formed graphically  to  simulate  the  other  family  of  curves  presented 
in  the  graph;  direct  comparisons  can  then  be  made.   Another  method 
of  transformation  requires  a  mathematical  equation  or  set  of  equations 
that  correlate  the  field  results  with  laboratory  responses.   Such 
a   transformation  would  be  complicated  and  difficult  to  perform. 

If  laboratory  and  field  compaction  produced  nearly  identical 
sets  of  compaction  curves,  then  the  laboratory  values  could  be  directly 
substituted  into  the  regression  equations  for  field  swell  pressure. 
Since  these  sets  of  curves  are  dissimilar,  the  laboratory  values  would 
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require  transformation  before  being  used  in  the  field  response  equa- 
tions.  This  would  also  be  a  complicated  mathematical  process. 

The  approach  used  in  this  study  is  based  upon  relationships 
that  were  found  between  laboratory  and  field  compaction  results  (viz. 
maximum  dry  densities  and  optimum  moisture  contents) .   These  relation- 
ships have  been  obtained  from  the  compaction  curves  shown  in  Figures 
2.3,  2.4,  2.11,  and  2.12  and  are  presented  in  the  "Application  of 
Results".   They  are  used  to  predict  the  field  maximum  dry  density 
and  optimum  moisture  content  from  the  results  of  Standard  Proctor 
laboratory  compaction  performed  upon  the  soil.   It  is  assumed  that 
these  relationships  are  also  valid  for  other  soils  in  the  A-6  classi- 
fication. 

The  regression  models  for  field  dry  density  are  used  to  pre- 
dict the  field  compaction  curves.   In  some  cases,  the  models  must 
be  adjusted  to  correspond  with  the  predicted  maximum  dry  density. 
Because  these  adjustments   are  based  upon  relationships  developed 
for  the  compaction  curves  (not  for  the  regression  mode  Is)  ,  discrepan- 
cies in  y    might  occur  near  OMC;  judgment  should  be  used  in  these 
cases. 

The  predicted  values  of  y    and  OMC  can  then  be  utilized 
to  write  a  compaction  specification  which  will  predict  and  control 
the  magnitude  and  variability  of  swell  pressure  that  subsequently 
develop  in  the  field.   If  the  compaction  specification  has  previously 
been  written  for  a  project,  the  predicted  values  of  y         and  OMC  can 
be  used  for  inspection  testing  and  also  to  predict  the  swell  pressure 
and  its  variability  that  will  develop  for  the  compaction  procedure 
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specified.   The  complete  prediction  procedure  is  described  in  the 
following  chapter. 
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5  -  APPLICATION  OF  RESULTS 

Design  Engineering 

This  study  presents  evidence  that  large  swell  pressures  can 
be  induced  by  field  compaction  methods.   These  tendencies  to  swell 
can  reduce  the  strength  of  the  compacted  soil,  as  well  as  cause  damage 
to  structures  placed  on  or  in  the  fill.   To  ensure  proper  performance 
of  a  compacted  fill  and  the  structures  it  supports,  strength,  compres- 
sibility, and  swelling  characteristics  must  be  accounted  for  in  the 
engineering  design.   If  the  borrow  soil  is  located  and  identified 
well  in  advance  of  its  placement,  proper  compaction  specifications 
can  be  prepared  to  account  for  these  three  factors  by  use  of  such 
relations  as  are  presented  in  this  report. 

This  section  is  presented  in  order  to  illustrate  how  the  rela- 
tionships developed  in  this  study  can  be  used  to  create  an  appropriate 
specification  to  limit  the  swell  pressures  induced  by  compaction. 
These  relationships  can  be  utilized  when  an  A-6  soil  (i.e.,  one  similar 
to  the  St.  Croix  clay  of  this  study}  is  compacted  in  the  field  by 
one  of  the  compactors  used  in  this  study.   To  begin  the  specification 
process,  the  engineer  must  decide  what  magnitude  of  maximum  swell 
pressure  can  be  tolerated  in  the  design;  this  decision  is  a  judgment 
based  upon  the  requirements  of  the  project  under  consideration. 
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The  rest  of  the  procedure  is  best  presented  in  the  form  of 
the  following  example  and  commentary.   Let  us  assume  that  a  maximum 
swell  pressure  of  20  psi  (138  kPa)  can  be  tolerated  in  the  design. 
As  an  upper  permissible  bound,  this  quantity  is  represented  by  the 
expected  swell  pressure  (mean  value)  plus  the  expected  variability 
in  swell  pressure  that  will  occur  with  a  particular  set  of  compaction 
conditions  (i.e.  ,  compactor  type,  E  ,  Aw,  and  half-range  variability 
of  Aw).   Figures  5.1  thru  5.4  can  be  used  directly  to  obtain  appropri- 
ate combinations  of  compaction  conditions  that,  if  used,  will  assure 
a  maximum  swell  pressure  of  20  psi  (138  kPa) . 

A  number  of  figures  are  presented  herein  to  convey  the  various 
data  developed  in  this  study.   Figure  5.1  is  developed  for  the  Cater- 
pillar compactor  over  all  energy  levels  and  Figures  5.2  thru  5.4  repre- 
sent the  three  energy  levels  for  the  Rascal,  roller.   The  expected 
maximum  swell  pressures,  for  different  values  of  the  half-range  vari- 
ability of  Aw,  are  shown  by  the  solid  curves  in  the  four  figures; 
again,  the  expected  maximum  SP  represents  the  sum  of  the  expected 
mean  SP  and  the  expected  variability  of  SP  for  the  given  conditions 
of  compaction.   The  dashed  curves  indicate  the  expected  mean  swell 
pressure  and  are  drawn  within  the  limits  of  the  gathered  data.   These 
limits  should  be  observed  when  using  all  curves  in  these  figures, 
as  values  outside  these  Aw  bounds  are  an  extrapolation  of  the  data 
and  are  subject  to  error.   Horizontal  dashed  lines  appear  in  Figures 
5.1  thru  5.3  wet  of  OMC.   These  are  believed  to  be  reasonable  limiting 
values  of  expected  swell  pressures,  rather  than  the  large  increases 
shown  for  the  wet  side.   Reasons  for  indicating  that  these  increases 
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are  not  to  be  expected  in  the  field  are  given  in  the  "Discussion  of 
Results"  section. 

Returning  to  the  example,  results  obtained  for  the  Caterpillar 
compactor  (see  Figure  5.1)  indicate  that  this  maximum  swell  pressure, 
20  psi  (138  kPa) ,  will  not  be  exceeded  if  the  following  conditions 
are  met:   (1)  the  number  of  passes  ranges  between  four  and  sixteen 
(i.e.,  1  <  E  <  4) ;  (2)  the  soil  is  compacted  at  least  3.3%  wet  of 
OMC  with  a  half-range  variability  of  Aw   less  than  0.2%;  or  (3)  the 
soil  is  compacted  at  least  3.5%  wet  of  OMC  with  a  half-range  variabil- 
ity of   Aw   less  than  2.0%.   Other  possible  combinations  also  exist 
for  intermediate  values  of   Aw  variability  which  can  be  located  by 
interpolation.   Dry  of  OMC,  maximum  swell  pressures  less  than  20  psi 
(138  kPa)  can  be  achieved  only  outside  the  limits  of  the  gathered 
data;  therefore,  no  attempt  is  made  here  to  denote  these  compaction 
conditions. 

For  the  Rascal  compactor,  the  maximum  swell  pressure  will 
not  be  exceeded  if  four  passes  of  the  compactor  are  used  between  Aw 
of  0.0%  to  8.0%  with  a  half-range  variability  of  Aw   less  than  0.2%. 
This  may  be  expressed  as  (1,  0.0,  8.0,  0.2),  i.e.,  (energy  ratio 
=1=4  passes;  0.0  =  minimum  value  of  Aw;  8.0  =  maximum  value  of 
Aw;  0.2  =  maximum  half-range  of  variability  of  Aw).   Using  this  man- 
ner of  expression,  and  referring  to  Figure  5.2,  the  maximum  swell 
pressure  will  not  be  exceeded  with  (1,  0.3,  8.0,  0.6),  (1,  0.5,  8.0, 
1.0),  (1,  0.7,  8.0,  1.4),  (1,  0.8,  8.0,  2.0).   Figure  5.3  shows  that 
the  20  psi  (138  kPa)  maximum  swell  pressure  will  always  be  exceeded 
if  this  soil  is  compacted  at  the  B  energy  level  (i.e.,  8  passes). 
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From  Figure  5.4  it  can  be  seen  that  the  maximum  swell  pressure  will 
not  be  exceeded  for  the  following  combinations  of  variables  (4,  -2.0, 
1.3,  0.2) ,  (4,  -0.7,  0.5,  1.0) . 

In  the  development  of  the  project's  specifications  for  compac- 
tion the  above  will  control  the  swell  pressure.   However,  strength, 
compressibility,  and  other  requirements  must  also  be  considered  as 
the  engineer  develops  "trade-offs"  for  the  in-service  fie  Id  behavior. 

Since  the  results  of  this  study  are  based  on  water  content  de- 
viation from  OMC,  it  is  important  to  know  where  optimum  moisture  con- 
tent and  maximum  dry  density  will  occur  in  the  field.  If  economies  (of 
time  and  money)  do  not  permit  the  construction  of  field  test  sections 
to  determine  OMC  and  y  directly,  the  following  procedure  can  be 

used  to  predict  these  quantities.   First,  a  laboratory  compaction 
curve  is  generated  for  the  soil  with  the  Standard  Proctor  energy  level; 
this  can  be  done  either  by  impact  or  kneading  compaction.   Compare 

these  values  of  OMC  and  v,    with  those  in  Table  5.1.   If  the  values 

dmax 

are  equal  or  nearly  so,  field  values  of  OMC  and   v ^     can  be  taken 
n  *  'dmax 

directly  from  Table  5.1. 

If  the  generated  y^    is  different,  then  Figure  5.5  can  be 
dmax 

used  to  adjust  the  tabular  value  of  maximum  dry  density  to  predict 
that  which  should  be  anticipated  in  the  field.   Figure  5.5  shows  the 
ratios  of  maximum  dry  densities  produced  in  this  study  (field  nuclear 
gage  to  laboratory  Standard  Proctor  values)  versus  field  energy  ratios. 
To  use  Figure  5.5  choose  the  curve  which  corresponds  to  the  compactor 
types  being  used  (e.g.,  impact  and  Rascal),  determine  the  energy  ratio 
at  which  the  field  compactor  will  be  used  (e.g.,  4,  which  means  16 
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Table  5.1   Optimum  Moisture  Contents  and  Maximum  Dry  Densities  Used 
for  Predictions 


Compactor 
Type 

Location  of 
Measurement 

Energy 
Level 

Energy 
Ratio* 

OMC 
(%) 

v 
dmax 

(pcf) 

Impact 

Compaction 
Mold 

B 
(Standard 
Proctor) 

1.67 

16.6 

113.8 

Kneading 

Compaction 
Mold 

B 
(Standard 
Proctor) 

1.75 

16.5 

114.9 

Caterpillar 

Field 

A 

1 

15.7 

105.5 

Nuclear 

B 

2 

15.6 

110.9 

Gage 

C 

4 

15.5 

111.5 

Rascal 

Field 

ft 

1 

16.2 

103.7 

Nuclear 

B 

2 

15.4 

107.2 

Gage 

C 

4 

15.0 

110.6 

To  obtain  the  number  of  passes  for  field  compactors,  multiply 
the  energy  ratio  by  4. 
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passes) ,  and  locate  their  point  of  intersection.   Obtain  the  ordinate 

of  this  point  (ratio  of  y  ;    e.g.,  0.972)  and  multiply  the  generated 

Standard  Proctor  maximum  dry  density  by  this  value.   This  will  yield 

the  predicted  value  of  maximum  dry  density  in  the  field  (y    ) • 

pmax 

To  obtain  the  other  values  of  dry  density  for  the  predicted 
field  compaction  curve  the  following  procedure  can  be  used.   First, 
pick  the  appropriate  field  regression  equation  (see  Table  3.3  or  Fig- 
ures 4.1  thru  4.4)  and  insert  appropriate  values  of   Aw  for  the  E 

R 

under  consideration.   Next,  multiply  each  y        so  obtained  by  the 
ratio : 

field  v     for  the  soil  under  consideration 
pmax 

appropriate  field  y,  from  Table  5.1 

dmax 

Then,  one  more  adjustment  of  this  y        is  necessary  because  the  regres- 
sion equations  for  field  dry  density  were  developed  for  test  specimen 
values.   This  final  adjustment  is  to  multiply  by  0.93  to  yield  the 
proper  predicted  field  dry  density  (y  )•   The  coefficient  of  0.93  is 
the  average  ratio  of: 

v     (field  nuclear  gage  values) 
'dmax . 

Y     (test  specimen  values) 
'dmax 

The  underlying  assumption  is  that  the  field  compaction  curves  are  of 
the  same  general  shape  as  the  regression  equation  obtained  for  the 
test  specimens. 

If  the  generated  Standard  Proctor  OMC  is  appreciably  different 
from  that  shown  in  Table  5.1,  the  appropriate  field  OMC  can  be  pre- 
dicted in  the  following  manner.   First,  calculate  the  ratio  of  the 
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generated  Standard  Proctor  OMC  to  that  obtained  in  this  study.   Then 

multiply  the  appropriate  value  of  field  OMC  (given  in  Table  5.1)  by 

this  quantity  to  obtain  the  magnitude  of  the  predicted  field  optimum 

moisture  content  (OMC  ) . 

P 

As  a  continuation  of  the  previous  example,  let  us  assume  that 
Standard  Proctor  impact  compaction  on  the  soil  under  consideration 
produced  a  maximum  dry  density  of  116  pcf  and  OMC  of  15.0%.   The  engi- 
neer has  accepted  the  Rascal  compactor  at  an  E  =  1.0  (4  passes)  and 

R 

Aw  =  +2.0%  (see  Figure  5.2)  to  produce  the  allowable  maximum  swell 
pressure;  the  half-range  variability  of   Aw  will  be  controlled  at  1.4%. 
Using  Figure  5.5,  the  appropriate  ratio  of 

y,     (field  nuclear  gage  value) 

dmax  _  „  , 

is  0.911. 


y,     (laboratory  Standard  Proctor  value) 
dmax 


Thus,  the  predicted  field  maximum  dry  density  is: 


y     =  116  pcf  x  0.911  =  105.7  pcf. 
'pmax  e 

To  predict  the  expected  field  dry  density,  insert  the  appropriate   Aw 
and  E   into  the  proper  regression  equation  (see  Table  3.3,  line  6): 

R 

Yd   =   -1.378   x   2.0   +    1.044   x    1.0   +    112.0   =   110.3   pcf. 


This  value  could  also  have  been  obtained  from  Figure  4.2.   Next,  multi- 
ply this  quantity  by  the  appropriate  ratio  of 

predicted  field  y,     (i.e.,  y     ) 
dmax pmax 


appropriate  field  Y,     from  Table  5.1 

dmax 


,  ,  .  ,  .   105.7  pcf  ,    .    .,  ,   „  __ 
(which  is  ,--  =-E—=   here)  and  by  0.93: 
103.7  pcf 


134 


Yp  =  110.3  pcf  x  ^3^  ggf  x  0.93  =  104.5  pcf. 


To  predict  the  field  OMC,  multiply  the  appropriate  field  OMC 
in  Table  5.1  (16.2%)  by  the  ratio 


generated  Standard  Proctor  OMC 15.0% 

Standard  Proctor  OMC  from  Table  5.1    (l-e-'  1e.6%) 


Thus,  the  predicted  field  OMC  is: 


0MC   =  ^'?!  x  16.2%  =  14.6%. 
p    16.6% 


For  the  soil  under  consideration,  the  oredicted  values  of  field  OMC  , 

P 

Y     ,  and  y   at  Aw  =  2.0%  are  14.6%,  105.7  pcf,  and  104.5  pcf,  respec- 
pmax       p  r 

tively . 

Application  of  the  results  obtained  in  this  study  is  dependent 

upon  the  OMC  and  v,    which  will  be  achieved  in  the  field.   Predictinq 
dmax  ^ 

the  field  OMC  by  methods  described  above  allows  the  engineer  to  write 
a  proper  compaction  specification;  i.e.,  once  the  allowable  average 
or  maximum  swell  pressure  and  field  OMC  are  determined,  Figrues  5.1 
thru  5.4  can  be  used  to  specify  the  allowable  absolute  values  of  field 
moisture  content. 

Quality  Control 
There  are  cases  where  the  borrow  soil  to  be  used  for  compac- 
tion is  not  located  and  identified  well  in  advance  of  construction. 
For  these  cases,  the  compaction  specification  may  not  contain  the  neces- 
sary provisions  that  will  control  the  induced  swell  pressures  to  a 
tolerable  level.   However,  it  is  still  desirable  to  be  able  to  predict 
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the  maximum  swell  pressures  for  the  resulting  product.   This  section 
outlines  a  method  for  predicting  swell  pressures  from  inspection  test 
results  of  compaction  projects  in  progress. 

Figures  5.1  thru  5.4  could  be  adapted  for  use  by  the  quality 
control  engineer,  but  they  do  not  provide  all  of  the  necessary  informa- 
tion.  Computer  printouts  generated  by  the  programs  in  Appendix  B  would 
be  better  suited  for  this  purpose.   Table  5.2  is  a  sample  of  the  output 
produced  by  the  Caterpillar  program;  it  can  be  used  for  quality  assur- 
ance testing  and  for  prediction  of  the  expected  swell  pressure  magnitude 
and  variability. 

Before  this  table  can  be  used,  the  field  optimum  moisture 

content  and  maximum  dry  density  must  be  determined.   If  this  cannot 

be  established  from  the  on-going  field  operation,  these  values  can 

be  predicted  by  methods  described  in  the  "Design  Engineering"  section; 

i.e.,  a  laboratory  Standard  Proctor  compaction  test  is  performed,  and 

Y     and  OMC  are  generated  using  the  orooer  charts, 
pmax        p 

Next,  obtain  the  number  of  passes  of  the  compaction  equipment 
and  the  fill  water  content  and  dry  density  for  the  lift  being  examined. 
A  minimum  of  five  to  seven  samples  for  w   and  y        must  be  taken  to 
accommodate  the  approach  used  in  this  study.   The  average  water  content 
and  dry  density  are  then  calculated  for  the  group  of  samples.   Also 
calculate  the  expected  variabilities  in  the  water  content  and  dry  densi- 
ty; these  variabilities  are  equal  to  the  half-range  variation  for  each 
parameter,  as  calculated  from  the  group  of  samples.   Then  Aw  is  calcu- 
lated as  the  difference  between  the  average  water  content  and  field 
OMC  previously  determined. 
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With  the  above  data,  proceed  into  the  appropriate  section 
of  the  table  using  the  measured  half-range  of  the  water  content,  V(Aw), 
as  the  variability  measurement.   This  value  appears  in  the  left-hand 
column  (1)  of  the  table;  Table  5.2  is  based  upon  the  average  water 
content  variability  found  for  the  St.  Croix  test  pad  test  specimens 
(1.4%).   Next,  locate  the  proper  value  of  Aw  in  column  two,  and  the 
appropriate  number  of  passes  in  the  third  column  from  the  left;  energy 
ratio  (E  )  is  presented  in  column  four. 

R 

On  the  appropriate  line  for  the  combination  of  these  compaction 
variables,  column  five  contains  the  expected  dry  density  (y,  in  pcf) 
for  test  specimens.   Contents  of  column  five  are  used  in  the  variabil- 
ity analyses  and  are  consistently  higher  than  values  measured  in  the 
field  by  nuclear  devices.   The  predicted  field  dry  density  (y   in  pcf) 
is  given  in  column  six  and  should  be  used  for  comparisons  with  inspec- 
tion test  results.   If  the  appropriate  y^     in  Table  5.1  (lab  or 

dmax 

field,  depending  on  how  it  is  determined  for  the  borrow  material)  is 

different  from  that  achieved  for  the  borrow  soil,  the  values  in  column 

six  should  be  adjusted  by  methods  given  in  the  "Design  Engineering" 

section  to  provide  a  more  reliable  basis  for  comparisons  with  inspection 

tests  results.   Column  seven  shows  the  expected  dry  density  variability 

(V(y  )  in  pcf) ;  this  is  appropriate  for  use  with  both  measures  of  dry 

density  presented  in  the  table. 

The  expected  swell  pressure  (SP  in  psi)  appears  in  the  eighth 

column,  and  its  expected  variability  (V(SP)  in  psi)  is  given  in  column 

nine.   The  last  column  in  the  table  presents  the  expected  maximum  swell 

pressure  (SP    in  psi),  which  is  the  sum  of  the  entries  in  columns 
max 
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eight  and  nine.   The  values  in  these  last  three  columns  should  be  com- 
pared with  the  wet  side  limiting  values  (horizontal  dashed  lines)  pre- 
sented in  Figures  5.1  thru  5.4.   The  limiting  values  should  be  used, 
in  their  applicable  ranges  of  Aw,  to  predict  the  field  swell  pressures. 

Of  particular  interest  are  columns  six,  seven,  eight,  and 
ten.   Columns  six  and  seven  can  be  used  as  a  check  for  the  compaction 
process,  i.e. ,  as  a  check  against  gross  errors  being  made  in  establish- 
ing the  compaction  variables.   If  large  discrepancies  exist,  the  pre- 
dictions will  be  questionable.   Columns  eight  and  ten  represent  the 
best  estimate  of  the  swell  pressure  and  its  variability  for  the  com- 
pacted soil.   If  the  measured  dry  density  variability  is  greater  than 
the  value  shown  in  column  seven,  the  values  in  columns  nine  and  ten 
may  be  in  slight  error.   Again,  swell  pressure  variability  is  relatively 
insensitive  to  small  changes  in  dry  density  or  its  associated  variabil- 
ity. 

The  tabular  form  of  the  regression  results  provides  the  quality 
control  engineer  with  a  convenient  means  of  evaluating  field  inspection 
tests,  predicting  subsequent  swell  behavior,  and  determining  when  cor- 
rective action  is  required  in  the  field.   It  can  also  be  used  to  assess 
the  application  of  these  results  to  other  projects  and  soil  types. 

The  computer  programs  in  Appendix  B  generate  tables  for  a 
wide  variety  of  compaction  conditions,  of  which  Table  5.2  is  a  small 
sample.   These  programs  can  easily  be  modified  to  assess  conditions 
not  presented  herein.   However,  the  results  should  be  used  only  within 
the  limits  of  the  data  (compactor  type,  Aw,  E  )  gathered  in  this  study. 
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6  -  CONCLUSIONS 

For  the  residual  medium  plastic  St.  Croix  clay  and  methods  of 
compaction  used  in  this  study,  the  following  conclusions  can  be  drawn: 

1.  Soil  compacted  by  the  Caterpillar  compactor  exhibits  the 
following  relationships : 

a.  The  magnitude  of  dry  density  is  governed  by  the  water  con- 
tent deviation  from  OMC  (Aw)  for  conditions  dry  of  OMC. 

Wet  of  OMC,  both  Aw  and  energy  ratio    (E  )   contribute 

R 

toward  the  resulting  magnitude. 

b.  Swell  pressure  magnitude  is  defined  by   Aw  dry  of  OMC, 

2 
while  both   Aw   and   Aw   are  important  wet  of  OMC. 

2.  For  soil  compacted  by  the  Rascal  compactor,  the  following  re- 
lationships exist: 

a.  Dry  density  magnitude  is  most  influenced  by  Aw  and  E 
both  wet  and  dry  of  OMC. 

b.  Swell  pressure  magnitude  is  divided  into  two  subsets;  the 

A   and   B   energy  levels  form  one  subset  and  the  C   level 

2 
comprises  the  other.   For  the  former  subset,   Aw  ,  E  ,  and 

R 

the  interaction  term  AWYH  produce  the  best  relationships, 

2 
even  though  Aw   and  E   are  not  statistically  significant. 

No  significant  relation  exists  for  the   C   energy  level, 

but  Aw  and  AwY,   generate  reasonable  magnitudes. 
a 
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3.  For  soil  compacted  by  the  laboratory  impact  and  kneading 
methods  used  in  this  study,  the  following  relationships  are  obtained: 

a.  The  two  compactor  types  produce  similar  compaction  curves, 
with  Aw  and  E  being  the  most  important  descriptors  of 
dry  density  magnitude  both  wet  and  dry  of  OMC. 

b.  The  variables  contributing  most  to  the  resultant  swell 

pressure  magnitude  are  Aw,  E  ,  and  Awy  .   Swell  pressure 

R  a 

response  curves  are  also  similar  for  the  two  compactor 
types,  with  impact  samples  generating  slightly  larger  swell 
pressures. 

4.  The  following  relationships  occur  for  the  variability  of  dry 
density  and  swell  pressure  for  field  compacted  samples: 

a.  The  magnitude  of  SP  variability  is  large  and  is  generally 
higher  for  the  Rascal  compactor. 

b.  The  magnitudes  of  the   y    and   SP   variabilities  of  both 
compactors  are  reduced  if  the  water  content  variability  is 
reduced. 

c.  Due  to  the  inherent  variability  of  a  compacted  soil  mass, 
the  average  of  a  number  of  samples  (a  minimum  of  5  to  7) 
must  be  used  to  compute  the  measured  values  of  field  water 
content  and  dry  density. 

5.  Field  test  specimens  exhibit  higher  and  less  variable  values  of 
dry  density  than  those  measured  in  the  field  by  nuclear  devices.   This 
is  probably  caused  by  the  occurrence  of  drying  during  storage  as  well 
as  differences  in  the  location  and  volume  of  the  material  within  the 
measurement,  rather  than  by  trimming  and  testing  procedures. 
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6.   A  method  has  been  devised  that  is  believed  to  reasonably  pre- 
dict the  field  optimum  moisture  content,  dry  density,  and  swell 
pressure  response  from  laboratory  compaction  data. 
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7  -  RECOMMENDATIONS  FOR  FUTURE  STUDIES 

1.  Verify  the  field  dry  density  prediction  relations  presented 
in  this  study.   This  should  be  done  for  St.  Croix  clay  as  well  as 
other  A-6  soils  at  future  compaction  project  sites. 

2.  Verify  the  field  swell  pressure  prediction  relations  pre- 
sented herein  for  this  and  other  A-6  soils.   Predictions  could  be 
made  for  future  embankments  or  those  now  in-service;  field  samples 
could  be  obtained  and  tested  to  verify  the  relations.   This  could 
also  include  the  development  of  methods  for  measuring  the  in  situ 
field  swell  pressure. 

3.  Develop  prediction  relations  for  field  dry  density,  swell 
pressure,  and  other  mechanical  properties   for  different  compactors 
and  soil  types. 

4.  Develop  better  methods  to  measure  and  control  the  sources 
of  variability  in  field  compaction  processes.   Application  of  these 
results  is  necessary  to  improve  the  overall  uniformity  of  a  compacted 
soil  mass. 

5.  Determine  if  the  associated  relationships  presented  in  this 
study  are  valid  for  other  soils  and  compactor  types. 

6.  Develop  and  standardize  swell  pressure  tests  that  will  produce 
the  specific  information  desired  for  a  particular  project.   These 
include : 
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a.  For  CVS  tests,  develop  a  rigid  testing  appratus  that 
will  ensure  constant  volume  and  provide  a  simpler  testing 
procedure  than  that  required  for  a  consolidometer. 

b.  Perform  swell  tests  in  a  triaxial  cell  and  determine  the 
loss  in  shear  strength  due  to  swelling  tendencies.   This 
apparatus  can  also  be  used  to  measure  three-dimensional 
(or  lateral)  swell  pressures  and  volume  changes. 

c.  Using  a  consolidometer,  apply  various  surcharge  loads  on 
the  soil  samples  before  inundation.   This  may  better  simu- 
late field  loading  conditions  and  allow  the  appropriate 
amount  of  volume  change  to  occur,  and  can  be  used  to  re- 
late volume  change  to  the  reduction  in  swell  pressure  as 
compared  to  the  CVS  condition. 

7.  Determine  the  effects  of  drying  on  susequent  swell  behavior, 
and  compare  with  results  from  samples  that  remain  at  their  as-compacted 
water  content. 

8.  Determine  if  there  is  a  relationship  between  the  as-compacted 
pore  size  distribution  and  the  swell  pressure  that  develops.   Also  com- 
pare pore  size  data  taken  from  samples  before  and  after  swelling  has 
occurred. 
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Appendix  A 
Presentation  of  Field  and  Laboratory  Compacted  Sample  Data 


FIELD  COMPACTED  SAMPLE  DATA 
-CATERPILLAR  COMPACTOR- 
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TEST 

TEST 

SPECIMEN 

SPECIMEN 

SAMPLE 

ENERGY 

LIQUID 

WATER 

DRY 

SWELL 

NO. 

RATIO 

LIMIT 

CONTENT 

DENSITY 

PRESSURE 

(PERCENT) 

(PERCENT) 

(PCF) 

(PSI) 

C1A3 

1.0 

45. a 

15. G 

111.7 

41.3 

C1A4 

1.0 

42.7 

15.7 

114.3 

30.9 

C1A6 

1.0 

39.1 

1G.0 

110. G 

17.3 

C2A4 

1.0 

32. S 

18.1 

111.7 

8.0 

C2A5 

1.0 

41.0 

14.0 

111.2 

19.1 

C3A2 

1.0 

42.5 

14.7 

111.5 

48.0 

C3A4 

1.0 

39. S 

15.2 

116.0 

20.2 

C3A5 

1.0 

32.9 

14.7 

11G.2 

13.8 

C4A1 

1.0 

38.3 

20.0 

10G.3 

5.6 

C4A2 

1.0 

43.9 

17.7 

110.5 

15.5 

C4A3 

1.0 

39.7 

19.1 

107.9 

6.6 

C5A1 

1.0 

4S.0 

17.0 

112.3 

41.2 

C5A2 

1.0 

44.9 

1G.0 

111.1 

26.3 

C5AG 

1.0 

41.9 

17.4 

112.2 

14.3 

C1B1 

2.0 

37.1 

13.4 

119.7 

45.6 

C1B2 

2.0 

42.8 

13.4 

116.9 

49.2 

C1B4 

2.0 

33.0 

13.5 

120.2 

49.2 

C1B5 

2.0 

42.0 

14.2 

117.2 

42.6 

C1BG 

2.0 

43. G 

13.4 

115.7 

36.2 

C2B2 

2.0 

38.3 

15.0 

115.8 

33.8 

C2B3 

2.0 

38.8 

14.5 

116.4 

39.9 

C2BS 

2.0 

32.5 

12.7 

119.4 

32.0 

C3B4 

2.0 

41.0 

12.0 

10G.9 

15.4 

C4B1 

2.0 

42.8 

18.5 

111.1 

12.5 

C4B4 

2.0 

41.0 

17.3 

111.5 

13.0 

C4B5 

2.0 

44.0 

15.4 

113.3 

30.6 

C5B2 

2.0 

4G.1 

15.7 

115.3 

30.0 

C5B4 

2.0 

42.5 

18. G 

109.1 

14.1 

C5B5 

2.0 

42.1 

1G.7 

109.2 

15.4 

C1C1 

4.0 

39. G 

13.4 

119.2 

53.2 

C1C2 

4.0 

42.3 

14. G 

118.1 

43.3 

C1C5 

4.0 

38.4 

12.7 

115. G 

17.3 

C2C1 

4.0 

43.0 

14.1 

117.1 

48.0 

C2C5 

4.0 

45.5 

13.1 

113.3 

24.5 

C3C2 

4.0 

39.7 

20.4 

107.0 

6.6 

C3C4 

4.0 

30.0 

15.7 

117.0 

7.2 

C3C5 

4.0 

3S.5 

14.0 

114.3 

19.7 

C4C3 

4.0 

44.5 

19.0 

108.5 

8.2 

C4C5 

4.0 

43.5 

1G.G 

112.8 

17.1 

C4CG 

4.0 

41.9 

15.4 

113.9 

24.2 

C5C1 

4.0 

45.9 

13.5 

108.3 

13.3 

C5C3 

4.0 

43.0 

23.7 

101.4 

13.3 

C5C4 

4.0 

33. G 

17.7 

111.4 

12.0 

C5C5 

4.0 

41.0 

1G.2 

112.4 

17.5 

FIELD  COMPACTED  SAMPLE  DATA 
-RASCAL  COMPACTOR- 
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TEST 

TEST 

SPECIMEN 

SPECIMEN 

SAMPLE 

ENERGY 

LIQUID 

HATER 

DRY 

SWELL 

MO. 

RATIO 

LIMIT 

CONTENT 

DENSITY 

PRESSURE 

( PERCENT) 

(PERCENT) 

(PCF) 

(PSI) 

R1A1 

1.0 

3G.S 

13.7 

111.6 

22.2 

R1A3 

1.0 

37.0 

14.1 

110.2 

11.4 

R1A4 

1.0 

37.2 

13.5 

110.9 

19.9 

RlftS 

1.0 

3G.2 

14.4 

110. G 

15.3 

R2A1 

1.0 

34.0 

14.1 

112. G 

12.5 

R2A5 

1.0 

3G.1 

15.4 

109. G 

7.0 

R2AS 

1.0 

35. G 

15.0 

111.4 

11.3 

R3A2 

1.0 

40.8 

15.9 

10G.0 

5.4 

R3A5 

1.0 

38.0 

1G.0 

111.6 

12.0 

R3AB 

1.0 

39.8 

14.4 

113.2 

33.6 

R4A2 

1.0 

51.8 

19.5 

108.1 

15.5 

R4A3 

1.0 

34.0 

15.9 

113.3 

2.6 

R4A4 

1.0 

43.5 

19.3 

108. G 

17.5 

R4A5 

1.0 

45. G 

20.9 

104.3 

7.5 

R4A6 

i.O 

43.4 

17.9 

108.7 

16.8 

R5A1 

1.0 

44. G 

20.2 

107.2 

5.2 

R5A3 

1.0 

3G.4 

21.0 

10G.O 

8.7 

R5A5 

1.0 

43.0 

20.5 

105.3 

24.0 

R5AG 

1.0 

41.2 

19.4 

109.7 

4.8 

R1B1 

2.0 

33. G 

14.2 

113.3 

25.3 

R1B2 

2.0 

3G.4 

13. G 

113.0 

24.7 

R1B4 

2.0 

37.0 

1G.3 

112.5 

12.5 

R2B1 

2.0 

38.4 

15.7 

114.1 

10.7 

R2B3 

2.0 

3G.5 

14.9 

116.9 

16.0 

R2BS 

2.0 

34.3 

15.3 

112.8 

11.7 

R3B4 

2.0 

39.4 

13.8 

115.8 

37.3 

R3B5 

2.0 

39.2 

14.0 

113.6 

26.0 

R4B2 

2.0 

44.0 

1G.1 

112.9 

9.4 

R4B3 

2.0 

42.5 

17.0 

111.8 

13.3 

R4B4 

2.0 

45.2 

18.4 

111.4 

19.5 

R5B2 

2.0 

41.4 

21.8 

103.0 

5.1 

R5B4 

2.0 

42.4 

18.9 

109.7 

15.4 

R5BG 

2.0 

4G.G 

18.0 

110. G 

17.3 

R1C3 

4.0 

34.8 

15.8 

114.8 

11.5 

R1C4 

4.0 

37.0 

15.0 

117.5 

17.0 

R1CS 

4.0 

34.2 

11.8 

111.4 

19.8 

R2C1 

4.0 

31.1 

12.9 

119. 6 

14.6 

R2C2 

4.0 

37.4 

12. G 

114.1 

11.6 

R2C3 

4.0 

33.0 

14.3 

114.5 

12.8 

R2C4 

4.0 

37.3 

14.1 

109.0 

6.9 

R2C5 

4.0 

37.2 

14.0 

117.3 

15.7 

R3C3 

4.0 

39.8 

18.9 

110.3 

26.9 

R3CG 

4.0 

37.0 

14.7 

119.0 

17.5 

R4C1 

4.0 

47.4 

17. G 

112.4 

19.4 

R4C4 

4.0 

40.7 

12.3 

114.8 

11.4 

R4C5 

4.0 

40.5 

18.8 

109.4 

8.8 

R5C3 

4.0 

42.4 

17.2 

110.3 

15.3 

R5C4 

4.0 

34.6 

22.4 

102.9 

11.6 

R5CG 

4.0 

4G.8 

23.1 

101.9 

5.4 
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LABORATORY  COMPACTED   SAMPLE   DATA 
-IMPACT    COMPACTOR- 


TEST 

TEST 

COMPACTION 

COMPACTION 

SPECIMEN 

SPECIMEN 

MOLD 

MOLD 

SAMPLE 

ENERGY 

WATER 

DRY 

HATER 

DRY 

SU'ELL 

ISO. 

RATIO 

CONTENT 

DENSITY 

CONTENT 

DENSITY 

PRESSURI 

(PERCENT) 

CPCF) 

(PERCENT) 

(PCD 

(PSD 

I1A1 

1.00 

12.9 

101.0 

13.2 

101.3 

11.4 

I1A2 

1.00 

12. S 

104.3 

12.2 

103.7 

19.7 

I2A1 

1.00 

14.  G 

107. S 

15.3 

10-1.9 

13.0 

I2A2 

1.00 

13.8 

10G.5 

14.2 

105. G 

16.0 

I2A3 

1.00 

14.5 

109.7 

14.3 

10G.2 

19.7 

I3A2 

1.00 

1G.G 

111.1 

1G.4 

110.1 

14.7 

I3A1 

1.00 

17.3 

110.2 

17,8 

109.1 

9.3 

I4A1 

1.00 

18.1 

109.5 

13. G 

108.5 

8.8 

I5A1 

1.00 

19. G 

107.0 

20.7 

107.0 

4.5 

IGA1 

1.00 

-0 

-0 

22.8 

102.2 

-0 

I1B1 

l.G? 

12.3 

107.2 

13.3 

108.1 

25.3 

I1B2 

1.G7 

12.2 

109.9 

12.3 

108.9 

29.3 

I2B2 

1.67 

14.2 

112.8 

14.5 

112.1 

24.0 

I2B3 

1.G7 

14.4 

109.5 

14.1 

111.6 

20.0 

I2B1 

1.G7 

15.0 

114. G 

15.4 

112.7 

25.3 

I3B2 

1.G7 

1G.5 

113.1 

1G.G 

113.8 

17.9 

I2'B1 

1.G7 

1G.8 

111.9 

17.1 

113. G 

10.1 

I3B1 

1.G7 

17.0 

112.3 

17.4 

113.2 

10.1 

I4B1 

1.G7 

17.9 

110.5 

18.4 

111.7 

7.7 

I5B1 

1.G7 

19.4 

107. G 

20.5 

108.4 

4.9 

IO'Cl 

7.58 

8.3 

118.4 

8.1 

120. G 

102.3 

I0C1 

7.58 

9.4 

120.2 

9.8 

121. 6 

112.3 

I1C2 

7.58 

7.5 

115.7 

7.5 

119.5 

80.  G 

I2C2 

7.53 

8.8 

117.0 

8.G 

118.3 

9G.3 

I3C2 

7.58 

10.2 

120.7 

9.9 

123.3 

123.6 

IO'Cl 

7.58 

10.8 

121.2 

10.8 

123. G 

89.0 

I1*C1 

7. S3 

12.1 

122.0 

12.2 

124.9 

57.6 

I1C1 

7.58 

12.7 

122. G 

13.4 

123.7 

43.9 

I2C1 

7.58 

15.2 

11G. 7 

15.3 

120.4 

18.7 

I3C1 

7.58 

1G.6 

113.0 

17.3 

115.0 

9.3 

I4C1 

7.58 

17.8 

110.3 

18.2 

113.1 

7.5 
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LABORATORY  COMPACTED  SAMPLE  DATA 
-KNEADING  COMPACTOR- 


TEST 

TEST 

COMPACTION 

COMPACTION 

SPECIMEN 

SPECIMEN 

MOLD 

MOLD 

SAMPLE 

ENERGY 

HATER 

DRY 

WATER 

DRY 

SWELL 

NO. 

RATIO 

CONTENT 

DENSITY 

CONTENT 

DENSITY 

PRESSURE 

(PERCENT) 

(PCF) 

(PERCENT) 

(PCF) 

(PSD 

K1A1 

1.00 

12.9 

97.2 

13.1 

99.1 

10. S 

K1A2 

1.00 

12.0 

95.6 

11.9 

98.8 

7.7 

K2A1 

1.00 

14.9 

103.4 

14.3 

102.7 

12.0 

K2A2 

1.00 

14.3 

100.0 

14.2 

102.5 

10.4 

K3A2 

1.00 

15.5 

105.0 

15.4 

103.6 

14.4 

K3A1 

1.00 

16.7 

108.6 

16.8 

107.2 

10.6 

K4A1 

1.00 

18.1 

109.7 

13.5 

110.6 

7.5 

K5A1 

1.00 

19.7 

106.6 

20.1 

108.5 

4.0 

KGA1 

1.00 

-0 

-0 

21.7 

104.7 

-0 

K1B1 

1.75 

12.8 

10S.0 

12.7 

108.9 

22.4 

K1B2 

1.75 

11.9 

106.3 

12.1 

105.9 

IS. 5 

K2B2 

1.75 

14.0 

110.4 

14.0 

109.0 

24.3 

K2B3 

1.75 

14.1 

109.2 

14.1 

110.0 

20.5 

K2B1 

1.75 

15.0 

111.2 

14.3 

112.8 

18. S 

K3B1 

1.75 

IS. 7 

112.9 

IS. 7 

114.7 

10.2 

K3B2 

1.75 

15.8 

113.6 

15.9 

113.1 

18.6 

K4B1 

1.75 

18.1 

110.2 

18.2 

113.2 

8.0 

K5B1 

1.75 

19.2 

107.7 

19.1 

110.1 

3.2 

KO*Cl 

G.50 

8.1 

113.5 

8.1 

114.9 

71.5 

K0C1 

G.50 

9.8 

116.1 

9.8 

117.6 

84.2 

KO'Cl 

S.50 

11.0 

117.5 

10.9 

119.7 

78. S 

KO'C2 

S.50 

10.4 

117.7 

10.2 

119.2 

74.9 

K1C2 

S.50 

7.6 

110.1 

7.2 

114.1 

50.5 

K2C2 

S.50 

9.0 

112.2 

8.8 

113.9 

77.0 

K3C2 

S.50 

11.5 

120.4 

11.3 

120.6 

91.3 

K1C1 

S.50 

12.5 

122. 6 

12. S 

124.1 

44.7 

K2C1 

6.50 

14.4 

118.3 

14.4 

121.6 

17.1 
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Appendix  B 

Computer  Programs  for  Evaluating 

Swell  Pressure  Variability 
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An  Example  of  the  Computer  Programs'  Use 
The  computer  programs  presented  in  this  appendix  are  used  to 
evaluate  the  magnitude  and  variability  of  the  dry  density  and  swell 
pressure  for  field  compacted  samples.   The  three  programs  presented 
herein  are  for  the  Caterpillar  compactor,  the  Rascal  A  and  B  energy 
levels,  and  the  Rascal  C  energy  level,  respectively.   Programs  for  labor- 
atory compacted  samples  were  developed  but  are  not  included  because 
they  are  not  necessary  to  implement  the  procedures  developed  in  this 
study. 

The  statistical  methods  and  parameters  used  in  these  three 
programs  are  discussed  and  explained  in  the  "Analysis  of  Data"  chapter. 
Comment  cards  in  the  programs  also  explain  the  names  of  variables  used 
and  steps  taken  during  the  evaluation  process.   The  following  is  an 
example  of  the  programs'  use. 
GIVEN:   Statistical  prediction  models  for  the  dry  density  and  swell 

pressure  of  a  medium  plastic  clay.   These  are  developed  for  the 
Caterpillar  and  Rascal  compactors  used  in  this  study. 
DESIRED:   Having  reviewed  the  requirements  of  the  project,  the  design 

engineer  specifies  that  the  soil  is  to  be  compacted  with  8  passes 
of  the  Rascal  compactor  (E  =  2.0)  1%  wet  of  the  field  OMC,  with 
a  maximum  half -range  variability  of  water  content  of  1.6%.   The 
engineer  wishes  to  predict  the  expected  and  maximum  swell  pres- 
sures to  be  encountered  in  the  embankment. 
ANALYSIS:   The  engineer  selects  the  computer  program  for  the  Rascal  A  and 
B  energy  levels.   This  program  may  be  used  as  presented,  or  the 
initial  values  may  be  changed  to  include  only  the  values  under 
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consideration    (i.e.,    Aw   =  +1.0,    V(Aw)    =   1.6,    and   E      =   2.0). 

R 

The  magnitudes  of  dry  density  and  swell  pressure  are  evaluated 
by  the  program  using  the  appropriate  regression  equations: 


Yd  =  -  1.378  Aw  +  1.044  E   +  112.0 

2 
SP  =  0.2919  Aw  +  4.593  E   -  0.02357  Aw  y,   +    9.908  . 

R  d 

To  evaluate  the  variabilities  of  y        and   SP,  the  program  utilizes  the 
values  of  [X'X]    derived  from  the  data  sets  of  this  study.   See  the 
"Analysis  of  Data"  chapter  for  a  complete  explanation  of  the  variabil- 
ity formula  and  the  statistical  parameters  used.   V  (y    )  is  calculated 

at   Aw  ±  V(Aw)  and  the  largest  value  of   V(v,)  is  chosen.   For  swell 

d 

pressure  variability,  the  following  conditions  are  assessed: 
A  Yd  V(SP) 


Aw  +  1.6  y,    +   V(yJ  V(SP)  , 

da  1 

Aw  +  1.6  y      -    V(y,)  V(SP). 

Q.        \j.  2. 

Aw  -  1.6  vd  +  V(Y(3)  V(SP) 

Aw  -  1.6  y^    -  V(y,)  V(SP) „ 

d      d  4 

and  the  largest  value  of  V(SP)  is  selected  as  the  appropriate  measure 

of  variability. 

SOLUTION:   From  the  computer  output: 

y\   =  112.7  pcf 
d 

V(yd)  =  1.5  pcf 
SP  =  16.7  psi 


This  value  of  y        may  need  to  be  adjusted  to  give  the  proper  value 
expected  in  the  field.   See  the  "Application  of  Results". 


158 


V(&P)  =  6.7  psi 

SP    =  SP  +  V(SP)  =  23.5  psi 
max  c 

Therefore,  for  the  conditions  of  compaction  specified,  the  average 
expected  swell  pressure  will  be  16.7  psi  and  the  maximum  expected  swell 
pressure  will  be  23.5  psi  in  the  compacted  fill. 

If  these  values  of  swell  pressure  are  not  desirable  for  the 
project  under  consideration,  the  engineer  can  change  the  imput  values  of 
Aw,    V(Aw),  or  E   until  a  more  suitable  swell  pressure  is  attained. 

R 

As  written,  the  programs  automatically  increment  Aw,  V(Aw) ,  and  E  ; 

R 

however,  if  the  programs  are  used  1:0  evaluate  a  single  set  of  Aw,  V(Aw), 
and  E  ,  these  values  must  be  changed  manually  for  each  run  of  the  pro- 

R 

gram.   In  this  manner  the  engineer  can  obtain  compaction  conditions  and 

values  of  swell  pressure  that  are  compatible  with  other  requirements  of 

the  project. 

If  the  engineer  wishes  to  control  the  expected  or  maximum 

swell  pressure  to  some  predetermined  value,  he  can  obtain  the  appro- 
priate compaction  conditions  by  using  the  methods  outlined  in  the 
"Design  Engineering"  section.   Figures  5.1  thru  5.4  in  that  section 
are  generated  from  the  data  produced  by  the  computer  programs  presented 

in  this  appendix.   If  he  wishes  to  use  the  programs  to  control  the 
allowable  swell  pressure  he  must  run  the  programs  in  their  entirety, 

select  the  desired  values  of  swell  pressure  from  the  output,  and  choose 

the  conditions  of  compaction  (also  presented  in  the  output)  which  are 

most  appropriate  for  the  project  at  hand. 

The  programs  presented  also  generate  data  which  are  applicable 

to  quality  control  work;  the  "Quality  Control"  section  explains  such 

usage.   Some  additional  notes  about  the  programs  are  now  included. 
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At  OMC  (Aw  =0.0)  all  parameters  are  evaluated  twice  -  once  using 
the  dry-side  regression  equations  and  once  using  those  for  the  wet 
side. 

For  the  Caterpillar  program,  the  values  of  SP  and  V(SP)  do 
not  change  for  different  values  of  E  ,  all  other  variables  remaining 

R 

constant.   This  is  because  the  swell  pressure  is  a  function  of  Aw  only. 
The  magnitude  and  variability  of  dry  density  wet  of  OMC  do  change 
with  E  ,  and  the  output  contains  these  values  for  quality  control  work. 
These  programs  should  produce  sufficient  information  for  most 
projects.   However,  if  intermediate  values  of   Aw,  V(Aw),  or  E   are 

R 

required,  the  programs  can  be  changed  to  include  them.   This  can  be 
done  by  changing  the  initial  values,  size  of  the  increments,  or  number 
of  iterations.  It  must  be  noted  that  use  of  these  parameters  outside 
the  ranges  given  in  each  program  constitutes  an  extrapolation,  and 
must  be  usad  with  caution. 
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COMPUTES  FROGRAM  FCR  EVALUATING  SWELL  PRESSURE  UflRIfiBILITY,  CATERPILLAR 

C 

C  WCO=WATER  CONTENT  DEUIATICN  FROM  CMC 

C  UUCO=EXPECTED  WATER  CONTENT  UARIA3ILITY  (HALF  RANGE) 

C  AD=(X'X)_1FOR  DRY  DENSITY,  DRY  OF  OMC 

C 

C  AW=CX'X)"1  FGR  DRY  DENSITY,  WET  OF  OMC 

C  XDD=UECTOR  OF  INDEPENDENT  UARIA3LE5  fCP.   DRY  DENSITY,  DRY  OF  OM^ 

C  XDW=UECTOR  OF  INDEPENDENT  UARIABLE5  FGR  DRY  DENSITY  WET  OF  OMC 

C  TDD=T-STATISTIC  FOR  DRY  DENSITY,  DRY  OF  CMC  l-'31'Y*  Ut'    UnC 

C  TDW=T-STATISTIC  FOR  DRY  DENSITY,  WET  OF  OMC 

C  MSEDD=MEAN  SQUARE  ERROR  FOR  DRY  DENSITY,  DRY  GF  CMC 

C  MSEDW=MEAN  SQUARE  ERROR  FCR  DRY  DENSITY,  WET  OF  OMC 

C  UDD=EXPECTED  DRY  DENSITY  UARIABILITY 

C  , 

C  BD=CX'X)_1  FOR  SWELL  PRESSURE,  DRY  OF  OMC 

C         _! 

C  BW=CX'X)    FOR  SWELL  PRESSURE,  WET  CF  OMC 

C  XSD=UECTCR  OF  INDEPENDENT  UARIA3LE3  FOR  3P,  DRY  CF  OMC 

C  XSW=UECTOR  CF  INDEPENDENT  UARIA3LZ5  FCR  SP,  WET  OF  OMr 

C  TSD=T-STATISTIC  FOR  SWELL  PRESSURE,  DRY  GF  CMC 

C  TSW=T-STATISTIC  FOR  SWELL  PRESSURE,  WET  GF  OMC 

C  MSESD=MEAN  SCUARE  ERROR  FCR  SWELL  PRESSURE,  DRY  OF  OMC 

C  MSESW=MEAN  SQUARE  ERROR  FGR  SWELL  PRESSURE,  WET  OF  OMC 

C  USP=EXPECTED  SWELL  PRESSURE  UARIABILITY 

C 

REAL  MSEDD, MSEDW. MSESD, MSESW 

DIMENSION   XDDC2),XDWC3),AD(2,2),AWC3.3),XSDC2), 
SXSUC3),BD(2,2),3WC3,3) 

1  FORMAT(SX,5F10.7/F10.4) 

2  FCRMAT(9X,GF10.7,F10.4/"F10.4) 

3  F0RMATO//2(24X,2F12.7W/24X,2F12.4) 

4  FORMAT <V//3(24X,3F12.?W/24X,2F12. 4"///- ) 

5  FORMATC/X.10F3.1) 
C 

C  READ  DATA 

C 

READ   1. ADC  1,1), ADC  1,2), AD C2, 1 ) , ADC  2, 2) , MSEDD, TDD 

READ  2,AWC1,1),AWC1,2),AWC1,3),AWC2,2),AWC2,3),AWC3,3),MSEDW,TDW 

READ   1 , ED ( 1 , 1 ) , 3D ( 1 , 2 ) , 3DC  2, 1 ) , BDC  2 , 2 ) , MSESD, TSD 

READ  2,EWC1,1),EWC1,2),EWC1,3).3WC2,2),3WC2,3),BWC3,3),MSESW,TSW 

C  ASSEMBLE  SYMMETRIC  MATRICIES 


C 


AWC2,1)=AWC1,2) 
AWC3,1)=AWC1,3) 
AU(3,2)=AW(2,3) 
BWC2,1)=BWC1.2) 
BW(3,1)=BWC1,3) 
BWC3,2)=BWC2,3) 
PRINT  3, AD, MSEDD, TDD 
PRINT  4,AW,MSEDW,TDW 
PRINT  3, 3D, MSESD, TSD 
PRINT  4,BW,MSESW,TSW 
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C 

C  SET  INITIAL  UALUE3 

C 

XDD(1)=1.0 

XDW(1)=1.0 

XSD(1)=1.0 

XSUC1)=1.0 

DO  30  1=1.3 

IF  (I  .EQ.  1)  ER=1.0 

IF  (I  .EQ.  2)  ER=2.0 

IF  (I  .EQ.  3)  ER=4.0 

PASS=4*ER 

UUCO=0.0 

DO  20  J=l,10 

UUCO=UUCO+0.2 

UC0=-2.5 

DO  10  N=l,27 

UCO=UCO+0.5 
C 

C  CALCULATE  DRY  DENSITY  UARIA3ILITY,  DRY  CF  CMC 
C 

IF  (WCO  .GT.  0.0)  GO  TO  100 

ADD1=ADD2=0.0 

XDD(2)=UC0-UUICO 

ADD1=XD3U)*ADU.1)+XDD(2)«ADC2. 15 

ADD2=XDD(1)*AD(1,2)+XDD(2)*ADC2,2) 

SDDS01=MSEDD*(ADD1»XDD( 1 )+ADD2*XDDC2) ) 

SDD1=SQRT(SDD501) 

UDD1=SDD1*TDD 

ADD1=ADD2=0.0 

XDD(2)=XDD(2)+(2.0*UWCO) 

ADD1=XDDC1)*AD(1,1)+XDD(2)*A0(2,1) 

ADD2=XDD ( 1 ) -AD (1,2) +XDD ( 2 ) * AD  ( 2 .  2 ) 

SDDS02=MSEDD->  ( ADD  1  *XDD  ( 1 )  +  ADD2»XDD  ( 2 ) ) 

SDD2=SORT(SDDSQ2) 

UDD2=SDD2»TDD 

UDD=l»DD2 

IF  (UDD1  .GT.  UDD2)  UDD=UDD1 
C 

C  EUALUATE  REGRESSION  ECUATION  FOR  DRY  DENSITY,  DRY  OF  CMC 
C 

DD=5.3S3*UC0+113.3 

DDF=0.33»DD 
C 

C  CALCULATE  SWELL  PRESSURE  UARIA3ILITY,  DRY  OF  OMC 
C 

BSB1=BSD2=0.0 

XSD(2)=UC0-UUC0 

BSD1=XSD(1)*BD(1,1)+XSD(2)*30(2,1) 

BSD2=XSD(I)*3D(1,2)+XSDC2)*3DC2,2) 

SSDS01=NSES0*(BSD1*XSB( 1 )+BSD2-XSD(2) ) 

SSD1=SQRT(SSDSQ1) 

USD1=SSD1*TSD 

BSD1=BSD2=0.0 

XSD(2)=XSO(2)  +  (2.0*UU!CO) 

BSD1=XSD(1)*BD(1,1)+XSD(2)*3D(2,1) 

BSD2=XSD(1)*30(1,2)+XSD(2)*3DC2,2) 

SSDS02=riSE5D*(ESDl*XSD(l)+3SD2*X5D(2)) 

SSD2=S0RT(S3DSQ2) 

USD2=SSD2*TSD 

USP=USD2 

IF  (USD1  .GT.  USD2)  US?=USD1 
C 

C  EUALUATE  REGRESSION  EOLATION  FOR  SUIELL  ?RZSS'C2E,    DRY  OF  CMC 
C 

SP=21.72*UC0+44.12 

SPMAX=SP+USP 

GO  TO  200 
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C 

C  CALCULATE  DRY  DENSITY  UARIABILITY,  UET  OF  OKC 

C 

100  IF  (N  .Ed.  S)  WCO=0.0 

ADU1=ADW2=ADU3=0.0 

XDU(25=UCO-UUCO 

XDWC3)=ER 

A0Wl=XDU(l)»AU(l»n+XI)U(2)»AlJC2,l)+XBW(3)*AUC3il) 

ADU2=XDU(1)*AW(1,2)+XDU(2:>*AW(2,2)+XDW(3)*AW(3>2) 

ADW3=XDWm«AW(1.3)+XDW(2)*AU(2,3:i+XDW(3)*AW(3,3) 

SDUS01=MSEDW»(ADW1»XDU(1)+ADW2*XDU(2:>+ADW3*XDW(3:>:) 

SDW1=SQRT(SDWSQ1) 

UDU1=SDU1»TDW 

ADW1=ADW2=ADW3=0.0 

XDUC2)=XDW(2)+(2.0*UWCO) 

ADWl=XDW(l)«AW(l.l)+XDW(23*AW(2.1)+XDU(3)»AW(3f 1) 

ADW2=XDUl(l)*AU(1.2)+XDU(2)*AW(2»2)+XEUC3)*ftU(3.2) 

ADW3=XDU(1)*AW(1,3:>+XDW(2)*AW(2,3)+XEW(3)*AW(3,3) 

SDWSQ2=MSEDW*(ADUIl*XmJ(l)+ADW2»XDW(2)+ADU3*XDUC3)) 

SDW2=SQRT(SDWSQ2) 

UDW2=SDW2»TDU 

UDD=UDU2 

IF  CUDU1  .GT.  UDW2)  UDD=UDW1 
C 

C  EVALUATE  REGRESSION  EQUATION  FOR  DRY  DENSITY,  UET  OF  ONC 
C 

DD=-1.G09»UCO+1.00S«ER+114.3 

DDF=0.93»DD 
C 

C  CALCULATE  SWELL  PRESSURE  UARIABILITY,  WET  OF  CMC 
C 

BSW1=BSW2=BSU3=0.0 

XSW(2)=WC0-UWC0 

XSW(3)=XSU(2)*»2 

BSW1=XSU(1)»BUU.1)+XSUK2)*3UC2.1)+XSUK3)«3W(3. 1) 

BSW2=XSW(1)*BW(1.2:+XSW(2)«BW(2,2)+XSW(3)*BW(3,2) 

3SU3=XSW(1)»BW(1,3)+XSW(2)*BW(2,3)+XSU(3)*3WC3,3) 

SSWS01=MSESW*(BSWl*XSW(l)+BSW2*XSW(2)+BSW3»XSUC3n 

SSU1=SQRT(SSWSQ1) 

USW1=SSW1*TSU 

BSW1=BSW2=BSU3=0.0 

XSW(2)=XSW(2)+(2.0*UWCO) 

XSU(3)=XSW(2)**2 

BSUl=XSW(l)»BW(l,n+XSW(2)*BW(2.  n+XSWC3)*3U(3.  1) 

BSW2=XSW(1)*BW(1,2)+XSW(2)*BW(2.2)+XSU(3)*3W(3.2) 

BSW3=XSU(l)»BUIClt3)+XSW(2)*BW(2,3)+XSW(3)*BW(3,3) 

SSWSQ2=MSESW*(BSUl»XSW(l)+B5U2*XSW(2^+3SW3*XSW(3n 

SSW2=SQRT(SSUSQ2) 

USW2=SSW2*TSU 

USP=USW2 

IF  (USW1  .GT.  USW2)  USP=USW1 
C 

C  EVALUATE  REGRESSION  EQUATION  FOR  SWELL  PRESSURE,  WET  OF  CMC 
C 


WC02=WC0*WC0 

SP=-1 0 . 87*WCO+0 . 7S93*WC02+43 . 03 

SPt1AX=SP+USP 
200  PRINT  5,UWC0.WC0,PASS,ER,DD,DDFiUDD.SP,USP,SPi1AX 
10  CONTINUE 
20  CONTINUE 
30  CONTINUE 

END 
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C  THE  PRINTED  DATA  MATRICIES  FOLLOU 


.1784300 
.1967400 


.1967400 
.4070600 


9.9970 


2.2280 


.1555000 
.0109300 
•038SS00 


.0109300 
. 00SSE80 
.0030630 


-.0386500 

-.0030G80 

•0192100 


2.7190 


2.0420 


.1784300 
.1967400 


.1967400 
.4070600 


9.9970 


2.2280 


.0837600 
.0385100 
.0033470 


.0385100 
.0311400 
.0033620 


,0033470 

-.0033620 

.0004353 


74.0580 


2.0340 
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COMPUTER  PROGRAM  FCR  EUALUATING  SUELL  PRESSURE  UARIABILITY,  RASCAL  A+S 

C 

C 

C  UCO=UATER  CONTENT  DEUIATION  FROM  OMC 

C  UUCO=EXPECTED  WATER  CONTENT  UARIABILITY  (HALF  RANGE) 

C  _n 

C  AD=CX'X)  l  FOR  DRY  DENSITY,  DRY  OF  OMC 

C  -1 

C  AU=(X'X)    FOR  DRY  DENSITY,  WET  OF  OMC 

C  XDD=UECTOR  OF  INDEPENDENT  UARIABLES  FOR  DRY  DENSITY,  DRY  OF  OMC 

C  XDW=UECTOR  OF  INDEPENDENT  UARIASLES  FOR  DRY  DENSITY,  WET  OF  OMC 

C  TDD=T-STATISTIC  FOR  DRY  DENSITY.  DRY  OF  OMC 

C  TDU=T-STATISTIC  FOR  DRY  DENSITY,  WET  OF  OMC 

C  MSEDD=MEAN  SQUARE  ERROR  FOR  DRY  DENSITY,  DRY  OF  OMC 

C  MSEDU=MEAN  SQUARE  ERROR  FOR  DRY  DENSITY.  WET  OF  OMC 

C  UDD=EXPECTED  DRY  DENSITY  UARIABILITY 

C        -1 

C  B=(X'X)  l    FOR  SUELL  PRESSURE 

C  XS=UECTOR  OF  INDEPENDENT  UARIABLES  FOR  SWELL  PRESSURE 

C  TSP=T-STATISTIC  FOR  SWELL  FRESSURE 

C  MSESP=MEAN  SQUARE  ERROR  FCR  SUELL  PRESSURE 

C  USP=EXPECTED  SWELL  PRESSURE  UARIABILITY 

C 

REAL  MSEDD.MSEDU.MSESP 

DIMENSION   XDD(3).XDU(3),AD(3.3),AU(3.3),XS(4),3(4,4) 

1  FORMAT(9X,GF10.7,F10.4/F10.4) 

2  FORMAT(9X,4F15.11/5F15.11/'F15.11,F10.4/'F10.4) 

3  F0RMAT(/v/3(24X.3F12.7/V/24X,£F12.4) 

4  FORMAT(///4(24X,4F15.10/)//24X,2F15.4/////) 

5  FORMATC/X.lOFS.n 
C 

C  READ  DATA 
C 

READ   l,AD(l»n.ADa,2).AD(l,3),AD(2,2),AD(2,3),AD(3,3),MSEI)D,TDD 

READ   l,AW(l.l).Ag(l,2),AU(l,3),AWC2,2),AU(2,3),AU(3,3),MSEDU,TDU 

READ  2,Ba,l).3(1.2).B(1.3).B(1.4),3(2,2),3(2,3).B(2,4).BC3.3). 
$B(3.4),B(4,4).MSESP.TSP 
C 

C  ASSEMBLE  SYMMETRIC  MATRICIES 
C 

AD(2,n=AD(l,2) 

AD(3.1)=AD(1,3) 

AD(3.2)=AD(2,3) 

AW(2,1)=AW(1,2) 

AU(3.1)=AU(1,3) 

AU(3.2)=AW(2,3) 

B(2,1)=B(1,2) 

B(3.1)=B(1.3) 

B(3.2)=B(2,3) 

B(4.1)=B(1.4) 

B(4,2)=B(2.4) 

B(4.3)=B(3.4) 

PRINT  3,AD,MSEDD,TDI3 

PRINT   3,AU,MSEDW,TDU 

PRINT  4,B,MSESP,TSP 
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C 

C  SET  INITIAL  UALUES 

C 

XTJDC1)=1.0 

XDU(1)=1.0 

XS(1)=1.0 

DO  30  1=1, a 

IF  (I  .EQ.  1)  ER=1.0 

IF  (I  .EQ.  2)  ER=2.0 

PASS=4*ER 

UUCO=0.0 

DO  20  J=l,10 

UUCO=UUICO+0.2 

WCO=-3.0 

DO  10  N=l,21 

UCO=UCO+0.5 
C 

C  CALCULATE  DRY  DENSITY  UARIABILITY,  DRY  OF  OMC 
C 

IF  CWCO  .GT.  0.0)  GO  TO  100 

ADD1=ADD2=ADD3=0.0 

XDD(2)=UIC0-UUIC0 

XDDC3)=ER 

ADDl=XDD(l)*AD(l,l)+XDDC2)*AD(2»l)+XDDi:3)»ADi:3.1) 

ADD2=XDD(1)*AD(1,2)+XDD(2)*AD(2.2)+XDD(3)»AD(3.2) 

ADD3=XDD(l)*AD(1.3)+XDD(2)»AD(2»3)+XDD(3)*AD(3i3) 

SDDSGl=t1SEDD»(ADDl*XDDCl)+ADD2*XDD(2)+ADD3*XDD(3)) 

SDDl=SGRTi:SDDSOn 

UDD1=SDD1»TDD 

ADD1=ADD2=ADE3=0.0 

XDD(2)=XDDC2)+(2.0»UUCO) 

ADD1=XDDC1)»AD(1, 1 )+XDD(2)*AD(2. 1 )+XDDC3)*AD(3, 1 ) 

ADD2=XDD(1)«AD(1.2)+XDD(2)*AD(2,2)+XDD(3)*AD(3,2) 

ADD3=XDDC1)«AD(1,3)+XDD(2)*ADC2,3)+XDD(3)*AD(3.3) 

SDDS02=MSEDD* ( ADD1*XDD ( 1 ) +ADD2*XDD ( 2 ) +ADD3»XDD  C  3 ) ) 

SDD2=SQRT(SDDSQ2) 

UDD2=SDD2»TDD 

UDD=UDD2 

IF  (UDD1  .GT.  UDD2)  UDD=UDD1 
C 

C  EUALUATE  REGRESSION  EQUATION  FOR  DRY  DENSITY,  DRY  OF  ONC 
C 

DD=1 . 90S»WC0+1 . 330*ER+1 12. 5 

DDF=0.33»DD 

GO  TO  200 
C 

C  CALCULATE  DRY  DENSITY  UARIABILITY,  WET  OF  OMC 
C 


100  IF  (N  .EQ.  7)   UCO=0.0 
ADUI1=ADW2=ADU3=0.0 
XDUC2)=UC0-UUC0 
XDW(3)=ER 

ADWl=XDWCl)«AW(l,l)+XDUI(2)»AUC2,l)+XDWC3)*Ag(3pl) 
ADW2=XDW(1)*AW(1,2)+XDU(2)»AU(2,2)+XDU(3)*AW(3.2) 
ADU3=XDU(1)«AU(1,3)+XDUK2)«AU(2,3)+XDUC3)*AU(3,3) 
SDUSQ 1 =MSEDU» < ADU 1 »XDU ( 1 ) + ADW2»XDU ( 2 ) + ADU3»XDU ( 3 )) 
SDW1=SQRT(SDUISQ1) 
UDU1=SDU1»TDW 
ADUI1=ADW2=ADW3=0.0 
XDW(2)=XDU(2)+(2.0*UUC0) 

ADUl=XDW(l)*AW(l,n+XDW(2)*AUI(2,l)+XDW(3)*AUI(3,n 
ADW2=XDU(1)«AU(1,2)+XDW(2)*AU(2,2)+XDUC3)*AW(3,2) 
ADU3=XDUI(1)*AW(1,3)+XDU(2)*AU(2,3)+XDW(3)*AU(3,3) 
SDUSQ2=MSEDU* ( ADW 1 *XDW ( 1 ) +ADU2»XDW  C  2 ) +ADW3»^DW ( 3 ) ) 
SDW2=SQRT(SDUSQ2) 
UDUI2=SDU2»TDW 
UDD=UDW2 
IF  (UDU1  .GT.  UDW2)  UDD=UDW1 
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C 

C  EVALUATE  REGRESSION  EQUATION  FOR  DRY  DENSITY.  WET  OF  GMC 


C 


DD=-1 . 373*WC0+1 . 044*ER+112. 0 
DDF=0.93*DD 


C 

C  EUALUATE  REGRESSION  ECUATION  FOR  SWELL  PRESSURE 

C 

200  UD=UCO*DD 

UC02=UC0*UC0 

SP=0 . 29 19*UC02+4 . 593*ER-0 . 02357»UD+9 .903 

C 

C  CALCULATE  SHELL  PRESSURE  UARIABILITY 

C 

Z1=Z2=Z3=Z4=0.0 
XS(2)=(UC0+UUC0)*»2 
XS(3)=(UC0+UUC0)*(DD+UDD) 
XS(4)=ER 

Z1=XS(1)«B(1»1)+XS(2)«3(2,1)+XSC3)*B(3.1)+XSC4)*B(4, 1) 
Z2=XSU)*B(1.2)+XS(2)»3(2,2)+XS(3)*B(3,2)+XS(4)»BC4,2) 
Z3=XS(1)»3(1,3)+XS(2)*3(2»3)+XS(3)*3(3,3)+XS(4)«B(4.33 
Z4=XS(1)*3C1,4)+XS(2)*B(2,4)+XSC3)*3(3.4)+XS(4)*BC4,4) 
SSPSQ1=MSESP»(Z1*XSC1)+Z2*XS(2)+Z3»XSC3)+Z4*XS(4» 
SSP1=SQRTCSSPSQ1) 
U1=SSP1»TSP 
Z1=Z2=Z3=Z4=0.0 
XS(2)=(UC0+UUC0)»»2 
XS  C  3 )  =  ( UCO+UUICO )  *  ( DD-UDD ) 

Zl=XS(l)*B(l,l)+XS(2)»3(2,l)+XS(3)*3(3.1)+XS(4)*3(4,n 
Z2=XSC1)»B(1.2)+XSC2)«3(2.2)+XSC3)»BC3.2)+XS(4)*B(4.2) 
Z3=XS(1)»B(1»3)+XS(2)»B(2,3)+XS(3)»3(3,3)+XS(4)*B(4,3) 
Z4=XS(1)«3(1.4)+XS(2)»B(2,4)+XSC3)*B(3.4)+XS(4)*3C4.4) 
SSPS02=MSESP*(Zl*XS(n+Z2*XS(2)+Z3*XS(3)+Z4*XS(4)) 
SSP2=SQRT(SSPSQ2) 
U2=SSP2»TSP 
Z1=Z2=Z3=Z4=0.0 
XS(2)=(UC0-UUC0)**2 
XS ( 3 ) = ( UCO-UUCO ) » ( DD+UDD ) 

Z1=XS(1)*B(1.1)+XS(2)»B(2»1)+XS(3)*B(3.1)+XSC4)*B(4. 1) 
Z2=XS(1)*B(1,2)+XS(2)»B(2.2)+XSC3)*B(3,2)+XS(4)»3(4,2) 
Z3=XS(1)*B(1.3)+XS(2)»B(2»3)+XS(3)»B(3»3)+XSC4)*3(4,3; 
Z4=XS(1)»B(1,4)+XS(2)»B(2,4)+XSC3)*3(3.4)+XS(4)*BC4,4) 
SSPSQ3=MSESP» ( Zl *XS ( 1 ) +Z2»XS ( 2 ) +Z3»XS ( 3 ) +Z4*XS ( 4 ) ) 
SSP3=SQRT(SSPSQ3) 
U3=SSP3»TSP 
Z1=Z2=Z3=Z4=0.0 
XS(2)=(UCQ-UUC0)»*2 
XS ( 3 ) = ( UCO-UUCO ) * ( DD-UDD ) 

Zl=XS(n»BCl.l)+XS(2)»B(2,l)+XS(3)*B(3.1)+XS(4)*3(4,l) 
Z2=XS(1)*B(1.2)+XS(2)»B(2,2)+XS(3)*3(3,2)+XS(4)*B(4,2) 
Z3=XS(l)*B(lt3)+XS(2)*B(2,3)+XS(3)*3C3,3)+XSC4)»BC4,3) 
Z4=XS(1)«BC1,4)+XS(2)»B(2.4)+XS(3)»3(3,4)+XS(4)»B(4,4) 
SSPSQ4=MSESP»(Z1»XS(1)+Z2*XS(2)+Z3»XSC3)+Z4»XS(4)) 
SSP4=S0RT(SSPSQ4) 
U4=SSP4»TSP 
USP=U4 

IFCV1.GE.U2.AND.U1.GE.U3.AND.U1.GE.U4)  USP=U1 
IFCU2.GE.U1.AND.U2.GE.U3.AND.U2.GE.U4)  USP=U2 
IFCU3.GE.U1.AND.U3.GE.U2.AND.U3.GE.U4)  U3P=U3 
SPMAX=SP+USP 

PRINT  5, UWCO. UCO. PASS. ER. DD, DDF, UDD, SP, USP. SPMAX 
10  CONTINUE 
20  CONTINUE 
30  CONTINUE 
END 
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C 

C 

c 


THE  PRINTED  DATA  MATRTCIES  FOLLOW 


.4553000 

.2148000 

-.0702000 


.2148000 
.2167000 
.0126500 


-.0702000 
.0126500 
.0385600 


2.3610 


2.1600 


.1634000 
-.0141700 
-.0416300 


.0141700 
.0043280 
.0001782 


-.0416300 

-.0001762 

.0173400 


5.4380 


2.0400 


.3136000000 
-.0048720000 

.0001175000 
-.1871000000 


.0048720000 
.0003815000 
.0000388300 
,0013670000 


.0001175000 
-.0000388300 

.0000013580 
-.0000740000 


-.1871000000 
.0013670000 

-.0000740000 
.1280000000 


55.4250 


2.G830 
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COMPUTER  PROGRAM  FOR  EVALUATING  SWELL  PRESSURE  USRIA3ILITY,  RASCAL  C 

I— 

C 

C  UC0=WATER  CONTENT  DEVIATION  FROM  OMC 

C  UWCO=EXPECTED  WATER  CONTENT  UARIA3ILITY  (HALF  RANGE) 

C         -1 

C  AD=(X'X)  x  FOR  DRY  DENSITY,  DRY  OF  OMC 

C  , 

C  AW=(X'X)  ■'■FOR  DRY  DENSITY,  WET  OF  OMC 

C  XDD=UECTOR  OF  INDEPENDENT  UARIABLES  FOR  DRY  DENSITY,  DRY  OF  OMC 

C  XDW=UECTOR  OF  INDEPENDENT  UARIA3LES  FOR  DRY  DENSITY,  WET  OF  OMC 

C  TDD=T-STATISTIC  FOR  DRY  DENSITY,  DRY  OF  CMC 

C  TDW=T-STATISTIC  FOR  DRY  DENSITY,  WET  OF  OMC 

C  MSEDD=MEAN  SQUARE  ERROR  FOR  DRY  DENSITY,  DRY  OF  OMC 

C  MSEDW=MEAN  SQUARE  ERROR  FOR  DRY  DENSITY,  WET  OF  OMC 

C  UDD=EXPECTED  DRY  DEMSITY  UARIABILITY 

C        -1 

C  B=(X'X)  x  FOR  SWELL  PRESSURE 

C  XS=UECTOR  OF  INDEPENDENT  UARIABLES  FOR  SWELL  PRESSURE 

C  TSP=T-STATISTIC  FOR  SWELL  PRESSURE 

C  MSESP=MEAN  SQUARE  ERROR  FOR  SWELL  PRESSURE 

C  USP=EXPECTED  SWELL  PRESSURE  UARIABILITY 

C 

REAL  MSEDD.MSEDW,f1SESP 

DIMENSION  XDD(3).XDU(3),AD(3.3),AW(3,3),XSC3),B(3.3) 

1  FORMAT(9X.6F10.7,F10.4/F10.4) 

2  FORMAT(9X,GF10.7,F10.4/F10.4) 

3  F0RMAT(///3(24X,3F12.?/')//24X,2F12.4) 

4  F0RMATC/O'3(24X.3F12.7/)"24X.2F12.4/VV) 

5  FORMATC/X.10F8.1) 
C 

C  READ  DATA 
C 

READ  1,AD(1,1),AD(1,2),ADC1,3),AD(2,2),AD(2,3),AD(3,3),MSEDD,TDD 

READ   1,AW(1,1),AW(1,2),AWU,3),AW(2,2),AW(2.3),AWC3,3),MSEDW,TDW 

READ  2,B(1,1),B(1,2),B(1,3),B(2,2),3C2,3),B(3.3),MSESP,TSP 
C 

C  ASSEMBLE  SYMMETRIC  MATRICIES 
C 

ADC2, 1)=AD(1,2) 

AD(3.1)=AD(1,3) 

AD(3,2)=AD(2,3) 

AW(2,1)=AW(1,2) 

AW(3,l)=flW(l,3) 

AW(3,2)=AW(2,3) 

B(2,1)=B(1,2) 

B(3,1)=B(1,3) 

B(3f2)=BC2i3) 

PRINT  3,AD,MSEDD,TTJD 

PRINT  3,AW,MSEDW,TDU 

PRINT  4,B,MSESP,TSP 
C 

C  SET  INITIAL  UALUES 
C 

XDD(1)=1.0 

XDW(1)=1.0 

XS(1)=1.0 

ER=4.0 

PASS=4»ER 

UWCO=0.0 

DO  20  J=l,10 

UWCO=UWCO+0.2 

WC0=-2.5 

DO  10  N=l,25 

WCO=WCO+0.5 
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C 

C     CALCULATE  DRV  DENSITY  LIABILITY,    DRY  OF  OMC 

IF   (UCO   .GT.    0.0)   CO  TO  100 
ADD1=ADD2=ADD3=0.0 
XDD(2)=UC0-UUC0 
XDD(3)=ER 

^°°i=^°(1)#f:iD(l'^+XDDC2)*ADC2,n+XDD(3)*AD(3,l) 
ADD2=XDD(n*AD(l,2)+XDDC2)»AD(2,2)+XDDC3   «AD   3   P 
ADD3=XDD(1)*AD(1,3KXDD(2)*AD(2  3  +XDD  3  ->AD  3  3 

UDD1=SDD1*TDD 

ADD1=ADD2=ADD3=0.0 

XDD(2)=XDD(2)+(2.0»UUCO) 

ADD1=XDD(1)»AD(1, 1 )+XDD(2)*AD(2, D+XDD(3)«ADf3. 1 1 

«D°2=XDDCl).flDCl,a)+XDDC2).AD<2,2  +XM  3  »2S  3  2 

ADD3=XDDC1)«ADU.3)+XDD<2)*AD(2,3)+XDDC3   *AD   3   3 

UDD2=SDD2*TDD 
UDD=UDD2 

IF  (UDD1  .GT.  UDD2)  UDD=UDD1 
C 

C  EUALUATE  REGRESSION  EQUATION  FOR  DRY  DENSITY,  DRY  OF  OMC 


C 


c 


DD=1.90G*LICO+1.390»ER+112.5 

DDF=0.93*DD 

GO  TO  200 


C  CALCULATE  DRY  DENSITY  UARIABILITY,  WET  OF  OMC 

100  IF  (N  .EQ.  6)  UCO=0.0 
ADU1=ADU2=ADU3=0.0 
XDU(2)=UC0-UUC0 
XDUC3)=ER 

ADUl=XDWCn»AU(l,l)+XDUC2)*AW(2,l)+XDUI(3)«AW(3,n 

ADU2=XDU(l)«AU(l,2)+XDU(2)»AUIC2.0)+XDU(3)eAU(3,2) 

ADW3=XDU(1)»AW(1,3)+XDU(2)«AU(2,3)+XDW(3)*AW(3,3) 

SDWSQ1=MSEDUI*(ADU1»XDW(1)+ADW2»XDU(2)+ADW3»XDU(3)) 
SDW1=SQRT(SDUSQ1) 

UDU1=SDW1»TDU 
ADW1 =ADU2=ADU3=0 . 0 
XDU(2)=XDU(2)+(2.0*UWCO) 

ADWl=XDW(l)*AW(l,l)+XDW(2)»AU:2,l)+XDUI(3)*AW(3,n 
ADU2=XDU(1)«AW(1,2)+XDWC2)»AW(2,2)+XDUI(3)»AW(3,2) 
ADU3=XDUCl)»Ag(l,3)+XDU(2)*AU(253)+XDU(3)*AUC3,3) 
SDUSQ2=MSEDUI» ( ADW1 »XDU ( 1 ) +ADU2»XDU ( 2 ) +ADW3*XDW C  3  » 
SDU2=SQRT(SDUSQ2) 
UDU2=SDW2»TDU 
UDD=UEUI2 

IF  (UDW1  .GT.  UDU2)  UDD=UDW1 
C 

0  EUALUATE  REGRESSION  EQUATION  FOR  DRY  DENSITY,  UET  OF  OMC 

DD=-1 .  378*UIC0+1 .  044»ER+1 12.0 
DDF=0.93«DD 
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C 

C  EVALUATE  REGRESSION  EQUATION  FOR  SWELL  FRESSURE 

C 

200  UD=UCO»DD 

SP=-1S.03»UCO+0.1501»UD+14.50 
C 

C  CALCULATE  SWELL  PRESSURE  UARIABILITY 
C 

Z1=Z2=Z3=0.0 

XS(2)=UC0+UUC0 

XS  C  3 )  =  C UCO+UUCO ) * ( DD+UDD ) 

Zl=XSU)»B(l.l)+XSC2)»B(2f  1)+XS(3)»B(3. 1) 

Z2=XSCn*B(1.2)+XSC2)*BC2,2)+XS(3)*BC3,2) 

Z3=XS(n*3(1.3)+XS(2)*B(2»3)+XSC3)»B(3»3) 

SSPSQ1=MSESP*(Z1*XSC1)+Z2»XS(2)+Z3*XSC3)) 

SSP1=SQRT(SSPSQ1) 

U1=SSP1»TSP 

Z1=Z2=Z3=0.0 

XS(2)=UC0+UWC0 

XS ( 3 )  =  ( WCO+UUCO  3  * ( DD-UDD  ) 

Zl=XS(i:>*BCl.n+XS(2)*B(2.1)+XS(3)*3(3,l) 

Z2=XSC1)*3(1,2)+XS(2)*B(2,2)+XS(3)«3(3,2) 

Z3=XS(1)*B(1,3)+XS(2)*B(2,3)+XS(3)»3(3.3) 

SSPS02=MSESP»(Zl«XSa)+Z2»XS(2)+Z3*X5(3» 

SSP2=S0RT(SSPSQ2) 

U2=SSP2*TSP 

Z1=Z2=Z3=0.0 

XS(2)=WC0-UWC0 

XS  ( 3 )  =  ( UCO-UUICO )  •  ( DD+UDD ) 

Z1=XS(1)*B(1,1)+XS(2)«B(2,1)+XS(3)*3(3,1) 

Z2=XS(1)*B(1.2)+XS(2)*B(2,2)+XS(3)»B(3.2) 

Z3=XS(1)»B(1.3)+XSC2)«B(2.3)+XSC3)*BC3.3) 

SSPSa3=MSESP»(Zl*XS(l)+Z2*XS(2)+Z3*X5(3)) 

SSP3=SQRT(SSPSQ3) 

U3=SSP3*TSP 

Z1=Z2=Z3=0.0 

XS(2)=WC0-UWC0 

XS ( 3 )  =  ( UCO-UUCO ) * ( DD-UDD  ) 

Zl=XS(l)»BCl.l)+XSC2)«B(2>n+XS(3)»BC3.1) 

Z2=XS(1)»3C1,2)+XS(2)*B(2»2)+XS(3)»3(3.2) 

Z3=XS(1)»BC1.3J+XS(2)*B(2.3)+XSC3)»B(3.3) 

SSPS04=f1SESP*(Zl»XS(l)+Z2»XS(2)+Z3*XS(3n 

SSP4=SQRT(SSPSQ4) 

U4=SSP4*TSP 

USP=U4 

IFCU1.GE.U2.AND.U1.GE.U3.AND.U1.GE.U4)  USP=U1 

IFCU2.GE.U1.AND.U2.GE.U3.AND.U2.GE.U4)  USP=U2 

IFCU3.GE.U1.AND.U3.GE.U2.AND.U3.GE.U4)  USP=U3 

SPMAX=SP+USP 

PRINT  5. UUCO. WCO, PASS. ER, DD. DDF. UDD. SP. USP. SPttAX 
10  CONTINUE 
20  CONTINUE 

END 
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C  THE  PRINTED  DATA  MATRICIES  FOLLOW 


.45B3000 

.2148000 

-.0702000 


.2148000 
.2157000 
.0125500 


-.0702000 
.012S500 
.0385500 


2.3810 


2.1800 


.1854000 
.0141700 
.0418300 


.0141700 
.0043280 
,0001782 


.0416300 
,0001782 
.0173400 


5.4380 


2.0400 


.0848000    .0131500   -.0002302 

.0131500   2.4035000   -.0228100 

-.0002302   -.0228100    .00021G5 


24.3200 


2.5580 


u 
DC 
o 


z 
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LU 
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